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Preface 


The  Earthquake  Engineering  and  Geosciences  Division  (GG),  Geotechnical 
Laboratory  (GL),  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES), 
conducted  the  workshop  described  in  this  report  during  the  period  27-28  June ' 
1995.  Funding  for  the  workshop  was  provided  by  the  Earthquake  Engineering 
Research  Program  (EQEN),  Work  Unit  No.  33014,  “Geophysical  Methods  for 
Site  Characterization.” 

Mr.  Robert  F.  Ballard,  Jr.,  is  the  Principal  Investigator  for  EQEN  Work  Unit 
No.  33014.  Dr.  Joseph  P.  Koester,  Earthquake  Engineering  and  Seismology 
Branch  (GG-H),  and  Ms.  Tina  L.  Holmes,  Soil  Mechanics  Branch,  Soil  and  Rock 
Mechanics  Division  (GS-S),  compiled  the  proceedings  of  the  subject  workshop. 
The  EQEN  Program  Manager  is  Dr.  Mary  E.  Hynes,  GG-H.  EQEN  research  was 
conducted  under  the  general  supervision  of  Dr.  Arley  G.  Franklin  Chief,  GG,  and 
Dr.  William  F.  Marcuson  III,  Director,  GL. 

This  report  is  a  summary  of  presentations  and  discussions  that  took  place 
during  the  subject  workshop,  as  transcribed  by  Dr.  Koester  and  Ms.  Holmes.  Not 
all  presenters  submitted  proceedings  papers;  hard  copies  of  presentation  materials 
and  available  papers  are  provided  in  appendices  for  those  materials  that  were 
available. 

Director  of  WES  during  the  conduct  of  this  study  and  preparation  of  the 
report  was  Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K  Howard 
EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes .  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Multiply 

By 
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kilograms  force  per  square  centimeter 
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pascals 

kilopounds  per  square  foot  (kips) 

47,880.26 

pascals 

tons  per  square  foot  (tsf) 

95,760.52 
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1  Introduction 


Purpose 

The  Earthquake  Engineering  Site  Characterization  -  Research  Needs 
Workshop  was  held  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  on  27-28  June  1995.  The  purpose  of  the  workshop  was  to  bring  together 
technical  specialists  from  WES,  Corps  of  Engineers  Districts,  the  US  Bureau  of 
Reclamation  (USER)  and  academia  to  solicit  input  for  initial  planning  of  a  multi¬ 
year  project  intended  to  develop  improved  geophysical  techniques  for  earthquake 
engineering  site  characterization. 

Problem  Definition 

The  ultimate  goal  of  the  sponsoring  Earthquake  Engineering  Research 
Program  (EQEN)  work  unit  is  to  develop  tools  to  assist  Corps  of  Engineers 
agencies  responsible  for  safe,  cost-effective  design  of  reservoir  dams  and  other 
critical  water  resource  and  navigation  facilities  in  their  determination  of  accurate 
engineering  site  condition  parameters.  The  focus  of  the  work  unit  is  on  in  situ 
technologies. 

Approach 


The  subject  workshop  was  convened  at  an  early  stage  in  the  investigation,  to 
identify  candidate  techniques  for  detailed  evaluation  and  to  facilitate  informed 
planning  of  project  direction.  Primary  activities  to  follow  included  selection  of  a 
field  test  site  and  performance  of  field  experiments  in  FY  1996.  Future  project 
activities  were  to  be  planned  based  on  results  from  the  initial  experimentation 
program. 

Results  (Anticipated) 

In  the  near  term,  the  objectives  to  be  realized  from  the  workshop  and 
subsequent  research  in  the  sponsoring  work  unit  included: 
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(1)  consensus  of  attendees  on  the  definition  and  causes  of  liquefaction  during 
earthquakes; 

(2)  selection  of  desirable  characteristics  of  a  test  site  for  evaluation  of 
geophysical  techniques;  and 

(3)  determination  of  criteria  for  test  site  credibility. 

Ultimately,  the  work  unit’s  goals  include,  but  are  not  lumted  to; 

(1)  reliable  determination  of  indicators  that  liquefaction  has  occurred  or  will 
occur  in  the  future; 

(2)  determination  of  any  correlations  between  geophysical  and  geotechmcal 
parameters; 

(3)  definition  of  seismic  methods  for  site  characterization  using  energy  loss 
and  other  parameters,  and,  generally; 

(4)  improvement  of  geophysical  site  characterization  procedures. 

Format  of  Proceedings  Report 

This  report  summarizes  issues  raised  by  each  presenter  during  the  lecture 
sessions  of  the  workshop,  in  order  of  presentation.  The  primary  techmcal  content 
of  each  lecture  is,  of  course,  found  in  papers  submitted  by  the  presenters  at  or 
following  the  workshop;  all  papers  submitted  are  provided  in  appendices  to  this 
proceedings  report;  a  transcript  of  notes  taken  by  the  authors  of  this  report  during 
the  workshop  is  also  provided  in  the  beginning  text.  It  was  also  the  intent  of  the 
proceedings  authors  to  reconstruct  audience  discussion  as  faithfully  as  possible. 
Electronic  media  versions  were  not  available  for  all  papers;  in  these  instances, 
copies  of  documents  provided  by  presenters  were  scanned  and  adapted  into  a 
consistent  format  for  publication.  Photographic  images  and  graphic  illustrations 
were  not  scanned;  these  were  photocopied  and  placed  at  the  end  of  each  paper.  A 
list  of  attendees  and  their  addresses  at  the  time  of  the  workshop  is  provided  in 
Appendix  A.  A  transcript  of  the  meeting  agenda,  developed  from  a  handout  made 
available  to  attendees  at  the  workshop,  is  provided  in  Appendix  B . 
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2  Presentations 


The  workshop  was  opened  with  a  welcoming  sBtement  by  Dr  William  F 
Mm^son  m.  Direaor  of  fte  Geotechmcal  Laborahny,  fol loL  by  a  ZLy 

DmS.  L  ,h  F  Ob  “f Engineering  and  Geosciences 
Division  on  the  Earthquake  Engineenng  Research  Program  (EQEN)  which 

provided  fondmg  support  for  the  workshop.  The  EQEN  is  itself  funded  by  direct 

^0  ment  from  Headquarters,  US  Army  Corps  of  Engmeers,  through  the  Civil 

^  w  Development  Program.  Dr.  Dwain  K.  Butler  and  Mr. 

•  +  ?  J'  ^  “troduced  the  workshop  participants  and  described  the 

diS2d^S)T^  technical  presentations  are 

Topic  I:  Geotechnical  Indicators  of  Liquefaction 
Potential  of  In  Situ  Materials 

Professor  W.  Liam  Finn,  University  of  British  Columbia,  presented  the 
openmg  lecture  of  the  workshop,  describing  the  fundamental  issues  of  liquefaction 
Aat  are  of  most  concern  to  engineers.  In  addition  to  his  presentations  at  fte 
workshop,  Professor  Finn  submitted  two  techmcal  papers,  which  are  mcluded  in 
Appendix  C  to  these  proceedmgs.  The  followmg  several  paragraphs  are  a 
transcnpt  of  notes  taken  durmg  Professor  Finn’s  lecture. 

Professor  Finn  prefaced  his  presentation  by  posing  the  following  three 
qu^tions  md  one  task  of  geotechnical  engineering  concern,  namely: 

•What  IS  the  potential  for  pore  water  pressure  buildup? 

•What  is  the  response  of  the  structure? 

•How  has  and  how  much  has  the  structure  moved,  deformed,  ...etc? 

Look  for  indicators  of  movement  and  overall  response  during  and  after  the 
earthquake. 

The  troublesome  areas,  as  Prof  Finn  saw  them,  were  (a)  measures  to  deal 
economcally  with  the  occurrence  and  consequences  of  liquefaction,  and  (b)  the 
cost  of  remediation.  ^  ^ 
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The  conventional  view  of  liquefaction,  according  to  Prof  Finn,  involves  a 
three  step  concept; 

(1)  Before  the  earthquake;  sand  may  be  loose,  and  sand  grains  are 

intercoimected  with  each  other; 

(2)  During  the  earthquake;  sand  grains  float  in  groundwater,  and  the  water  and 

sand  are  forced  to  the  surface,  where  sand  and  water  are  observed  to  boil  from 

the  ground;  u  a  a 

(3)  After  the  earthquake;  sand  is  more  compacted,  the  ground  has  subsided, 

and  water  has  been  removed  from  within  the  sand  structure. 

Relating  to  the  use  of  geotechnical/geophysical  techniques  to  indicate  potential 
for  liquefaction  occurrence.  Prof.  Finn  observed  that  liquefaction  is  accompanied 
by  notable  modifications  in  frequency  and  signature  in  tlie  ground  motion  record; 
he  suggests  that  p-wave  velocity  may  be  a  stable  indicator.  He  noted  that  blast- 
induced  liquefaction  induces  surface  settlements  on  the  order  of  10  percent, 
whereas  earthquake-induced  liquefaction  generally  causes  settlements  nearer  to 
about  2  percent.  Blasting  appears  to  break  down  contacts  between  particles, 
allowing  the  soil  to  achieve  a  denser  state. 

At  shallow  depths  (less  than  about  30  m),  characterization  of  the  effects  of 
(previous)  earthquakes  on  a  site  may  be  a  good  measure  of  liquefaction  potential. 
These  effects  may  be  characterized  by  visual  observation  of  surface 
manifestations,  or  by  use  of  in  situ  investigation  tools,  such  as  Standard  or  Cone 
Penetration  Tests  (SPT’s  or  CPT’s,  respectively).  These  tests  are  strongly 
affected  by  density,  and  Prof  Finn  maintained  that  liquefaction  does  not  depend 
solely  on  density;  stress-versus-strain  behavior,  as  tracked  by  shear  modulus,  may 
be  a  prime  indicator  of  liquefaction  (at  least  in  theory).  Prof  Finn  emphasized 
that  surface  mamfestations  of  liquefaction,  though  most  obvious,  are  not  always 
of  engineering  sigmficance;  rather,  high  excess  pore  pressures  and  concomitant 
shear  strength  reduction  may  be  more  important.  He  cautioned  attendees  agamst 
restricting  their  concern  only  to  100%  pore  pressure  response,  as  damage  may 
result  from  pore  pressure  development  at  lower  levels.  His  concluding  question  to 
the  audience  was  ‘What  is  the  pore  water  pressure  doing  if  we  have  saturated 
material,  where  boundary  conditions  and  permeability  are  conducive  to  sustaining 
excess  pore  water  pressures  caused  by  earthquake  shaking? 

Prof  Finn  challenged  the  attendees  to  consider  changing  contemporary 
defensive  design  of  dams  to  restrict  deformations  at  key  points,  rather  than  by 
limiting  allowable  global  deformation  involving  the  footprint  of  the  embankment. 
He  contended  that  a  factor  of  safety  (against  liquefaction)  equal  to  umty  does  not 
necessarily  denote  complete  loss.  He  recommends  designing  the  structure  to  limit 
damage  to  reparable  levels.  Prof  Finn  also  stated  that  geologists  and  geologic 
mapping  information  (e.g.,  percent  fines,  how  much  sedimentation,  etc.)  should  be 
brought  into  an  investigation  to  aid  understanding  of  a  site  pnor  to  detailed  field 
visits  and  exploration. 
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Finally,  Prof.  Finn  discussed  the  use  of  blast-induced  liquefaction  to  judge 
whether  a  site  would  likely  liquefy  during  earthquake  shaking.  He  agrees  that 
blasting  may  indeed  cause  liquefaction  of  susceptible  deposits,  depending  on 
geologic  structure  and  layering,  but  he  maintains  that  much  of  the  response  is 
dependent  on  the  blasting  pattern,  which  could  induce  an  artificially  damapiTig  set 
of  conditions.  He  advised  the  group  to  be  careful  in  the  use  of  dgta  on  re¬ 
liquefaction. 

As  mentioned  earlier.  Professor  Firm  submitted  two  companion  technical 
papers  for  publication  in  these  workshop  proceedings,  titled,  respectively:  “A 
Basis  for  Selecting  Indicators  of  Liquefaction  Potential,”  and  “Evaluation  of 
Liquefaction  Potential  by  In  Situ  Methods.”  Both  papers  are  provided  in 
Appendix  C,  adapted  to  the  format  of  choice  for  the  proceedings. 

Topic  II:  Geophysical  Indicators  of  Liquefaction 
Potential 

Investigation  of  Liquefiable  Soils  Using  Shear  Wave  Velocity  -  Ken 
Stokoe 

Professor  Ken  Stokoe  II,  University  of  Texas,  and  Dr.  Ronald  Andrus, 
National  Institute  of  Standards  and  Technology,  co-presented  a  lecture  on 
“Investigations  of  Liquefiable  Soils  Using  Shear  Wave  Velocity  and  Associated 
Correlations,”  for  which  lecture  graphics  were  reproduced  for  inclusion  in  these 
proceedings.  The  lecture  developed  the  topic  of  seismic  velocity  surveying,  by 
both  borehole  and  surface  means,  and  its  specific  application  to  evaluate 
liquefaction  potential  of  various  soils.  Prof  Stokoe  observed  during  his 
presentation  that  currently  applied  procedures  cannot  determine  whether  a  soil  is 
as  much  as  90%  saturated,  and  that  perhaps  the  ratio  of  compression  wave 
velocity,  Vp,  to  shear  wave  velocity,  V^,  offers  a  reliable  correlation.  He  pointed 
out  that  Vj  is  directly  relatable  to  shear  modulus  and  has  been  correlated  to  SPT 
blowcounts,  thus  it  follows  that  liquefaction  potential  may  be  inferred  from  the 
parameter,  as  well. 

Prof  Stokoe  listed  a  series  of  engineering  parameters  and  behavior  indices  that 
may  be  directly  or  indirectly  evaluated  in  terms  of  seismic  velocity,  and  presented 
field  occurrence  data  on  liquefaction  to  prove  his  point.  Shear  wave  velocity  is 
gaining  acceptance  rapidly  in  the  geotechnical  earthquake  engineering  field,  based 
on  successful  application  to  observed  soil  response  at  well-studied  sites. 

Professor  Stokoe  and  Dr.  Andrus  presented  the  theory  and  application  of  the 
Spectral  Analysis  of  Surface  Waves  (SASW)  technique  (developed  by  Prof 
Stokoe  and  his  colleagues  at  the  University  of  Texas),  and  described  recent 
adaptation  of  the  technique  to  liquefaction  potential  evaluation.  Copies  of  lecture 
graphics  are  provided  in  Appendix  D  for  detailed  examination  of  these  data  by  the 
reader. 
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Seismic  Velocity,  Attenuation  and  Liquefaction  -  Philip  Sirles 

Mr.  Philip  Sirles,  US  Bureau  of  Reclamation,  gave  a  background  discussion 
of  the  purposes  and  techniques  for  geophysical  investigations  applied  to  evaluate 
engineering  sites  for  liquefaction  potential.  He  did  not  submit  a  formal  techmcal 
paper  to  the  workshop;  hard  copies  of  his  lecture  graphics  and  a  copy  of  an  article 
Mr.  Sirles  published  in  a  newsletter  for  the  US  Committee  on  Large  Dams  titled 
“Geophysics  for  Soil  Characterization  at  New  and  Existing  Dams”  are  included  in 
the  text  of  the  proceedings,  immediately  following  this  synopsis. 

Mr.  Sirles’  presentation  emphasized  safety  evaluation  of  reservoir  dams,  in 
keeping  with  the  focus  of  the  sponsoring  EQEN  research  program.  He  discussed 
the  entire  range  of  field  investigation  techniques  available  to  engineers,  including 
those  employed  at  the  ground  surface,  such  as  seismic  refraction  and  reflection 
surveys,  SASW  methods,  electrical  property  measurement  and  electromagnetic 
survey  approaches.  Mr.  Sirles  also  presented  the  principles  and  practices 
associated  with  borehole-based  techniques.  His  presentation  was  organized  to 
then  integrate  the  suite  of  available  tools  into  a  method  for  liquefaction  potential 
investigation,  culminating  in  case  studies  of  two  existing  dams  where  the  process 
has  been  employed.  His  presentation  graphics  are  provided  in  Appendix  E  and 
contain  substantial  detail  to  convey  the  salient  features  of  his  lecture. 

A  Non-Destructive  Method  of  Estimating  Properties  of  Soil 
Subjected  to  Earthquake  Strong  Motion  -  Steven  Glaser 

Prof  Steven  Glaser,  Colorado  School  of  Mines,  Golden,  Colorado,  presented 
the  subject  lecture  and  submitted  a  companion  technical  paper  with  the  same  title. 
His  topic  addressed  the  concept  and  application  of  system  identification,  based  on 
recorded  site  response  during  earthquakes,  to  evaluate  dynamic  properties.  Prof 
Glaser  presented  the  scheme  used  to  acquire  and  analyze  data  at  the  instrumented 
Wildlife  Liquefaction  Site  in  southern  California  during  two  recent  earthquake 
events.  Prof  Glaser  discussed  the  interpretation  of  motion  records  between  in  situ 
accelerometers,  leading  to  parametric  modeling  of  time-varying  soil  properties. 
Complete  pore  pressure  histories  were  also  obtained  by  transducers  implanted  at 
the  site.  The  system  identification  process  allows  the  use  of  a  reed  geologic  site 
for  experimentation,  whereas  boundary  conditions  and  loads  in  existence  are  not 
well  simulated  in  an  elemental  representation  of  the  actual  site  by  means  of 
laboratory  tests.  Prof  Glaser’s  presentation  graphics  and  technical  paper  are 
provided  in  Appendix  F  and  provide  additional  detail  on  the  concept  and  the  case 
studies  just  described. 

Surface  Wave  Method  -  A  Tool  for  Lifeline  Earthquake  Engineering  - 
Ronald  Andrus 

Dr.  Ronald  Andrus,  National  Institute  of  Standards  and  Technology,  presented 
and  submitted  a  technical  paper  titled  as  above,  to  describe  application  of  the 
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SASW  technique  to  liquefaction  potential  evaluation  in  general,  and  assessment  of 
the  effects  of  liquefaction  on  lifelines  specifically.  The  SASW  technique  is 
attractive  in  this  application  due  to  its  ability  to  assess,  non-intrusively,  large 
areas  quickly  and  economically  without  requiring  boreholes.  Dr.  Andrus 
described  the  use  of  the  SASW  procedure  to  discriminate  liquefied  and  non- 
liquefied  soils  in  the  areas  affected  by  the  1983  Borah  Peak,  Idaho  earthquake. 

Few  discussion  comments  were  made  by  attendees  during  the  workshop;  the  brief 
techmcal  paper  is  provided  in  Appendix  G  and  provides  some  details  for  the 
reader  to  understand  the  concepts  and  case  studies  presented. 

The  Biot  Model  and  Waveform  Tomography  -  Richard  Rechtien 

Prof  Rechtien,  University  of  Missouri  -  Rolla,  presented  a  background  lecture 
on  the  Biot  Model  applied  to  soil-water  continua  and  the  hydrodynamic 
interactions  between  the  soil  solid  and  liquid  phases,  for  the  interpretation  of 
seismic  waveforms  transmitted  through  geologic  media  between  source  and 
receiver  boreholes.  Waveform  tomography  provides  multidimensional 
interrogation  and  modeling  of  the  ground,  and  involves  collection  and  processing 
of  a  very  large  quantity  of  high-resolution  motion  data.  The  result  of  seismic 
surveys  conducted  in  this  manner  may  be  thought  of  as  a  computer-aided 
tomographic,  or  CAT  scan,  of  the  subsurface.  The  parameters  of  the  particular 
scan  may  vary,  to  the  point  that  dynamic  soil  properties  may  be  mapped  and 
perhaps  liquefaction  potential  may  be  inferred  in  a  two-dimensional  fashion.  It 
may  also  be  possible  to  extract  permeability  data,  through  adaptation  of  the  Biot 
Model.  In  any  case,  two-  or  three-dimensional  measurement  of  shear  wave 
velocities  is  possible  applying  the  tomography  procedure,  and  shear  wave 
velocity,  as  has  already  been  discussed,  is  well  correlated  to  liquefaction  potential. 

Prof  Rechtien  submitted  no  presentation  visuals  or  technical  papers  to  the 
workshop  for  inclusion  in  these  proceedings.  Prof  Rechtien  published  a  WES 
techmcal  report.  Contract  Report  No.  GL-96-1,  “In  Situ  Seismic  Investigation  of 
Liquefaction  Potential  of  Soils,”  after  the  workshop,  which  provides  excellent 
descriptions  of  theory,  field  data  acquisition,  and  the  state  of  the  art  for 
application  of  the  waveform  tomography  method  to  the  interpretation  of 
liquefaction  potential. 

Summary  of  Topics  I  and  II  -  Dwain  Butler 


Dr.  Dwain  Butler,  Geotechnical  Laboratory  presented  a  summary  of  the  key 
issues  discussed  in  the  previous  portions  of  the  workshop,  begmning  by  echoing 
back  to  attendees  what  he  understood  from  earlier  presentations.  Dr.  Butler 
outlined  the  liquefaction  concept  for  a  simple,  level-ground,  three-layer  deposit 
using  the  following  illustration; 
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In  this  diagram  the  soils  are  assumed  to  be  sandy  or  silty  deposits  of  the  types 
conventionally  held  to  be  susceptible  to  liquefaction.  It  is  generally  recognized 
that  only  saturated  soils  can  liquefy,  or,  in  an  extreme  sense,  soils  that  could 
become  saturated  by  migration  of  excess  pore  water  pressures  from  adjacent  soils. 
If  the  soils  in  layer  three  have,  either  in  the  laboratory  or  in  the  field,  ever 
experienced  liquefaction  when  subjected  to  a  given  earthquake  loading 
(characterized  by  a  peak  acceleration,  a„,ax,  a  fundamental  period,  T,  and  a  number 
of  cycles,  N),  then  the  soil  is  considered  potentially  liquefiable  in  future  events. 

If,  on  the  other  hand,  the  same  soil  has  not  liquefied  during  real  events  or  in  the 
laboratory  when  subjected  to  any  set  of  earthquake  motions,  the  soil  is  considered 
non-liquefiable. 

Engineering  properties  to  be  considered  in  characterizing  the  liquefaction 
potential  of  a  soil  include,  but  may  not  be  limited  to:  density  (bulk  or  relative), 
shear  or  compression  wave  velocity,  attenuation  quality  factors  for  shear  or 
compression,  and  the  in  situ  stress  conditions  in  effect  prior  to  the  earthquake. 

The  earthquake  ground  motions  may  be  characterized  by  the  parameters 
mentioned  earlier,  which  are  themselves  functions  of  such  parameters  as 
magnitude,  focal  depth,  epicentral  distance,  and  geometric  constraints  such  as  the 
type  of  fault  and  the  energy  radiation  pattern  from  the  source.  The  earthquake 
characteristics  are  generally  independent  of  the  properties  of  the  soil  in  layer  three 
as  idealized,  and  the  soil  properties  in  layer  three  are,  in  general,  independent  of 
the  earthquake  shaking  at  its  onset.  Of  course,  the  effects  of  ground  shaking  on 
the  response  of  the  soil  are  strain  dependent  during  shaking. 

Dr.  Butler  reiterated  the  task  for  earthquake  scientists  concerned  with 
liquefaction  potential  evaluation,  namely;  to  characterize  the  liquefaction  potential 
of  soil  layer  three  using  techniques  that  determine  properties  of  the  layer, 
independent  of  the  earthquake,  then  relate  them  to  liquefaction  potential  by  means 
of  a  plot  or  chart  (based  on  substantial  accumulation  of  data).  It  may  be  possible, 
for  example,  to  develop  a  relationship  between  geophysical  and  geotechnical 
measurements  and  liquefaction  potential.  Shear  wave  velocity  or  peak 
acceleration  may,  for  instance,  be  relatable  to  liquefaction  potential  using  strength 
parameters  such  as  may  be  interpreted  from  standard  penetration  test  blowcounts 
(for  which  a  large  body  of  data  exists).  The  evolution  of  in  situ  test  procedures 
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may  allow  for  correlation  of  additional  types  of  data  with  liquefaction  potential. 
Dr.  Butler  listed  a  variety  of  seismic  survey  techniques  for  obtaining  velocity  data, 
including,  crosshole  or  surface  refraction  surveys  and  their  adaptation  for 
tomographic  interpretation,  uphole  or  downhole  seismic  surveys,  the  SASW 
approach,  and  reflection  seismic  surveys.  Electrical  surveying  techniques  are  also 
available  to  interrogate  the  subsurface,  including  electromagnetic  measurements, 
resistivity,  and  ground  penetrating  radar,  where  applicable. 

Topic  ill:  Field  Test  Sites 

Mr.  Ballard  led  the  concluding  workshop  session  on  selection  of  documented 
test  sites  for  research  on  site  characterization  techniques,  by  stating  the  primary 
objectives  of  the  session,  namely:  (1)  to  examine  potential  test  sites  for  relevancy 
to  development  of  geophysical  techniques;  (2)  to  estimate  test  site  criteria;  and  (3) 
to  select  and  prioritize  sites  in  accordance  with  accepted  criteria.  He  followed  his 
summary  of  the  objectives  by  introducing  the  session  speakers. 

Organization  and  Operation  of  National  Geotechnical  Test  Sites  - 
Richard  Woods 

Prof  Richard  D.  Woods,  University  of  Michigan,  presented  the  National 
Geotechnical  Experimentation  Sites  (NGES),  which  he  manages  and  which  are 
intended  to  accommodate  researchers  conducting  experiments  on  in  situ 
investigation  tools  and  procedures.  The  NGES  consists  of  five  primary  sites  that 
vary  in  their  type  and  amount  of  instrumentation  and  background  data  and  which 
are  located  on  property  belonging  to:  the  University  of  Houston,  the  University  of 
Massachusetts,  Northwestern  University,  Texas  A&M  University,  and  the 
Treasure  Island  Naval  Station  offshore  from  San  Francisco,  California.  The 
National  Science  Foundation  (NSF)  and  the  Federal  Highway  Administration 
(FHWA)  provide  research  funding  in  support  of  the  five  primary  sites.  The  Texas 
A&M  and  Treasure  Island  sites  are  categorized  as  Level  I  sites,  and  receive  the 
greatest  amount  of  supportive  funds  for  site  improvement  and  user  access  and 
activities,  while  the  three  remaining  primary  sites  are  funded  at  lesser  amounts 
(designated  Level  II).  Thirty-five  additional  NGES  sites  are  included  in  the 
overall  system  database,  but  are  not  currently  funded  for  instrumentation  or 
maintenance.  Summary  descriptions  of  each  of  the  five  primary  sites  are  given  in 
the  technical  paper  submitted  by  Prof  Woods  to  the  workshop,  which  is  provided 
in  Appendix  H. 


Prof  Woods  also  listed  several  international  sites  for  consideration  as 
experimentation  facilities,  in  the  following  countries:  Armenia,  Brazil,  Canada, 
France,  Italy,  Japan,  the  Netherlands,  and  the  United  Kingdom.  He  mentioned' 
that  a  catalog  has  been  published  on  diskettes,  with  a  user’s  manual.  Interested 
persons  should  contact  Prof.  W^oods  for  copies.  An  online  access  feature  was 
under  development,  but  was  not  ready  at  the  time  of  the  workshop. 
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Effect  of  Stress  State  and  Soil  Fabric  on  Seismic  Wave  Velocities 
Measured  by  Borehole  Methods  -  Ken  Stokoe 

Prof.  Stokoe  returned  to  present  the  subject  topic  in  the  session  relating  to  test 
site  selection,  and  presented  research  results  from  large-scale  laboratory  triaxial 
chamber  tests  on  fabric  and  soil  structure  effects.  His  presentation  concluded 
with  crosshole  test  program  that  included  a  case  study  on  comparisons  of 
horizontally  and  vertically  polarized  shear  waves  from  crosshole  tests  at  the 
Treasure  Island  NGES  facility.  No  technical  paper  was  submitted  to  accompany 
this  presentation;  Prof  Stokoe’s  lecture  graphics  for  this  topic  are  provided  in 
Appendix  I. 

Success  Dam  -  John  Nickell 

Mr.  John  Nickell,  US  Army  Engineer  District,  Sacramento  (CESPK-ED-D) 
made  an  impromptu  presentation  on  Success  Dam,  California,  for  consideration 
by  attendees  as  an  additional  field  test  site.  WES  has  assisted  CESPK  in  an 
evaluation  of  seismic  stability  of  Success  Dam,  where  a  substantial  in  situ 
investigation  program  has  been  conducted.  Mr.  Nickell  described  conditions  and 
geology  of  the  Success  Dam  site,  as  well  as  preliminary  stability  analysis  results. 
The  length  of  study  areas  involving  the  embankment  of  this  dam  ranges  from  3000 
ft  to  4000  ft,  depending  on  inclusion  of  auxiliary  dikes  around  the  reservoir.  The 
main  embankment  is  140  ft  tall,  and  is  constructed  of  compacted  earth  and  rock 
fill.  Recent  alluvium  deposits  left  in  place  beneath  the  shells  of  the  dam  range  in 
thickness  from  about  1 1  ft  to  25  ft.  Mr.  Nickell  described  the  suite  of  field  and 
laboratory  tests  conducted  to  assess  liquefaction  potential  of  the  dam  and  its 
foundation,  and  related  plans  for  additional  experiments  to  further  refine  site 
characterization.  In  presenting  the  features  and  test  programs  at  Success  Dam, 
Mr.  Nickell  offered  access  and  support  to  research  personnel  for  an 
experimentation  site  that  involved  an  actual  reservoir  dam.  Persons  interested  in 
additional  details  on  the  investigations  sponsored  by  CESPK  at  Success  Dam 
should  contact  Mr.  Nickell  at  the  address  given  in  the  list  of  attendees  provided 
earlier.  (Note:  Success  Dam  was  chosen  as  the  first  test  site.  Field  investigations 
were  conducted  during  the  summer  of  1996.) 

Test  Sites  Discussion  -  Robert  Ballard 

Mr.  Ballard  presided  over  an  open  discussion  of  the  relative  merits  of  various 
test  sites  that  may  be  available  for  site  characterization  field  research,  and  invited 
the  inclusion  of  additional  sites  that  had  not  yet  been  discussed  in  the  workshop. 
This  discussion  was  directed  to  aid  the  principal  investigators  in  selection  of 
candidate  test  sites  for  experiments  to  commence  in  FY  96.  Attendees  were  polled 
to  secure  their  informed  opinions  on  potential  test  sites  discussed  to  this  point,  and 
to  add  any  new  candidate  sites  that  might  prove  useful  beyond  those  originally 
considered.  Mr.  Ballard  circulated  a  blank  table  to  attendees  early  m  the 
workshop  to  solicit  opinions  on  the  pros,  cons,  and  priority  of  each  test  site  name 
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submitted.  The  consensus  results  obtained  during  the  closing  session  of  the 
workshop  are  summarized  in  Table  1 . 


3  Epilogue 


In  the  course  of  time  since  these  workshop  proceedings  were  compiled,  certain 
events  worthy  of  note  have  occurred  durmg  the  progression  of  this  Work  Unit.  A 
chronology  of  those  events  follows; 

1 .  A  lecture  /  workshop  was  conducted  at  WES  on  the  subject  of  seismic 
tomography  and  Biot  theory.  The  featured  lecturer  was  Professor  Tokuo 
Yamamoto,  University  of  Miami. 

2.  Dr.  Richard  Rechtien’s  report  “In  Situ  Seismic  Investigation  of 
Liquefaction  Potential  of  Soils”  was  published. 

3 .  Applications  of  the  Biot  theory  as  it  pertains  to  the  objectives  of  this  Work 
Unit  were  documented. 

4.  Field  tests  were  planned  and  conducted  at  Success  Dam,  California.  The 
borehole  pattern  was  an  equilateral  triangle  twenty-five  feet  long  on  each  side.  A 
single  boring  ninety  feet  deep  was  placed  at  each  apex  and  ten  foot  deep  borings 
were  placed  along  each  leg  at  2.5  foot  spacing.  Thus,  a  total  of  thirty  borings 
were  used  durmg  testing.  High  quality  data  were  collected  in  a  fashion  suitable 
for  tomographic  analysis.  These  data  will  be  analyzed  during  FY  97. 

5.  Dr.  Barbara  Luke,  University  of  Nevada,  Las  Vegas,  was  hired  under  the 
ARO  summer  faculty  program  to  plan  and  assist  with  potential  Spectral  Analysis 
of  Surface  Waves  (SASW)  applications  at  future  test  sites.  She  worked  at  WES 
during  the  time  June-August  1996. 

6.  The  Work  Unit  was  successfully  defended  during  field  review  and  the 
concept  of ‘Tlon-Invasive  Characterization  of  Upstream,  Water-Covered  Dam 
Embankments  and  Foundations”  was  proposed  as  either  a  new  start  or  inclusion  in 
the  existing  Work  Unit.  (The  proposal  was  later  included  as  a  part  of  the  existine 
Work  Unit). 

7.  Lastly,  the  U.S.-Japan  Workshop  on  Advanced  Research  on  Earthquake 
Engineering  for  Dams  was  held  at  WES  on  12-14  November  1996.  A  paper 
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“Geophysical  Methods  for  Site  Characterization  and  Measurement  of  In  Situ 
Properties”  by  Robert  F.  Ballard,  Jr.  was  presented  which  is  available  separately 
and  summarizes  the  features  and  progress  of  the  Work  Unit. 
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Table  1 
Potential  Field  Test  Sites 


T reasure  Island 
(NGES) 


University  of 
Houston 
(NGES) 


Texas  A&M 
(NGES) 


Northwestern 

University 

(NGES) 


Unity  Dam 
(USSR) 


Conconullay  Dam' 
(USSR) 


Bradbury  Dam 
(USBR) 


Jackson  Lake  Dam 
(USBR) 


Success  Dam'* 
(CESPK) 


Well  documented 

Earthquake 

emphasis 


Man-made  loose 
sand  deposit 
Bay  mud 

Small,  salt  water  site 


Doesn’t  meet 
requirements 


Not  good  for  long-  Low 
term  study 
Denser  than 
required 


Limited  access 


Time  factor 
Complex  geology 


Partially  de-watered 


Remediated 
Limited  access 


Accessible  Borehole  placement  High/Med 

Sands  and  gravels  Affected  by  four 

Documentation  major  faults 

Ongoing  WES  work  Cemented  lenses 

(possible  additional 

funding) 


Loose  sands 


Clean  sand 

Documentation 

Liquefiable 


Documentation 
Few  fines,  clean 
sand, 

some  gravels 


Gravelly  site 
Instrumented  for 
strong  motions 
(SMIP)^ 


High 

POC:  Mike  Stevens^ 


Liquefiable  Boreholes  not 

Earthquake  likely  to  accessible 
occur  Difficult  site  access 

Other  test  run  there 


High/Med 


Second  choice  of  group 

2 

Mr .  Stevens  was  POC  at  the  time  of  the  workshop  -  he  has  since  left  the  USBR 

3 

SMIP  -  Corps  of  Engineers'  Strong  Motion  instrumentation  Program,  managed  by  Robert  F  Baliard  Jr 
CEWES-GG  ’  ’ 


First  choice  of  group 


Appendix  A  -  List  of  Attendees 


The  following  participants  attended  from  outside  WES  and  either  presented 
technical  lectures  on  workshop  topics,  or  prepared  and  submitted  technical  papers 
for  the  proceedings  report,  or  both  (addresses  given  below  were  correct  at  the  time 
of  the  workshop): 


Name 

Organization 

Address/Phone  (fax) 

Dr.  Ronald  D.  Andrus 

National  Inst,  for  Standards  and 
Technology 

NIST 

Bldg.  226,RmB158 
Gaithersburg,  MD 

20899 

301-975-6051  (fax  301- 
869-6275) 

Prof.  W.  D.  Liam  Finn 

Univ.  of  British  Columbia 

Dept,  of  Civil 

Engineering 

2324  Main  Mall 
Vancouver,  BC 

CANADA  V6T  1W5 
604-822-4938  (fax  604- 
822-6901) 

Prof.  Steven  Glaser 

Colorado  School  of  Mines 

Division  of  Engineering 
Golden,  CO  80401 
303-236-4196  (fax  303- 
273-3602) 

Prof  Richard  Rechtien 

Consultant,  Univ.  of  Missouri- 
Rolla 

Dept,  of  Geology  and 
Geophysics 

Rolla,  MO  65401 
573-364-3886 
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Mr.  Phil  Sirles 

U.S.  Bureau  of  Reclamation 

D-2230 

PO  Box  25007 

Bldg.  67,  Denver  Federal 
Center 

Denver,  CO  80225 

Prof  Keimeth  H.  Stokoe  II 

University  of  Texas  at  Austin 

ECJ  9.227 

Austin,  TX  78712 
512-471-4929  (fax  512- 
471-6548) 

Prof  Richard  Woods 

University  of  Michigan 

Dept,  of  Civil 

Engineering 

2360  G.  G.  Brown  Lab 
Ann  Arbor,  MI  48109- 
2125 

313-764-4303  (fax  313- 
764-4292) 

The  following  participants  from  outside  WES  attended  lectures  and  discussions, 
but  did  not  make  formal  presentations  or  prepare  technical  papers  for  the 
proceedings; 


Name 

Organization 

Address/Phone  (fax) 

Mr.  William  E.  Hancock 

US  Army  Engineer  District, 

Seattle 

CENPS-EN-GT 

PO  Box  3755 

Seattle,  WA  98124-2255 
206-764-3711  (fax  206- 
764-6795) 

Dr.  Gregory  L.  Hempen 

US  Army  Engineer  District,  St. 
Louis 

CELMS-ED-GG 

1222  Spruce  Street 

St.  Louis,  MO  63103- 
2833 

314-331-8441  (fax  314- 
331-8244) 

Mr.  William  Joyner 

US  Geological  Sxirvey 

Mail  Stop  977 

345  Middlefield  Road 
Menlo  Park,  CA  94205 
415-329-5640  (fax  415- 
329-5163) 
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Mr.  Kenneth  Klaus 

US  Army  Engineer  District, 
Vicksburg 

CELMK-ED-GT 

4155  Clay  Street 
Vicksbm-g,  MS  39180- 
3435 

601-631-7122  (fax  601- 
631-5627) 

Mr.  John  Nickell 

US  Army  Engineer  District, 
Sacramento 

CESPK-ED-D 

1325  J.  Street 

Sacramento,  CA  95814- 
2922 

916-557-6614  (fax  (916- 
557-7846) 

The  following  WES  participants  either  presented  lectures  or  attended 
discussions  during  the  workshop; 


Name 

Organization 

Address/Phone  (fax) 

Mr.  Robert  F.  Ballard,  Jr. 
(Principal  Investigator) 

U.S.  Army  Engineer  Waterways 
Experiment  Station 

CEWES-GG 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-2201  (fax - 
3153) 

Dr.  Dwain  K.  Butler 
(Co-Principal  Investigator) 

WES 

CEWES-GG 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-2127  (fax - 
3453) 

Dr.  Arley  G.  Franklin 

WES 

Chief,  CEWES-GG 

3909  Halls  Ferry  Road 
Vicksbm-g,  MS  39180- 
6199 

601-634-2658  (fax - 
3453) 

Dr.  Joseph  P.  Koester 

WES 

CEWES-GG-H 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-2202  (fax  - 
3453) 
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Mr.  Richard  H.  Ledbetter 

WES 

CEWES-GG-H 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-3380  (fax - 
4894) 

Dr.  John  Peters 

WES 

CEWES-GS-GC 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-2590  (fax - 
4656) 

Mr.  Ronald  Wahl 

WES 

CEWES-GS-GC 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-3632  (fax- 
4656) 

Mr.  Donald  E.  Yule 

WES 

CEWES-GG-H 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180- 
6199 

601-634-2964  (fax- 
3453) 
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EARTHQUAKE  ENGINEERING 
SITE  CHARACTERIZATION  WORKSHOP 
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Introduction 


One  of  the  more  significant  factors  leading  to  ground  failure  during  earth¬ 
quakes  is  the  liquefaction  of  loose-  to  medium-dense  sands  below  the  water 
table.  Attention  was  focused  on  this  problem  for  the  first  time  as  a  result  of  the 
widespread  ground  failures  during  the  1964  earthquake  in  Niigata,  Japan  (Kawa- 
sumi  1968).  Most  of  the  damage  in  Niigata  attributable  directly  to  the 
earthquake  was  associated  with  liquefaction  and  such  damage  has  been  a  signif¬ 
icant  factor  in  most  major  earthquakes  since  then.  Some  recent  examples  are  the 
extensive  damage  caused  by  liquefaction  in  the  Marina  District  in  San  Francisco 
during  the  Loma  Prieta  earthquake  of  1989  (Seed  et  al.  1990)  and  the  destruction 
of  port  facilities  on  Port  and  Rokko  Islands  in  Osaka  Bay  off  Kobe  City,  during 
the  1995  Hyogo-ken-Nanbu  earthquake  (CAEE  1995). 

Figure  1  illustrates  the  response  of  granular  soils  to  cyclic  loading  in  triaxial 
tests.  The  test  specimens  were  prepared  to  relative  densities  that  covered  the 
range  of  sand  states  from  highly  contractive  to  highly  dilative  under  shearing 
deformation  (Vaid  and  Chem  1985).  All  test  specimens  carry  initial  static  shear 
stresses,  This  simulates  the  most  common  situation  in  earth  structures  and 
under  foundations.  Contractive  stress-strain  response  is  shown  in  Figure  1(a). 

At  some  stage  in  the  cyclic  loading,  strain  softening  occurs  leading  to  a  large 
drop  in  undrained  strength  which  then  remains  constant  over  a  large  range  in 
strain.  Large  deformations  occur  when  is  greater  than  the  steady  state 
strength  (or  residual  strength,  S„r).  Deformations  may  also  be  significant 
when  Tj,  <  Tjs,  if  Tj,  +  x^y  >  x^^,  where  is  the  cyclic  shear  stress.  The 

phenomenon  of  large  flow  deformations  occurring  at  constant  undrained  strength 
is  termed  liquefaction. 
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At  the  other  extreme,  the  sand  is  strongly  dilative  and  deformation  is  the 
result  of  progressive  stiffness  degradation  in  the  sand  due  to  the  buildup  in 
porewater  pressure  with  cycles  of  loading.  This  phenomenon  is  called  cyclic 
mobility  and  is  illustrated  in  Figure  1(b). 

Between  these  two  clearly  defined  states  is  a  transitional  region  combining 
features  of  both  contractive  and  dilative  response.  This  type  of  behavior  is 
characterized  by  a  loss  in  undrained  strength  due  to  strain  softening  over  a 
limited  range  followed  by  a  regain  of  undrained  strength.  Further  deformations 
result  from  cyclic  mobility.  This  behavior  is  illustrated  in  Figure  1(c).  The 
effective  stress  paths  for  these  tests  are  shown  in  Figures  1(d)  through  1(f).  Note 
the  pick  up  in  shearing  resistance  (o'l  -  o' ^12  with  increasing  strain  in 
Figures  1(e)  and  1(f).  The  strain  histories  of  sands  in  these  three  different  states 
are  shown  together  in  Figures  1(g)  through  1(1). 

The  instability  of  loose  saturated  sands  and  the  residual  strength  after  struc¬ 
tural  collapse  may  also  be  determined  by  static  triaxial  or  simple  shear  tests  in 
the  laboratory.  The  stress-strain  curve  derived  from  a  static  triaxial  test  is  shown 
in  Figure  2.  It  has  the  same  basic  characteristics  as  the  stress-strain  curves  in 
Figure  1. 

In  engineering  practice  the  word  liquefaction  is  usually  associated  with  the 
development  of  high  porewater  pressures  due  to  seismic  loading  irrespective  of 
whether  large  deformations  occur  or  not.  Indeed  the  empirical  procedures  for 
assessing  liquefaction  potential  are  all  based  primarily  on  surface  manifestations 
of  high  porewater  pressure  such  as  sand  boils  without  consideration  of  whether 
the  soil  would  behave  in  a  contractive  or  dilative  mode  during  shearing. 
However,  the  consequences  of  these  high  porewater  pressures  clearly  depend  on 
whether  the  soil  is  contractive  or  not,  and  there  is  great  need  for  a  reliable  in-situ 
indicator  of  these  conditions. 
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Objectives 


In  this  paper  we  will  first  explain  the  basic  phenomenon  of  the  development 
of  porewater  pressure  under  cyclic  loading  which  is  the  fundamental  cause  of 
liquefaction  and  cyclic  mobility.  We  will  examine  what  variables  control  the 
potential  for  high  porewater  pressures  and  consider  what  measures  of  in-situ 
properties  might  identify  this  potential.  Our  primary  interest  here  will  be  in  the 
use  of  in-situ  methods.  In  a  companion  paper  entitled  "Evaluation  of  Liquefac¬ 
tion  Potential  by  In-Situ  Methods  ”  (Finn  1995),  the  state-of-the-art  used  in  cur¬ 
rent  engineering  practice  is  described  in  considerable  technical  detail.  In  this 
paper,  we  will  discuss  the  basic  links  between  the  various  indicators  of  lique¬ 
faction  potential  and  the  soil  parameters  that  control  the  development  of  seismic 
porewater  pressures.  This  framework  gives  a  better  understanding  of  the  indi¬ 
cators  in  current  use  and  provides  the  basis  for  developing  alternative  concepts 
of  howto  evaluate  liquefaction  potential. 

There  are  three  basic  requirements  for  evaluating  liquefaction  potential  of  a 
soil  deposit  during  an  earthquake: 

a.  Characterizing  the  intensity  of  the  earthquake  shaking  of  the  soil. 

b.  Characterizing  the  in-situ  state-of-the-soil. 

c.  Developing  a  criterion  for  liquefaction  occurrence  based  on  indicators  of 
shaking  intensity  and  the  in-situ  condition  of  the  soil. 

We  will,  first  of  all,  characterize  the  intensity  of  shaking. 


Objectives 


3 


Characterizing  Earthquake 
Shaking 


Seismic  shear  stresses  play  a  major  role  in  the  development  of  liquefaction. 
Time  histories  of  shear  stresses  generated  by  earthquakes  are  usually  very  non- 
uniform  and  cannot  be  used  directly  in  empirical  assessments  of  liquefaction. 
Seed  (1983)  suggested  replacing  the  irregular  time  history  by  an  equivalent  num¬ 
ber  of  uniform  cycles  and  normalizing  the  uniform  shear  stress  by  dividing  by 
the  effective  overburden  stress.  Seed  proposed  that  this  uniform  cyclic  shear 
stress  ratio  (CSR)  be  determined  by: 


t  o 

CSR  =  — f-  =  0.65  -7^ 
o  o 

VO  VO 


a 


max 


(1) 


where 


T(,  =  equivalent  average  shear  stress 
=  total  overburden  stress 
o\g  =  effective  overburden  stress 
a„^  =  maximum  surface  acceleration  in  units  of  g 
g  =  acceleration  due  to  gravity 

r^  =  a  reduction  factor  to  account  for  soil  flexibility  and  depth 

In  Japanese  practice,  the  reduction  factor  rd  is  often  approximated  using  the 
following  expression, 
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rj  =  {\  -  0.015Z) 


(2) 


where  Z  =  depth  in  meters. 

Equation  (1)  defines  an  equivalent  level  of  shaking  only.  However,  the 
development  of  porewater  pressure  also  depends  on  the  number  of  cycles  of 
CSR  which  are  equivalent  in  their  effects  to  the  complex  time  history  of  shear 
stresses  developed  during  the  earthquake.  The  equivalent  number  of  cycles 
depends  on  the  duration  of  shaking,  which  in  turn  depends  on  the  magnitude  of 
the  earthquake.  Therefore,  the  intensity  of  shaking  will  be  characterized  by  both 
the  critical  stress  ratio,  CSR,  and  by  the  magnitude,  M,  of  the  earthquake. 
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Mechanism  of  Porewater 
Pressure  Development 


Consider  a  small  cube  of  saturated  sand  as  shown  in  Figure  3(a).  Assume 
that  all  drainage  is  prevented  when  the  sample  is  sheared  in  direct  simple  shear. 
In  one  cycle  of  small  shear  strain,  y,  an  increment  in  porewater  pressure,  Au, 
develops.  Now  consider  a  sample  with  exactly  the  same  granular  structure  and 
density  but  of  completely  free-draining  (Figure  3(b)).  When  subjected  to  one 
cycle  of  the  same  shear  strain,  it  undergoes  a  vertical  strain  Ae^^j.  Since  at  the 
beginning  of  the  cycle,  both  samples  had  the  same  structure,  they  both  have  the 
potential  to  develop  the  same  vertical  strain  increment  Ae^^. 

Since  the  sample  in  Figure  3(a)  did  not  change  in  volume,  the  potential  com- 
pactive  strain  Ae,,^  must  have  been  neutralized  in  some  way.  This  can  only  be 
accomplished,  if  the  sample  rebounds  the  same  amount,  Ae^j.  For  this  to  occur, 
the  effective  stresses  between  the  grains  must  be  reduced  by  the  appropriate 
amount  Ap,  such  that  Ap/E,  =  -Ae  where  is  the  effective  constrained  rebound 
modulus  of  the  soil  for  simple  shear  conditions.  The  decrease  in  Ap  is 
accomplished  by  an  increase  in  the  porewater  pressure  by  an  increment  Au  = 
-Ap.  Therefore,  the  increment  in  porewater  pressure  is  given  by 

A«  =  £  Ae  , 

This  is  the  fundamental  mechanism  of  porewater  pressure  development  during 
earthquakes.  There  is  very  strong  field  evidence  supporting  this  concept  of  how 
seismic  porewater  pressures  that  may  lead  to  liquefaction  are  generated. 

Examples  from  Post-Earthquake  Studies 


Mechanism  of  Porewater  Pressure  Development 


Water  pumping  out  of  the  ground  under  the  gradients  from  the  porewater 
pressures  generated  during  the  1983  Nihon-Kai-Chiba  earthquake  in  Japan  is 
shown  in  Figure  4.  The  speed  with  which  these  pressures  can  develop  and  their 
magnitude,  sometimes  results  in  an  explosive  ejection  of  the  soil  as  shown  by 
the  huge  crater  in  Figure  5.  Sometimes  the  pressures  are  relieved  in  a  more 
controlled  fashion  over  a  wide  area  giving  the  familiar  sand  volcano  effects 
shown  in  Figure  6. 


Settlement  of  the  ground,  similar  to  that  occurring  in  Figure  3(b)  is  shown  in 
Figure  7.  The  large  flow  deformations  that  can  occur  when  very  loose  sand 
liquefies  and  the  undrained  shearing  resistance  drops  to  low  values  of  residual 
strength,  are  exemplified  by  data  from  the  1995  Hyogoken-Nanbu  earthquake  in 
Kobe.  Figure  8  shows  the  large  lateral  movements  of  a  heavy  sea  wall  caisson 
that  occurred  when  the  foundation  soil  and  backfill  liquefied.  The  large 
settlements  are  shown  by  the  fleet  of  vans  trapped  in  a  water  filled  depression 
created  when  the  sea  wall  moved  outwards  up  to  3  m  (Figure  9). 

Large  quantities  of  sand  were  ejected  during  liquefaction  in  Kobe  at  level 
sites  such  as  the  parking  lot  of  the  amusement  park  on  Port  Island  shown  in 
Figure  10.  Widespread  settlement  occurred  all  over  Port  and  Rokko  Islands. 
Typical  settlement  is  shown  in  Figure  1 1  by  the  movement  of  the  ground  relative 
to  the  pile  foundation  of  the  elevated  railway  on  Port  Island.  These  examples 
stress  the  key  role  played  by  compaction  potential  in  causing  liquefaction. 


Mechanism  of  Porewater  Pressure  Development 


Potential  Indicators  of 
Liquefaction 


Characterization  of  In-Situ  State  of  Soil 

For  a  given  intensity  of  earthquake  shaking,  liquefaction  potential  is  con¬ 
trolled  by  the  tendency  for  compaction  and  the  capacity  for  drainage.  The  faster 
the  soil  can  drain,  the  lower  the  potential  for  high  porewater  pressures  and 
liquefaction  to  develop.  The  extent  of  potential  natural  drainage  during  an 
earthquake  is  difficult  to  evaluate  in  practice.  Therefore  a  reliable  indicator  of 
the  in-situ  permeability  would  be  a  major  improvement  in  practice.  We  will 
return  to  the  drainage  concept  later.  First,  we  will  look  at  indicators  of  other  in- 
situ  soil  characteristics. 


Compaction  Potential 


Cyclic  shearing  strain  is  recognized  as  the  prime  agent  for  generating  the 
volume  changes  in  saturated  cohesionless  soils  which  can  lead  to  liquefaction. 
The  small  cyclic  strains  resulting  from  earthquake  shaking  prior  to  liquefaction 
generally  result  in  the  compaction  of  both  loose  and  dense  sands  and  hence  in 
the  generation  of  porewater  pressures.  For  a  given  level  of  equivalent  cyclic 
shear  loading,  the  shear  strain  depends  on  the  in-situ  shear  modulus,  G.  There¬ 
fore,  for  a  given  intensity  of  loading,  the  shear  modulus  could  be  a  useful 
indicator  of  liquefaction  potential.  The  shear  modulus  in-situ  is  best  evaluated 
by  directly  measuring  the  in-situ  shear  wave  velocity,  W,  (G  =  where  p  is 
the  mass  density).  V,  is  measured  using  downhole  or  crosshole  methods,  or  by 
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Potential  Indicators  of  Liquefaction 


using  a  seismic  cone  (Robertson  et  al.  1986).  It  may  also  be  determined 
indirectly,  by  spectral  analysis  of  surface  waves  (Nazafian  and  Stokoe  1984). 

We  are  now  in  a  position  to  characterize  the  earthquake  shaking  at  a  site  in 
terms  of  earthquake  magnitude,  M,  and  cyclic  stress  ratio,  and  the  site 

condition  in  terms  of  shear  wave  velocity,  V^.  Since  is  dependent  on  the 
mean  normal  effective  stresses,  we  normalize  all  shear  wave  velocities  to  a 
standard  pressure,  P^,  of  the  1  tsf  (100  kPa),  giving. 


V.  =  V(P  la'  ) 

SI  s'^  a  VO' 


0.25 


(4) 


where  o',,;,-  effective  overburden  pressure  in  the  same  units  as  P^. 

Available  field  data  is  now  plotted  in  terms  of  the  state  parameters, 
and,  Vji,  and  a  boundary  curve  drawn  between  liquefied  and  non-Iiquefied  sites 
for  various  earthquakes  as  shown  in  Figure  12.  Liquefaction  potential  assess¬ 
ment  by  the  shear  wave  velocity  does  not  appear  to  be  very  dependent  on 
earthquake  magnitude.  This  is  surprising  because  porewater  pressure  develop¬ 
ment  depends  on  the  duration  as  well  as  the  intensity  of  shaking  and  magnitude 
is  a  good  indicator  of  duration.  Although  field  data  for  a  wide  range  of  shear 
wave  velocities  and  earthquake  magnitudes  is  not  yet  available,  shear  wave 
velocity  is  being  used  in  practice  as  an  indicator  of  liquefaction  for  materials 
which  are  difficult  to  probe  by  other  means  such  as  gravelly  or  cobbly  soils. 

The  liquefaction  assessment  chart  in  Figure  12  shows  the  essential  require¬ 
ments  of  any  indicator  of  liquefaction  potential;  it  must  be  capable  of  being 
plotted  on  a  chart  such  as  Figure  12  and  it  should  be  possible  to  refer  all  mea¬ 
surements  to  some  basic  reference  state.  Finally,  there  must  be  sufficient  direct 
field  data  to  develop  a  criterion  based  on  the  indicator  unless  there  is  a  strong, 
well-defined  correlation  between  the  proposed  indicator  and  another  indicator 
which  is  itself  supported  by  an  extensive  data  base. 


Potential  Indicators  of  Liquefaction 
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Bierschwale  and  Stokoe  (1984)  describe  some  of  the  early  procedures  for 
using  shear  wave  velocity  to  evaluate  liquefaction  potential.  Robertson  (1990) 
describes  the  more  recent  state-of-the-art. 

The  shear  wave  velocity  is  an  indicator  that  is  directly  linked  to  the  mecha¬ 
nism  by  which  porewater  pressure  is  developed  in  a  soil  by  seismic  loading.  It  is 
a  direct  measure  of  the  shear  modulus  which  controls  the  shear  strains  that  will 
develop  under  the  equivalent  cyclic  loading  described  by  Equation  1 ,  and  hence, 
controls  the  development  of  the  volumetric  strains  and  porewater  pressures. 

This  direct  connection  between  an  indicator  and  the  cause  of  liquefaction  is  the 
exception  rather  than  the  rule.  Other  in-situ  indicators  used  in  practice  measure 
the  combined  effects  of  a  number  of  in-situ  states  such  as  density,  stiffness, 
structure  and  fabric. 

An  important  requirement  for  liquefaction  is  a  very  high  degree  of  saturation. 
Otherwise,  some  compaction  can  be  accommodated  without  significant 
porewater  pressure  development.  It  may  be  possible  to  assess  the  degree  of 
saturation  using  P-waves  which  travel  at  different  velocities  through  the  soil  and 
water.  A  reliable  indicator  for  this  state  is  very  needed. 


Density 


The  shearing  resistance  of  sand  depends  on  the  relative  density,  that  is  the 
relationship  between  the  field  density  and  the  loosest  and  densest  states  of  the 
material.  Relative  density  is  defined  as, 

£)  —  ^  d(mm) 


(5) 


where  Yd  is  the  dry  density. 
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The  evaluation  of  the  loosest  and  densest  states  is  controlled  by  ASTM  Stan¬ 
dards  and  is  conducted  on  reconstituted  samples  of  the  soil.  Determining  the 
in-situ  density  is  a  difficult  problem.  In  some  critical  cases  in  which  expense 
was  justified,  frozen  samples  have  been  retrieved  (Byrne  et  al.  1994).  In  other 
cases,  in-situ  densities  have  been  determined  by  nuclear  probe.  This  procedure 
was  used  in  assessing  the  liquefaction  potential  of  Dashihe  Dam  in  China  by 
Morgenstem  and  Finn  (1994)  and  showed  promise.  There  is  a  lot  of  laboratory 
data,  chiefly  from  the  late  1960's  and  early  1970's  relating  liquefaction  potential 
of  reconstituted  sand  samples  from  cyclic  loading  tests  with  relative  density  but 
there  is  no  data  base  of  field  data  that  would  lead  directly  to  a  criterion  based  on 
relative  density.  However,  a  knowledge  of  the  in-situ  relative  density  would  be 
very  useful  to  an  experienced  engineer  in  making  judgements  about  liquefaction 
potential  and  its  consequences. 


Electrical  Resistivity,  Ground  Radar 


Any  in-situ  procedure  that  can  assess  the  fabric  and  structure  of  a  saturated 
deposit  has  the  capacity  to  be  an  indicator  of  liquefaction  potential.  Fabric  and 
structure  are  reflections  of  density,  shearing  stiffness,  shearing  resistance  and 
porosity.  Electrical  techniques  and  ground  radar  may  be  useful  tools  for  assess¬ 
ing  fabric  and  structure  and,  hence,  in  indicating  liquefaction.  Soil  fabric  and 
structure  are  important  for  other  reasons.  To  a  major  extent,  they  control 
permeability  which  controls  ease  of  drainage  and  hence  liquefaction  potential. 
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Methods  Used  in  Practice 


The  methods  used  in  current  practice  are  discussed  fully  in  the  companion 
paper  (Finn  1995)  referred  to  earlier.  They  are  listed  here  for  convenience  only. 
The  cited  paper  should  be  consulted  for  details  and  references. 


SPT-N  Method 


In  engineering  practice,  the  most  wisely  used  indicator  of  liquefaction 
potential  is  the  Standard  Penetration  Resistance  N  in  blows/ft  standardized  to 
1  tsf  and  60  percent  of  the  free  fall  energy  of  the  hammer  used  to  drive  the 
penetration  probe,  (NOeo-  The  largest  direct  data  base  on  liquefied  and  non- 
liquefied  sites  is  associated  with  (N,)6o. 


CPT:qci  Method 


In  recent  years,  the  cone  penetration  test  is  being  used  more  and  more  as  an 
indicator  of  liquefaction  because  of  the  reliability  and  reproducibility  of  mea¬ 
surements  of  cone  bearing,  q;.,  normalized  to  a  reference  pressure  of  1  tsf,  q^. 
There  is  a  good  correlation  between  N  and  q^  which  allows  the  data  base  for 
(N,)6o  to  support  q,i  as  an  indicator  of  liquefaction  potential  in  addition  to  the 
more  limited  field  data  directly  connecting  q,,  and  the  occurrence  of  liquefaction. 

Becker  Hammer  Test 


12 


Methods  Used  in  Practice 


This  test  is  used  to  establish  the  penetration  resistance  (NJ^o-  of  gravelly  and 
cobbly  soils  which  cannot  be  probed  by  the  cone  or  the  SPT  sampler.  A  large 
diameter  steel  tube  is  driven  by  blows  from  a  pile  driver.  These  blows  are 
calibrated  for  a  given  rig  and  hammer  to  (N,)6o  for  sands  by  running  Becker  and 
SPT  tests  in  sand,  side  by  side.  Liquefaction  potential  is  then  assessed  using  the 
liquefaction  assessment  chart  based  on  (N,)go. 


Shear  Wave  Velocity 


This  method  is  being  used  increasingly  in  the  same  soils  for  which  the  Becker 
test  was  developed.  In  these  soils,  the  shear  wave  velocity  is  derived  from  a 
spectral  analysis  of  surface  waves  (SASW).  This  method  requires  no  drill  holes 
and  no  penetration  testing.  It  has  become  a  fast  and  more  reliable  test  in  recent 
years. 


Methods  Used  in  Practice 
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Geophysical  Versus 
Geotechnical  Indicators 


The  geotechnical  approach  to  assessing  the  in-situ  state  of  a  potentially 
liquefiable  soil  is  very  location  specific.  The  state-of-the-soil  is  probed  at 
distinct  locations  and  a  great  many  soundings  are  necessary  in  order  to  define 
even  crudely  the  3-D  distribution  of  soil  states.  Furthermore,  the  methods 
require  drill  holes  or  the  driving  or  pushing  of  probes.  Geophysical  methods 
provide  data  on  a  much  larger  scale  and  fill  out  the  3-D  picture  much  better  and 
more  cheaply.  For  example,  the  details  of  soil  layering  can  be  defined  globally 
by  ground  penetrating  radar. 

The  geophysical  methods  such  as  ground  radar  and  the  spectral  analysis  of 
surface  waves  do  not  require  drill  holes  or  probes  and  theoretically  they  can 
function  in  any  type  of  soil.  This  feature  can  make  them  very  cost-effective. 
The  role  of  ground  radar  and  other  geophysical  (or  electrical)  measures  of  soil 
type  or  soil  state  in  detecting  liquefiable  soils  is  worthy  of  research  because  of 
the  potential  savings  in  costs  and  perhaps  an  increase  in  reliability  of  estimating 
liquefaction.  This  research  would  be  best  conducted  jointly  by  professionals 
from  both  groups. 
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Conclusions 


Any  in-situ  test  that  provides  a  measure  of  density,  shear  modulus,  fabric 
and/or  soil  structure  has  the  potential  to  be  a  useful  indicator  of  liquefaction 
potential.  However,  to  be  a  clear  discriminant  between  liquefiable  and  non- 
liquefiable  sites,  it  must  be  quantitatively  linked  to  site  performance  in  past 
earthquakes,  either  by  direct  measurement  or  through  a  reliable  correlation  with 
an  indicator  supported  by  such  a  data  base.  The  engineering  profession  is 
strongly  attached  to  the  indicators  now  being  used  to  assess  liquefaction.  One 
area  where  a  new,  reliable,  simple  indicator  would  be  very  welcome  is  gravely, 
cobbly  soils.  Existing  methods  are  effective  here  but  their  implementation  is  not 
always  easy,  and  interpretation  is  still  occasionally  controversial. 

There  is  scope  for  the  development  for  geophysical  (and  electrical)  methods 
for  the  assessment  of  liquefaction  potential.  The  attraction  of  such  methods  is 
that  they  usually  do  not  require  drill  holes,  are  more  global  in  their  assessments 
and  can  be  very  cost-effective.  Indications  of  potential  methods  were  presented 
in  the  paper,  but  the  main  function  of  the  paper  was  presentation  of  the  funda¬ 
mentals  of  liquefaction  to  provide  a  framework  of  understanding  that  would 
facilitate  the  development  of  new  indicators  for  engineering  practice. 


Conclusions 
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Figure  2.  Behavior  of  saturated  sands  under  static  loading  in  a  triaxial  test 
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Figure  6.  Sand  volcanoes  caused  by  ejection  of  liquefied  sand  by  seismic 
porewater  pressures 


Figure  7.  Settlement  of  ground  due  to  widespread  liquefaction 


Figure  10.  Widespread  ejection  of  sand  on  Port  Island  due  to  liquefaction 


Figure  1 1 .  Typical  settlement  of  level  ground  on  Port  Island  due  to  liquefaction 


Figure  12.  Liquefaction  assessment  chart  based  on  normalized  shear  wave  velocity 
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Introduction 


One  of  the  more  significant  factors  leading  to  ground  failure  during  earth¬ 
quakes  is  the  liquefaction  of  loose-  to  medium-dense  sands  below  the  water 
table.  Attention  was  focused  on  this  problem  for  the  first  time  as  a  result  of  the 
widespread  ground  failures  during  the  1964  earthquake  in  Niigata,  Japan 
(Kawasumi  1968).  Most  of  the  damage  in  Niigata  attributable  directly  to  the 
earthquake  was  associated  with  liquefaction  and  such  damage  has  been  a  signif¬ 
icant  factor  in  most  major  earthquakes  since  then.  Some  recent  examples  are 
the  extensive  damage  caused  by  liquefaction  in  the  Marina  District  in  San 
Francisco  during  the  Loma  Prieta  earthquake  of  1989  (Seed  et  al.  1990)  and 
the  destruction  of  port  facilities  on  Port  and  Rokko  Islands  in  Osaka  Bay  off 
Kobe  City,  during  the  1995  Hyogo-ken-Nanbu  earthquake. 

The  mechanics  of  liquefaction  are  now  well  understood  and  the  potential  for 
occurrence  can  be  estimated  with  a  reasonable  degree  of  confidence.  During 
shaking  the  sand  tends  to  compact.  The  water  in  the  pores  cannot  escape 
quickly  enough,  at  least  in  the  finer  sands,  to  accommodate  the  compaction 
instantaneously.  Therefore,  stresses  are  thrown  on  the  water  which  increase 
the  porewater  pressure  and  reduce  the  effective  or  intergranular  stresses 
between  the  sand  particles.  Sand,  a  frictional  material,  depends  on  the  effec¬ 
tive  stresses  between  the  grains  to  mobilize  shear  strength  and  resistance  to 
displacement. 


The  resistance  to  shearing  strain  or  deformation  is  also  reduced  by  increas¬ 
ing  porewater  pressure.  The  shear  modulus  G  which  controls  the  shearing 
strains  is  a  function  of  the  mean  normal  effective  stress.  Therefore,  as  the 
seismic  porewater  pressures  increase  the  shear  modulus  decreases.  In  the 
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extreme  case,  such  as  in  Niigata,  where  nearly  all  shear  strength  and  shear 
stiffness  was  lost  at  some  sites,  the  sand  behaved  like  a  liquid,  with  disastrous 
consequences  for  structures. 

Figure  1  illustrates  the  response  of  granular  soils  to  cyclic  loading  in  triaxial 
tests.  The  test  specimens  were  prepared  to  relative  densities  that  covered  the 
range  of  sand  states  from  highly  contractive  to  highly  dilative  under  shearing 
deformation  (Vaid  and  Chem  1985).  All  test  specimens  carry  initial  static 
shear  stresses,  t^j.  This  simulates  the  most  common  simation  in  earth  struc¬ 
tures  and  under  foundations.  Contractive  stress-strain  response  is  shown  in 
Figure  la.  At  some  stage  in  the  cyclic  loading,  strain  softening  occurs  leading 
to  a  large  drop  in  undrained  strength  which  then  remains  constant  at  the 
undrained  residual  strength  of  the  soil  over  a  large  range  in  strain. 

Large  deformations  occur  when  is  greater  than  the  steady  state  strength 
Tjs  (or  Sus).  Deformations  may  also  be  significant  when  tj,  <  if  Xj,  +  x^y  > 
Xjj  where  is  the  cyclic  shear  stress.  The  phenomenon  of  large  deformations 
occurring  at  constant  undrained  strength  is  termed  liquefaction. 

At  the  other  extreme  the  sand  is  strongly  dilative  and  deformation  is  the 
result  of  progressive  stiffness  degradation  in  the  sand  due  to  the  build  up  in 
porewater  pressure  with  cycles  of  loading.  This  phenomenon  is  called  cyclic 
mobility  and  is  illustrated  in  Figure  Ic. 

Between  these  two  clearly  defined  states  is  a  transitional  region  combining 
features  of  both  contractive  and  dilative  response.  This  type  of  behavior  is 
characterized  by  a  loss  in  undrained  strength  due  to  strain  softening  over  a  lim¬ 
ited  range  followed  by  a  regain  of  undrained  strength.  Further  deformations 
result  from  cyclic  mobility.  This  behavior  is  illustrated  in  Figure  lb.  The 
deformation  patterns  of  these  sands  are  shown  together  in  Figure  1(d). 

A  flow  chart  for  assessment  of  the  effects  of  various  cyclic  loading  condi¬ 
tions  is  given  in  Figure  2. 
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The  magnitude  of  potential  deformations  at  a  liquefied  site  with  contractive 
behavior  depends  on  the  extent  to  which  the  driving  shear  stresses  exceed  the 
residual  strength.  The  residual  strength  therefore  is  a  key  parameter  control¬ 
ling  the  extent  and  cost  of  remedial  measures  to  limit  deformations  to  tolerable 
levels. 

In  engineering  practice  the  word  liquefaction  is  usually  associated  with  the 
development  of  high  porewater  pressures  due  to  seismic  loading  irrespective  of 
whether  large  deformations  occur  or  not.  Indeed  the  empirical  procedures  for 
assessing  liquefaction  potential  are  all  based  primarily  on  surface  manifesta¬ 
tions  of  high  porewater  pressure  such  as  sand  boils  without  consideration  of 
whether  the  soil  would  behave  in  a  contractive  or  dilative  mode  during 
shearing. 

In  this  paper,  procedures  will  be  presented  for  determining  the  potential  for 
triggering  liquefaction  and  for  estimating  residual  strength. 
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Empirical  Assessment  of 
Liquefaction 


The  potential  for  triggering  liquefaction  is  evaluated  using  procedures  based 
on  the  standard  penetration  test,  cone  tip  penetration  resistance  and  normal 
shear  wave  velocity.  The  procedures  involves  three  steps: 

a.  Characterizing  the  dynamic  effects  of  the  earthquake. 

b.  Characterizing  the  in  situ  state  of  the  soil. 

c.  Application  of  a  criterion  for  the  incidence  of  soil  liquefaction. 


Characterizing  the  Earthquake 


Seismic  shear  stresses  play  a  major  role  in  the  development  of  liquefaction. 
Time  histories  of  shear  stresses  are  usually  very  non-uniform  and  are  difficult 
to  apply  in  empirical  methods.  Seed  (1983)  suggested  replacing  the  irregular 
time  history  by  an  equivalent  number  of  uniform  cycles  and  normalizing  the 
shear  stresses  by  dividing  by  the  effective  overburden  stress.  Seed  proposed 
that  this  uniform  cyclic  shear  stress  ratio  (CSR)  be  determined  by: 


CSR 


O  .  Oi  u 

0.65  ^  r 


/Ij. 


A> 
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where 


T(.  =  equivalent  average  shear  stress 
Oy  =  total  overburden  stress 
a'yo  =  effective  overburden  stress 

~  maximum  surface  acceleration  in  units  of  g 
g  =  acceleration  due  to  gravity 

rj  =  a  reduction  factor  to  account  for  soil  flexibility  and  depth 

In  Japanese  practice  the  reduction  factor  r^,  is  often  approximated  using  the 
following  expression: 

/■^  =  (1  -  0.0152) 


where  Z  =  depth  in  meters. 


Characterizing  the  In  Situ  State  of  the  Soil 

The  in  situ  state  of  the  soil  can  be  characterized  by  in  situ  testing  or  by  test¬ 
ing-  undisturbed  samples  obtained  from  boreholes  or  test  pits.  Granular  soils 
are  extremely  difficult  to  sample  without  disturbance.  Hence,  methods  to  char¬ 
acterize  the  soil  in  the  ground  rely  heavily  on  in  situ  tests.  The  most  widely 
used  methods  are  the  standard  penetration  test  (SPT)  and  the  cone  penetration 
test  with  pore  pressure  measurement  (CPTU). 


Standard  Penetration  Test  (SPT) 

Seed  and  his  colleagues  (Seed  1979,  Seed  et  al.  1983,  and  Seed  et  al.  1985) 
developed  correlations  between  the  SPT  N-value  and  the  cyclic  stress  ratio  to 
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cause  liquefaction  (CSR<.rit)  during  earthquakes  of  magnitude  M  =  7.5.  The 
correlations  which  are  presented  in  Figure  3,  were  based  on  the  observed 
response  of  sites  during  earthquake  loading.  Sites  were  considered  to  have 
liquefied  based  on  observed  surface  features,  such  as  sand  boils.  Lower  bound 
curves  separating  liquefied  from  non-liquefied  sites  are  shown  in  Figure  3  cor¬ 
responding  to  various  fines  contents  of  the  sands. 

To  compare  the  ground  conditions  at  one  site  with  those  of  another,  it  is 
necessary  to  standardize  the  measured  penetration  values  to  a  standard  driving 
energy  and  effective  overburden  pressure.  Seed  normalized  the  SPT  to  an 
energy  level  of  60  percent  of  the  free-fall  potential  energy  of  the  hammer  and 
an  effective  overburden  pressure  of  100  kPa  (1  tsf).  Hence,  the  correlations 
presented  in  Figure  3  show  the  normalized  SPT  N  value,  (N,)6o. 

The  correction  for  energy  level  may  be  made  using  the  data  in  Table  1  or  by 
direct  measurement  of  the  energy  imparted  to  the  rods. 

The  correction  factors  Cn  for  overburden  pressure  may  be  determined  using 
the  curves  in  Figure  4  which  allow  for  the  effects  of  relative  density  or  by 
using  the  equation  (Liao  and  Whitman  1985), 


(3) 


where  o\o  is  expressed  in  tons/per/sq  ft. 

The  use  of  equation  (3)  which  ignores  any  effects  of  relative  density  is 
becoming  common.  Judgment  is  required  in  limiting  the  maximum  value  of  the 
correction  factor.  Generally  the  upper  limit  is  taken  as  2.0. 
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Table  1 


Summary  of  Energy  Ratios  for  SPT  Procedures  (Seed  et  al.  1985) 


Hammer 

Type 


Est.  Rod  Correction  Factor  for 

Hammer  Release  Energy  (percent)  60  percent  Rod  Energy 


Free-Fall 
Rope  &  Pulley 
with  Special 
Throw  Release 


Rope  &  Pulley 
Rope  &  Pulley 


Rope  &  Pulley 


Free-Fall^ 
Rope  &  Pulley 


Japanese  SPT  results  have  additional  corrections  for  borehole  diameter  and  frequency 
effects. 

^  Prevalent  method  in  the  United  States  today. 

^Pilcon-type  hammers  develop  and  energy  ratio  of  about  60  percent. 


The  correlations  shown  in  Figure  3  are  representative  of  earthquakes  with  a 
magnitude  M  =  7.5.  The  critical  correlations  for  earthquakes  of  other  magni¬ 
tudes  may  be  established  by  multiplying  the  critical  cyclic  stress  ratios  by  the 
magnitude  dependent  correction  factors  in  Table  2  (Seed  and  Harder  1990). 


There  are  two  important  limitations  associated  with  Figure  3.  The  field  data 
correspond  to  level  ground  conditions  with  no  initial  static  shear  stresses  on 
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horizontal  planes  and  to  effective  overburden  pressures  less  than  150  kPa. 

Seed  (1983)  outlined  procedures  for  making  corrections  when  these  conditions 
are  violated. 

The  first  estimate  of  the  liquefaction  resistance  of  a  soil  element  in  a  dam  or 
slope  is  determined  using  the  in  situ  (Ni)6o  penetration  resistance  and  the  appro¬ 
priate  curve  for  critical  conditions  in  Seed’s  liquefaction  assessment  chart. 

This  resistance  must  then  be  corrected  for  deviations  from  the  standard  condi¬ 
tions  of  the  database  underlying  the  chart.  A  typical  element  in  a  slope,  for 
example,  will  carry  a  static  shear  stress  Tj,  on  the  horizontal  plane  and  therefore 
has  an  initial  shear  stress  ratio  =  «.  A  correction  factor  K„  is  estab¬ 
lished  for  various  values  of  «  by  laboratory  tests.  As  an  illustration,  the  K„ 
factor  for  Folsom  gravels  is  shown  in  Figure  5  (Hynes  et  al.  1988).  Note  that, 
as  is  commonly  assumed  in  practice,  initial  static  shear  increases  liquefaction 
resistance  substantially.  However,  this  increase  applies  to  resistance  to  cyclic 
mobility  rather  than  to  liquefaction,  that  is  to  non-contractive  materials. 

Vaid  and  Chem  (1985)  have  shown  for  contractive  materials  that  initial 
static  shear  decreases  the  resistance  to  cyclic  loading.  This  suggests  that  for 
contractive  soils,  the  cyclic  strength  should  be  reduced  not  increased  for  the 
effects  of  static  shear.  K„  factors  less  than  unity  have  been  reported  by  Mar- 
cuson  et  al.  (1990)  and  Seed  and  Harder  (1990).  Typical  values  of  K„  may  be 
found  in  Figure  6  for  a  range  in  relative  densities. 

For  effective  overburden  pressures  other  than  in  the  range  100-150  kPa,  a 
correction  factor  is  used.  Values  of  K,,  for  many  soils  are  shown  in  Fig¬ 
ure  7.  Increasing  confining  pressure  can  lead  to  a  substantial  reduction  in 
resistance  to  cyclic  loading  (Seed  and  Harder  1990).  However,  Vaid  and 
Thomas  (1994)  have  shown  (Figure  8)  that  for  some  sands  the  reductions  are 
much  smaller  than  those  given  by  the  generalized  average  reduction  curve  in 
Figure  7  which  is  commonly  used  in  practice.  data  derived  from  tests  on 
sands  from  the  foundations  of  Duncan  Dam  confirmed  that  the  average  curve 
in  Figure  7  may  be  very  conservative  in  some  cases  (Byrne  et  al.  1993).  These 
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findings  suggest  that  for  important  projects,  site  specific  values  of  K^,  and  K 
should  be  determined. 


Cone  Penetration  Test  (CPI) 


Electric  cone  tip  penetration  resistance  (Q^)  has  been  measured  at  few  sites 
as  yet  where  the  occurrence  or  nonoccurrence  of  liquefaction  during  actual 
earthquakes  have  been  documented.  Therefore  the  data  base  on  cone  penetra¬ 
tion  resistances  at  liquefied  sites  has  been  extended  by  converting  the  well  doc¬ 
umented  correlation  between  liquefaction  potential  and  standard  penetration 
resistance  to  cone  penetration  resistance. 

There  has  been  a  number  of  studies  on  the  correlation  between  cone  resis¬ 
tance  and  standard  penetration  resistance  (Douglas  et  al.  1981,  Robertson  et  al. 
1983).  Figure  9  indicates  that  the  ratio  of  cone  resistance  to  blowcount 
increases  with  increasing  mean  grain  diameter  (Robertson  et  al.  1983). 

The  liquefaction  assessment  chart  in  Figure  3  based  on  (Ni)^,)  can  now  be 
converted  to  a  chart  based  on  Q^.,  as  shown  in  Figure  10,  where  Q^j  is  the  tip 
penetration  resistance  normalized  to  100  kPa  (1  tsf)  using  either  Figure  4  or 
Equation  3.  Cyclic  strength  must  be  corrected  for  initial  static  shear  stress  (KJ 
and  overburden  pressure  (KJ  in  the  same  manner  as  (NJ^p. 


Fines  Content 


Seed  et  al.  (1985)  reviewed  sites  that  did  and  did  not  liquefy  during  earth¬ 
quakes  where  the  fines  content  was  greater  than  5  percent.  They  found  that  for 
the  same  penetration  resistance  the  liquefaction  resistance  increased  with 
increasing  fines  content  (Figures  3  and  10). 
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The  field  data  in  Figure  3  shows  a  strong  dependency  of  liquefaction  resis¬ 
tance  on  fines  content.  But,  surprisingly,  there  appears  to  be  considerable  con- 
fiision  about  the  effect  of  fines  on  liquefaction  resistance.  Part  of  the  problem 
depends  on  the  manner  on  which  the  comparison  is  made.  In  the  case  of  the 
field  data,  for  the  same  normalized  penetration  resistance  (Ni)go  the  liquefac¬ 
tion  resistance  increases  with  increasing  fines  content.  What  the  soils  being 
compared  have  in  common  is  the  (Nj)^,  value.  If  the  basis  for  comparison  is 
changed  the  conclusion  may  well  be  different.  For  example,  Troncoso  (1990) 
compared  the  cyclic  strength  of  tailings  sands  with  different  silt  contents  rang¬ 
ing  from  0  to  30  percent  at  a  constant  void  ratio.  On  this  basis  of  comparison 
he  found  that  cyclic  strength  decreased  with  increased  silt  content  (Figure  11). 
Kuerbis  and  Vaid  (1989)  studied  the  effects  of  fines  using  samples  deposited  in 
a  manner  that  replicated  field  deposition  conditions.  They  found  for  the  sand 
tested  that  up  to  20  percent  fines  could  be  accommodated  within  the  sand  skele¬ 
ton  and  that  samples  with  the  same  sand  skeleton  void  ratio  had  the  same  cyclic 
strength  for  fines  content  less  than  20  percent. 

The  different  methods  for  considering  the  effects  of  fines  are  summarized  in 
Table  3. 


Table  3 

Effects  of  Fines  on  Cyclic  Strength 

Criterion  of  Equivalency  Clean  and  Silty  Sand 

Effect  of  Fines  on  Cyclic  Strength 

1.  Same  (N,  Iso 

Increase 

2.  Same  void  ratio 

Decrease 

3.  Same  void  ratio  in  sand  skeleton 

No  change  while  fines  content  can  be 

accommodated  in  sand  voids 

Another  way  of  looking  at  the  field  data  is  to  note  that  the  penetration  resis¬ 
tance  of  silty  sands  is  somewhat  lower  than  for  clean  sands  at  the  same  cyclic 
strength.  Corrections  based  on  the  curves  in  Figure  3,  can  be  applied  to  the 
(N,)«,  values  measured  in  silty  sands  and  then  the  curve  for  clean  sands  may  be 
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used  ,0  determine  cyclic  strenstt.  for  the  corrected  (N,)„  value.  Such  correc- 
tions  have  been  widely  used. 


Corrections  to  the  normalized  cone  penetration  resistance  with  fines  (Q  ,) 
can  be  established  from  Figure  10.  These  may  be  used  to  adjust  cone  penlLl 
tion  resistance  data  from  test  sites  to  cone  penetration  resistances  of  clean  sands 
with  similar  liquefaction  potential.  The  correction  is  given  by 


AO 


(4) 


where  AQ„  are  given  in  Figure  12,  from  Robertson  (1990). 

It  is  possible  to  estimate  the  fines  content  (FC)  directly  from  CPT  and/or 
CPTU  (CPT  with  porewater  pressure  measurements)  data.  Figure  13  presents 
the  latest  soil  classification  chart  based  on  normalized  CPT  and  CPTU  data. 
Soils  that  fall  in  zone  6  are  generally  clean  sands  or  silty  sands  with  a  small 
amount  of  fines.  Soils  that  fall  in  zone  5  are  silty  sands  and  sandy  silts  that 
generally  have  fines  contents  greater  than  about  15  percent. 

Based  on  data  collected  in  the  Richmond  area  of  the  Fraser  Delta  of  British 
Columbia  a  correlation  has  been  established  by  Finn  et  al.  (1989)  between  fines 
content  and  the  time  for  50  percent  dissipation  (T^^)  of  the  pore  pressures 
developed  during  penetration  (Figure  14).  The  results  in  Figure  14  suggest 
that  for  Tso  >  50  sec  the  fines  content  is  greater  than  40  percent.  For  T^q 
between  10  sec  and  50  sec  the  cone  penetration  process  is  partially  drained  and 
there  is  a  poor  correlation  between  T^o  and  fines  content. 

If  the  fines  are  plastic,  criteria  based  on  field  data  from  Chinese  earthquakes 
must  be  taken  into  account  (Wang  1979)  to  determine  if  the  soil  will  liquefy. 
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Liquefaction  Potential  of  Plastic 
Fine-Grained  Soils 


Liquefaction  potential  of  plastic  fine-grained  soils  are  determined  using  the 
Chinese  criteria  developed  by  Wang  (1979).  These  criteria  are: 

a.  Percent  finer  than  0.005  mm  15  percent. 

b.  Liquid  limit,  LL  ^  35  percent. 

c.  Natural  water  content  ^  0.9  LL. 

d.  Liquidity  index,  ^  0.75. 

Soils  which  satisfy  all  four  criteria  are  judged  vulnerable  to  liquefaction  or 
significant  strength  loss.  In  addition,  any  fine  grained  soils  for  which  the  stan¬ 
dard  penetration  resistance  N  4  are  considered  to  be  vulnerable  to  liquefac¬ 
tion  or  to  significant  strength  loss  under  cyclic  loading  whether  or  not  they 
satisfy  the  Chinese  criteria  (Woodward  Clyde  Consultants  1989).  The  Chinese 
criteria  are  usually  applied  strictly  with  no  account  taken  of  uncertainties  in  the 
measurements  of  the  parameters  in  the  criteria. 

In  coimection  with  investigations  of  the  liquefaction  potential  of  Sardis  Dam 
in  Mississippi  (Firm  et  al.  1990b),  Woodward  Clyde  Consultants  (1989)  sug¬ 
gested  that  allowances  should  be  made  for  uncertainties  in  the  measured  values 
of  the  parameters  in  the  criteria.  Following  a  study  of  Wang’s  database,  they 
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recommended  ignoring  the  liquidity  index  and  malting  the  following  changes  in 
the  measured  soil  properties  before  applying  the  criteria: 

a.  Decrease  per  cent  fines  by  15  percent. 

b.  Decrease  LL  by  5  percent. 

c.  Increase  water  content  by  3  percent. 

These  changes  increased  significantly  the  extent  of  the  soils  vulnerable  to 
liquefaction,  and  strength  loss  so  that  almost  the  entire  length  of  dam  required 
remediation.  Therefore,  engineers  of  the  Vicksburg  District  of  the  U.S.  Army 
Corps  of  Engineers,  who  are  responsible  for  the  safety  assessment  and  remedi¬ 
ation  of  Sardis  Dam,  reviewed  reports  on  the  scatter  in  measured  index 
properties  in  U.S.  Corps  of  Engineers’  laboratories  over  the  last  30  years  to 
determine  the  likely  range  of  variations  in  test  data.  In  addition  they  conducted 
tests  on  samples  of  standard  soils  of  low  to  medium  plasticity  in  their  own  labo¬ 
ratory  to  establish  the  scatter  in  their  own  data.  These  standard  soils  are  used 
to  check  comparability  of  testing  procedures  and  results  between  different 
Corps  of  Engineers,  laboratories  and  private  laboratories.  As  a  result  of  these 
studies  the  following  changes  in  measured  properties  were  adopted  before 
applying  the  Chinese  criteria  (again  ignoring  the  liquidity  index): 

a.  Decrease  the  fines  content  by  5  percent. 

b.  Decrease  the  liquid  limit  by  2  percent. 

c.  Increase  the  water  content  by  2  percent. 

This  change  reduced  the  length  requiring  remediation  to  about  1220  m 
(4000  ft). 

The  impact  of  the  Chinese  criteria  on  the  extent  of  remediation  necessary 
for  stability  appeared  to  be  so  critical  that  an  investigation  of  Chinese  proce- 
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dures  was  undertaken  by  Koester  (1990).  The  Chinese  determine  the  liquid 
limit  using  a  fall  cone  rather  than  the  Casagrande  device  generally  used  in 
North  America.  Using  a  standard  Chinese  fall  cone  and  following  Chinese 
standard  SD  128-007-84,  Koester  (1990)  showed  that  the  fall  cone  gives  a 
liquid  limit  about  3  percent  to  4  percent  greater  than  the  Casagrande  device. 

He  continued  the  study  of  the  Chinese  criteria  for  liquefaction  and  reported  his 
final  conclusions  in  Koester  (1992). 

On  the  basis  of  all  the  above  studies  the  following  changes  in  measured 
index  properties  were  finally  adopted  to  account  for  uncertainty  before 
application  of  the  Chinese  criteria: 

a.  Decrease  the  fines  content  by  5  percent. 

b.  Increase  the  liquid  limit  by  1  percent. 

c.  Increase  the  water  content  by  2  percent. 

These  changes  reduced  the  length  requiring  remediation  to  about  926  m 
(3000  ft). 

Clearly  the  Chinese  criteria  for  evaluating  the  potential  for  liquefaction  or 
significant  strength  loss  in  clayey  soils,  based  on  liquid  limit,  water  content  and 
percent  fines  <  0.005  mm,  can  have  a  major  impact  on  the  extent  of  remedial 
measures  necessary  to  achieve  stability  in  earth  structures  with  potentially 
liquefiable  fine  grained  materials. 

Before  applying  the  Chinese  criteria  the  uncertainties  in  the  measured  soil 
properties  should  be  taken  into  account  in  a  reasonably  conservative  manner. 
This  may  be  done  by  adjusting  the  measured  water  content,  liquid  limit,  and 
the  fines  content  before  applying  the  criteria.  The  amount  of  these  adjustments 
should  be  based  preferably  on  the  estimated  variability  in  data  appropriate  for 
the  laboratory  conducting  the  tests.  In  the  absence  of  such  specific  informa¬ 
tion,  the  adjustments  noted  above  of  increasing  the  liquid  limit  by  1  percent. 
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the  water  content  by  2  percent,  and  decreasing  the  fines  content  by  5  percent 
may  be  considered.  These  adjustments  reflect  conservative  estimates  of  the 
variability  to  be  expected  from  very  experienced  persormel  operating  under 
high  standards  of  quality  control. 
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Liquefaction  Assessment  from 
Shear  Wave  Velocity 


Empirical  methods  have  been  developed  to  evaluate  liquefaction  resistance 
directly  from  shear  wave  velocity  (Bierschwale  and  Stokoe  1984).  Over  the 
past  15  years,  significant  advances  have  been  made  in  measuring  shear  wave 
velocities  in  the  field.  Accurate  and  detailed  profiles  can  be  determined  with 
conventional  crosshole  and  downhole  seismic  methods  (Stokoe  and  Hoar  1978, 
Woods  1978).  Shear  wave  velocity,  Vj  is  influenced  by  many  of  the  variables 
that  influence  liquefaction  potential,  such  as,  soil  density,  confinement,  stress 
history  and  geologic  age.  Thus,  Vj,  has  promise  as  a  field  index  in  evaluating 
liquefaction  susceptibility.  The  most  significant  advance  in  recent  years  has 
been  the  development  of  the  seismic  cone  penetration  test  (Robertson  et  al. 
1986)  which  allows  the  measurement  of  seismic  velocity  in  conjunction  with 
penetration  resistance. 

The  major  advantage  of  using  shear  wave  velocity  as  an  index  of  liquefac¬ 
tion  resistance  is  that  it  can  be  measured  in  soils  that  are  hard  to  sample  such  as 
cohesionless  silts  and  sands  or  hard  to  penetrate  such  as  gravels. 


Direct  Shear  Wave  Velocity  Correlations 


The  limiting  shear  wave  velocities  separating  liquefied  from  non-liquefied 
sites  for  a  given  intensity  and  duration  of  shaking  must  be  determined  from 
field  data.  So  far  the  data  base  is  quite  limited  but  it  clearly  shows  that  shear 
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wave  velocity  may  be  a  useful  index  of  liquefaction  potential.  Data  from  sites 
in  the  Imperial  Valley,  California  which  liquefied  during  the  1979  Imperial 
Valley,  1981  Westmorland  and  1987  Superstition  Hills  earthquakes  suggest  that 
the  limiting  shear  wave  velocity  separating  liquefiable  from  nonliquefiable  sites 
is  about  140  m/s  for  earthquakes  of  local  magnitude  =  6.5  generating  peak 
ground  accelerations  of  about  0.17  g  (Holzer  et  al.  1988,  Youd  and  Wieczorak 
1984). 


Shear  wave  velocity  is  a  function  of  void  ratio  and  effective  confining 
stress.  Hence,  for  a  sand  of  constant  void  ratio  (constant  density)  the  shear 
wave  velocity  will  increase  with  increasing  depth.  Therefore,  a  correlation 
between  V^  and  CSR  to  cause  liquefaction  should  be  based  on  normalized 
effective  overburden  stress,  similar  to  the  manner  in  which  penetration  resis¬ 
tance  is  normalized  with  depth.  Shear  wave  velocity  being  proportional  to  the 
square  root  of  the  shear  modulus  is  a  function  of  the  effective  overburden 
stress  to  the  power  of  0.25, 


(5) 


Therefore,  a  normalized  shear  wave  velocity  can  be  established  using  the 
relationship 


^eL  -  yoipj^'Ar- 


(6) 


where 


Pa  =  reference  stress,  typically  100  kPa 
a'vo  =  effective  overburden  stress  in  same  units  as 

The  proposed  correlation  between  normalized  shear  wave  velocity  and  the 
critical  cyclic  stress  ratio  necessary  to  prevent  liquefaction  is  shown  in 
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Figure  15  (Robertson  1990).  The  data  base  is  still  very  small  so  it  is  still  too 
soon  to  use  the  shear  wave  velocity  criterion  alone  in  practice. 
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Liquefaction  Potential  of 

Gravels  by  Becker  Penetration 
Test  (BPT) 


The  presence  of  large  gravel  and  cobble  particles  preclude  the  use  of  the 
SPT  in  the  evaluation  of  liquefaction  potential.  Misleadingly  high  SPT  blow- 
counts  would  be  recorded  as  the  2  in  sampler  bounced  off  the  large  particles. 

A  much  larger  and  heavier  penetration  tool  is  required  to  penetrate  gravels  and 
provide  a  continuous  record  of  penetration  resistance.  Harder  and  Seed  (1986) 
developed  the  Becker  Penetration  Test  (BPI)  for  this  situation. 


In  the  BPT  a  double-walled  casing,  generally  6.6  in.  O.D.  in  the  U.S.,  is 
driven  into  the  ground  with  a  diesel  pile  hammer  (Figure  16).  The  casing  can 
be  driven  with  an  open  bit  and  reversed  air  circulation  to  obtain  disturbed  sam¬ 
ples  or  with  a  plugged  bit  and  driven  as  a  solid  penetrometer.  The  blows 
required  for  each  foot  of  penetration  are  recorded  and  provide  a  record  of 
penetration  resistance  for  the  entire  profile.  Two  types  of  rigs  are  in  use,  the 
B-180  rig  (also  known  as  HAV-180)  and  the  AP-1000  developed  later  in  the 
1970' s.  Both  rigs  use  the  same  model  of  diesel  hammer,  an  ICE  Model  180. 
Harder  and  Seed  (1986)  showed  that  the  older  B-180  rig  was  more  efficient 
with  the  result  that  the  blowcounts  from  AP-1000  rig  were  about  50  percent 
higher  than  those  from  the  B-180  rig.  Harder  and  Seed  (1986)  developed  the 
BPT  procedures  using  the  AP-1000  so  blowcounts  obtained  by  the  B-180  rig 
must  be  adjusted  to  equivalent  AP-1000  blowcounts  unless  correlations  are 
developed  specifically  for  the  B-180  rig.  The  adjustment  factor  is  1.5. 
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The  first  step  in  evaluating  liquefaction  potential  using  BPT  blowcounts 
is  to  correct  the  blowcounts  to  standard  combustion  conditions,  using  the 
procedures  prescribed  by  Harder  and  Seed  (1986).  Then  using  a  correlation 
developed  by  Harder  and  Seed  (1986)  and  shown  in  Figure  17,  the  Nt,.  are  con¬ 
verted  to  equivalent  SPT  N^q.  Finally,  the  SPT  blowcounts  are  converted  to 
(N,)6o  for  use  in  the  liquefaction  assessment  chart  in  Figure  3  or  Figure  5. 

The  BPT  procedure  has  been  used  to  evaluate  the  liquefaction  potential  of 
foundation  gravels  in  the  Mormon  Island  Dam  near  Sacramento,  California. 

An  excellent  description  of  the  use  of  the  BPT  on  this  project  may  be  found  in 
Hynes  et  al,  (1988). 

The  correlation  between  Nb^  and  Nbo  is  established  by  conducting  both  tests 
in  the  same  sand  profile.  Is  this  correlation  still  valid  when  BPT  tests  are  run 
in  deposits  of  gravel  and  cobbles?  Even  in  relatively  fine-grained  soils,  the 
correlation  between  SPT  and  CPT  has  been  found  to  be  dependent  on  mean 
grain  size  as  shown  in  Figure  9.  A  larger  database  on  the  liquefaction  of 
gravels  is  needed  to  resolve  this  question. 

As  the  use  of  the  Becker  penetration  test  spread,  two  major  problems 
emerged.  The  Harder-Seed  (1986)  correlation  appeared  not  to  be  as  generally 
applicable  as  originally  thought  and  the  effects  of  casing  friction  on  blowcounts 
at  deep  penetration  became  a  major  concern.  Fundamental  research  on  both  of 
these  questions  have  been  conducted  in  British  Columbia,  Canada  by  Sy 
(1993a)  using  different  Becker  rigs  and  a  variety  of  sites. 


Energy  Standard  for  the  Becker  Penetration  Test 


Sy  (1993a)  has  instrumented  Becker  penetration  tests  to  measure  the  energy 
transferred  through  (ENTHRU)  the  anvil  and  capblock  to  the  Becker  casing. 
He  showed  that  the  average  transferred  energy  was  about  30  percent  of  the 
rated  energy  of  the  hammer  for  the  various  types  of  rigs.  Sy  (1993a)  proposed 
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adopting  this  energy  level  as  the  standard  for  establishing  the  Becker 
blowcount/ft,  Ni,3o. 

N^(ENTHRU/30) 


Using  this  standard,  he  has  been  able  to  show  that,  in  the  same  deposit,  the 
HAV-180  and  the  AP-1000  rigs  deliver  the  same  Nb3o. 


Blowcount  data  for  HAV-180  and  AP-1000  rigs  are  shown  in  Figure  18a 
from  tests  BPT-2  and  BPT-5,  respectively,  at  the  same  site.  The  corrected 
blowcounts  based  on  the  Harder-Seed  correction  procedure  are  shown  in  Fig¬ 
ure  18b.  Note  that  the  blowcounts  from  the  more  energy  efficient  HAV-180 
rig  are  less  than  the  blowcounts  from  the  AP-1000  rig  even  after  the  correc¬ 
tion.  However,  when  both  sets  of  blowcount  data  are  corrected  to  the  same 
transferred  energy  level  of  30  percent  of  the  rated  hammer  energy,  both  sets  of 
data  plot  together  (Figure  18c). 


Effect  of  Friction  on  Becker  Blowcount 


The  friction  between  the  Becker  casing  and  the  ground  increases  with  depth 
of  penetration.  This  friction  is  not  present  in  the  Standard  Penetration  Test. 
Clearly,  it  is  not  possible  to  maintain  a  stable  correlation  between  BPT  and 
SPT  soundings.  The  original  Harder-Seed  (1986)  correlation  includes  the 
effects  of  friction  for  depths  up  to  50  ft  only.  One  solution  to  the  friction  prob¬ 
lem  may  be  to  use  different  correlations  over  different  depth  ranges. 

Tests  were  conducted  by  Sy  (1993a)  at  the  site  of  a  large  development  in 
Richmond,  B.C.,  Canada,  to  determine  the  effects  of  casing  friction  on  Becker 
blowcounts.  The  BPT  was  conducted  in  two  ways.  First  the  test  was  con¬ 
ducted  in  the  standard  way  using  an  uncased  hole.  Then  a  BPT  was  conducted 
using  modified  procedures.  The  bottom  of  the  casing  was  the  standard  size 
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with  a  6.6  in.  O.D.  over  a  length  of  1  ft.  Above  this,  the  casing  O.D.  was 
5.5  in.  During  penetration  drilling  mud  was  pumped  into  the  casing,  exited 
through  holes  near  the  bottom  of  the  5.5  in.  O.D.  casing  and  was  forced 
upwards  between  the  casing  and  the  ground.  This  procedure  ensured  virtually 
friction  free  driving  over  the  entire  penetration  depth.  SPT  tests  were  also 
conducted  at  the  site.  The  results  are  shown  in  Figure  19a  in  which  both  BPT 
and  SPT  energy  corrected  blowcounts  are  plotted  with  depth.  The  comparison 
between  the  two  BPT  blowcounts  are  for  the  same  level  of  transferred  energy. 
Therefore,  the  different  casing  frictions  are  primarily  responsible  for  the  differ¬ 
ences  in  penetration  resistances.  The  effect  of  friction  on  the  blowcounts  is 
very  substantial  but  difficult  to  see  clearly  because  of  the  small  blowcount 
scale.  The  result  is  much  clearer  in  Figure  19b  where  the  BPT  blowcounts  are 
normalized  by  the  SPT  values.  Clearly  the  effect  of  friction  is  an  order  of 
magnitude  and  increases  with  depth. 

Recently  at  the  Terzaghi  Dam  in  British  Columbia,  BC  Hydro  decided  to 
use  Becker  blowcounts  from  cased  holes,  although  the  data  were  recognized  to 
be  conservative,  rather  than  the  data  from  the  uncased  holes  which  incorpo¬ 
rated  frictional  effects  that  could  not  be  evaluated  at  the  time  (Sy  1993b). 

Sy  (1993a)  at  the  University  of  British  Columbia  has  developed  correlations 
between  SPT  and  BPT  Nb3o  for  various  average  friction  forces  on  the  cas¬ 
ing,  The  friction  was  evaluated  using  CAPWAP  analyses  (Goble  et  al.  1980). 
Therefore,  by  using  the  energy  correction  approach,  it  is  possible  to  take  fric¬ 
tional  effects  into  account  directly  in  determining  the  in-situ  Ngo  from  Nb3o. 


Spectral  Analysis  of  Surface  Waves 


The  Spectral  Analysis  of  Surface  Waves  (SASW),  has  emerged  as  a  poten¬ 
tial  alternative  to  the  BPT  for  assessing  the  liquefaction  potential  of  gravels, 
cobbles  and  rockfill.  The  distribution  of  shear  wave  velocity  in  the  coarse 
materials  required  for  use  in  Figure  15  is  determined  by  analysis  of  surface 
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waves  only.  Therefore,  the  method  does  not  require  a  borehole.  The  method 
was  pioneered  by  Nazarian  and  Stokoe  (1984). 

Recent  developments  in  the  method  have  reduced  costs  and  accelerated  the 
procedure  so  that  it  has  become  a  viable  tool  for  routine  use  in  practice 
(Robertson  et  al.  1992).  Since  the  method  is  not  widely  used  yet,  a  brief 
explanation  and  some  typical  results  will  be  given  here. 

In  the  SASW  method,  Rayleigh  waves  are  generated  in  the  soil  by  steady 
state  harmonic  excitation  over  a  range  of  frequencies  or  by  hammer  impact  on 
the  ground  surface.  The  propagation  of  these  surface  waves  is  measured  by 
receivers  as  shown  in  Figure  20.  The  field  dispersion  curve  which  gives  the 
relationship  between  phase  velocity  and  wave  length  is  then  determined 
(Figure  21a). 

The  field  dispersion  curve  is  used  to  establish  the  shear  wave  velocity  pro¬ 
file  for  the  site.  The  process  of  going  from  dispersion  curve  to  velocity  profile 
is  an  inverse  process  and  is  basically  nonlinear.  An  iterative  forward  process 
is  commonly  used.  The  dispersion  curve  for  an  assumed  initial  velocity  profile 
(shear  wave  velocities  and  layer  thicknesses  over  a  half  space)  is  determined 
and  compared  with  the  field  curve.  The  initial  profile  is  modified  using  some 
form  of  optimization  until  the  theoretical  dispersion  curve  agrees  with  the  field 
curve  within  prescribed  limits.  Computed  and  measured  dispersion  curves  are 
shown  in  Figure  21a,  and  the  associated  shear  wave  velocity  is  shown  in 
Figure  21b. 

The  processing  of  the  receiver  data  previously  required  the  use  of  a  spectral 
analyzer.  The  spectral  decomposition  is  now  based  entirely  on  the  personal 
computer  using  improved  algorithms.  This  reduces  the  capital  cost  of  equip¬ 
ment  and  gives  a  faster  site  evaluation. 

The  SASW  method  is  fast.  Tokimatsu  et  al.  (1991)  quote  times  of  30  min 
to  1  hr  to  establish  a  profile  to  a  depth  of  50  m.  Robertson  et  al.  (1992)  quote 
times  of  1  hr  to  2  hr  for  a  profile  depth  of  30  m.  The  SASW  has  other 
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advantages.  It  does  not  require  a  borehole.  This  reduces  execution  time  and 
allows  the  method  to  be  used  at  gravel  and  rockfill  sites  at  which  is  difficult  to 
conduct  cross-hole  methods  or  downhole  methods  and  in  which  it  is  impossible 
to  conduct  seismic  cone  (SCPT)  measurements.  A  disadvantage  is  the  limited 
profile  depth  that  can  be  explored  with  common  convenient  energy  sources. 
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Figure  1.  Definition  of  liquefaction,  limited  liquefaction  and  cyclic  mobility 
under  cyclic  loading  (Vaid  and  Chern  1985) 


Figure  2.  Flow  chart  for  assessment  of  cyclic  loading  conditions 
(Robertson  1 990) 


Figure  3.  Relationship  between  stress  ratio  causing  liquefaction  and  (N,)6o 
values  for  silty  sand  for  M  =  7.5  (Seed  et  al.  1985) 


Figure  6.  Ranges  in  Kq  factors  (Seed  and  Harder  1990) 
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Figure  7.  Relationship  between  effective  vertical  stress  (o',,)  and  K-  (Seed  and 
Harder  1990) 


Figure  8.  -  Factors  for  clean  sand  as  a  function  of  relative  density  and 

sand  type  (Vaid  and  Thomas  1994) 
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Figure  9.  Relationship  between  Qc/NgpT  and  medium  grain  size  D50  numbers 
next  to  data  points  indicate  various  sources  of  data  (Robertson 
et  al.  1983) 
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Figure  10.  Relationship  between  critical  cyclic  stress  ratio  causing 

liquefaction  and  cone  tip  resistance  for  sands  and  silty  sands 
(Seed  and  De  Alba  1986) 


Cyclic  Strength  R  =  Oj/o^' 


Figure  1 1 .  Variation  in  cyclic  strength  with  fines  content  (Troncoso  1990) 


Figure  12.  Correction  for  fines  content  to  CPT  data  (Robertson  1990) 
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Figure  13.  Proposed  soil  classification  chart  based  on  normalized  CPT  and  CPTU 
data  (Robertson  1990) 


Figure  1 5.  Proposed  correlation  between  normalized  shear  wave  velocity  and  cyclic  stress 
ratio  (CSR)  to  cause  liquefaction  {Robertson  et  al.  1992) 
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Figure  17.  Correlation  between  corrected  Becker  and  SPT  blowcounts  (after 
Harder  and  Seed  1986) 


Figure  19.  (a)  Becker  blowcounts  N^so  for  mudded  hole  (BPT-66  FOB)  and  uncased  hole  (BPT-66 
RIO);  and  (b)  effect  of  friction  on  normalized  Becker  blowcounts 


Figure  20.  University  of  Alberta  SASW  field  equipment  (after  Addo  and 
Robertson  1992) 


Phase  Yeiccity(m/s)  Shear  Wave  Vdcaty(m/s) 
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Figure  21 .  Measured  and  computed  dispersion  curves  and  resulting  shear 
wave  velocity  (Tokimatsu  et  al.  1991) 


Appendix  D  -  Lecture  Slides 
Investigation  of  Liquefiable 
Soils  Using  Shear  Wave 
Velocity  -  Prof.  Ken  Stokoe 
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.  Discuss  Field  Seismic  Tests 
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FIELD  INVESTIGATIONS: 

•  SEISMIC  TESTING 

Crosshole,  Downhole,  SCPT  and 
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Simplified  Schematic  of  Crosshole  Seismic  Test 
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Seismic  Cone  Penetrometer  (SCPT)/Downhoie  Test 


SASW  Field  Testing  Configuration 
(SASW  =  Spectral  Analysis  of  Surface  Waves) 


Approximate  Distribution  of  Verticai  Particle  Motion  with 
Depth  for  Two  Rayleigh  Waves  with  Different  Wavelength: 
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Waves  Propagating  in: 

(a)  Uniform  Half-Space,  and 

(b)  Softer  Layer  over  a  Stiffer  Half-Space 


Wavelength, 


Rayleigh  Wave  Velocity,  Vr,  m/s 


Composite  Experimental  Dispersion 
Curve  from  Testing  with  Multiple 
Receiver  Spacings  at  Andersen  Bar 


Rayleigh  Wave  Velocity,  Vr,  m/s 


Match  Between  Experimental  and 
Theoretical  Dispersion  Curves 
at  Andersen  Bar 


Andersen  Bar 


150  300 

Shear  Wave  Velocity,  V  ,  m/sec 


Pence  Ranch 


Larter  Ranch 


riability  Study 


oon,  //  <  I  vV/v  \ 
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Variation  in  Normalized  Modulus  with  Strain  for  Imperial  Valley  Soils 
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Liquefaction  Potentiai  Charts  Based  on  and  T>  0.3 
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Liquefaction  Sites:  Gravelly  Soils 
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where  y  =  total  unit  weight,  and 

g  =  acceleration  of  gravity. 


o 

CM 


Velocity  Data  and  Values  from 
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SUMMARY 

4.  Liquefaction  Potentiai  Depends 
the  Site  and  Strong-Motion 
Characteristics: 
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Appendix  E  -  Lecture  Slides 
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GEOPHYSICS  IN  DAM 
SAFETY  INVESTIGATIONS 


I TOUEFACTION  INVESTIGATIONS 


*  PURPOSE  AND  OBJECTIVES 

A.  EMBANKMENT 

-  Old  Embankments 

-  New  Embankments 

B.  FOUNDATION 

-  Homogeneous  Geologic  Conditions 

-  Heterogeneous  Geologic  Conditions 

-  Unknown  Geologic  Conditions 


FIELD  METHODS  USED  FOR 
IDENTIFICATION  OF  LIQUEFACTION 
POTENTIAL 

A.  SURFACE  TECHNIQUES 

-  Refraction  Surveys 

-  Reflection  Surveys 

-  SASW  Surveys 

-  Electrical  Surveys 

-  Electromagnetic  Surveys 

B.  BOREHOLE  TECHNIQUES 

-  Crosshole  Seismic  Surveys 

-  Downhole  Seismic  Surveys 

-  Geophysical  Borehole  Logging 


GEOPHYSICS  IN  DAM 
SAFETY  INVESTIGATIONS 


LIQUEFACTION  INVESTIGATinN5;r 

(continued) 


C.  INTEGRATION  OF  TECHNIQUES 

-  Surface  and  Borehole  Surveys 

-  Geophysics  and  Geology 

-  Engineering  Applications 


*  CASE  HISTORY  AND  EXAMPLE 
INVESTIGATION 

A.  GRANBY  DAM  -  DIKE  No.  3.  COLORADO 

B.  McDonald  dam,  Montana 


*  PITFALLS  TO  WATCH  FOR  WHEN 
ASSESSING  LIQUEFACTION 
POTENTIAL  UTILIZING 
GEOPHYSICAL  METHODS 


GEOPHYSICS 


LIQUEFACTION  - > 

PURPOSE 

A.  EMBANKMENT  -  Determine  seismic 

STABILITY  OF  THE  STRUCTURE. 


•  Method  of  construction  and  equipment? 

Hydraulic  Fill/Sluiced  Core 
Homogeneous  Compacted  Fill 
Zoned  Earthfill 

•  Physical  Size? 

Structural  Height,  Length,  Crest  Width, 
Appurtenant  Features 


B.  FOUNDATION  -  Determine  seismic 

STABILITY  OF  CRITICAL  LAYERS 
WITHIN  FOUNDATION  SOIL  DEPOSITS. 

•  Homogeneous  Geologic  Settings  (Valley  Dams)? 

Flat  lying,  stratified,  and  continuous  deposits 
Typically  fine-grained  materials  (low  energy  systems) 

•  Heterogeneous  Geologic  Settings  (Mountain  Dams)  ? 

Unconformable,  unstratified,  and  discontinuous  deposits 
Typically  coarse-grained  materials  (high  energy  systems) 
Possibly  fine-  and  coarse-grained  materials  juxtaposed 
★Vertically  and  laterally  heterogeneous 


GEOPHYSICS. . . 


I .  LIQUEFACTION  - > 

*  OBJECTIVES: 

A-..  EMBANKMENT  -  Determine  latfpai  jmh 

y_ERTICAL _ SHEAR _ WAVE _ VELOCITY 

DISTRIBUTION! 


•Maximum  section  (through  cutoff) 

•  Downstream  area  of  critical  failure  - 

MIDSLOPE  OF  EMBANKMENT 

•  Upstream  area  of  critical  failure  -  midslope 

OF  EMBANKMENT 

•  Total  coverage  OR  selected  locations 

(27,000'  CREST  length) 


^ — foundation _ -  Determine  presence. 

MSTRIBUTION.  AND  CONTINUITY ^ 

OF  LIQUEFIABI  F  MATERIAL?;! 

•  Directly  beneath  embankment  crest 

•  Beneath  UPSTREAM  or  downstream  shells 

•  Along  downstream  or  upstream  toe 

•  Total  coverage  OR  selected  locations 

AGREED  UPON  BY  FET  AND  DESIGN  ENGINEERS  * 

-  "needs  vs.  wants"  ($)  - 


ELEVATION  (feel) 


GEOPHYSICS... 
I.  LIQUEFACTION  - > 


*  FIELD  METHODS  CURRENTLY  USED: 


A.  SURFACE  TECHNIQUES  -  Maximum  coverage  for 

MINIMUM  cost: 


1.  P-  AND  S-UAVE  REFRACTION:  Determine  lateral 

VELOCITY  VARIATIONS,  VELOCITY  LAYERING,  AND  DEPTH 
TO  SATURATED  MATERIALS.  OVERLAPPING  AND 
CONTIGUOUS  SPREADS  WITH  TIGHT  GEOPHONE  AND  SHOT 
SPACING  III. 
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FIGURE  1.  P-WAVE 
REFRACTION  PROFILE. 


FIGURE  2.  S-WAVE 
REFRACTION  PROFILE. 


T  iT«,.ir^.  geophysics... 

1.  LIQUEFACTION  - > 

*  FIELD  METHODS  CURRENTLY  USED: 

2.  P-WAVE  REFLECTION:  Determine  lateral  variability 

AND  CONTINUITY  OF  SUBSURFACE  LAYERS;  PARTICULARLY 
alluvium/bedrock  INTERFACE  OR  DENSE  LAYERS  UITHIN 
ALLUVIAL  DEPOSITS.  COMMON  OFFSET  GATHERS  (COG) 

ARE  WELL  SUITED  FOR  THIS  PURPOSE,  AS  WELL  AS 
DETECTION  OF  TRANSVERSE  CRACKS  IN  EMBANKMENTS  [2]  . 

•  Generally  no  velocity  information  or  cross-section 

properties;  therefore,  used  for 
^ALITATIVE  I NTERPRETATIONS . 


Shot 

No. 


Outance  (feet) 


COMMON  OFFSET  SURVEY 


GEOPHYSICS . . . 

I .  LIQUEFACTION  - > 

*  FIELD  METHODS  CURRENTLY  USED: 


3.  SPECTRAL  ANALYSIS  OF  SURFACE  WAVES  (SASW): 

Determine  Surface-wave  and  S-wave  velocity 

VARIATIONS  WITH  DEPTH.  REQUIRES  FREQUENCY 
SPECTRAL  ANALYZERS  AND  INVERSE  MODELING 

CAPABILITY  [31.  Surveys  are  expanded  about  a 

CENTER  POINT  YIELDING  APPARENT  (AVERAGE)  VELOCITY 
INFORMATION  BELOW  THAT  LOCATION. 


1 
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a.  Inversion 
Layering 


b.  Shear  Wave  Velocity  Profile 


GEOPHYSICS. . . 

I.  LIQUEFACTION  - > 

*  FIELD  METHODS  CURRENTLY  USED: 


4a.  electrical  SURVEYS:  Dipole-Dipole  Profiling  - 

GEOPHYSICAL  TECHNIQUE  REFERRED  TO  AS  "ELECTRICAL 

trenching"  (two  dimensional).  Used  to  identify 

ANOMALOUS  AREAS  OF  HIGH-  OR  LOW- RESISTIVITY  FOR 
l^YER  CONTINUITY  AND  MATERIAL  IDENTIFICATION. 

Changes  IN  resistivity  are  attributed  to  changes 

IN  AVERAGE  GRAIN  SIZE,  MOISTURE  CONTENT,  CLAY 
CONTENT,  POROSITY,  AND  PORE  WATER  CONDUCTIVITY. 

Shallow  depth  of  investigation,  approximately  40 

FEET,  FOR  SOD  FIELD  EXPLORATIONS  [43. 
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FIGURE  3,  DIPOLE-DIPOLE 
RESISTIVITY  PROFILE. 
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GEOPHYSICS... 

I.  LIQUEFACTION - > 

*  FIELD  METHODS  CURRENTLY  USED: 


4b.  electrical  SURVEYS:  Vertical  Electrical  Sounding  - 

VES:  GEOPHYSICAL  TECHNIQUE  REFERRED  TO  AS 

"electrical  drilling"  (one-dimensional).  Used  to 

IDENTIFY  CHANGES  IN  MATERIAL  PROPERTIES 
VERTICALLY.  EXPANSION  IS  ABOUT  A  CENTER  POINT  AND 
DEPTH  OF  INVESTIGATION  CAN  BE  100 'S  OF  FEET. 
Several  VES's  may  be  performed  to  develop  a  cross- 

section  FROM  THE  INDIVIDUAL  EARTH  MODELS  [5]. 

Battery  Current  meter 


Current  flow  through  earth 


VES  SOUNDING  STATIONING  9+00  27+50 


RMS  ERROR  %  2.53  3.79 


^  GEOPHYSICS... 

I .  LIQUEFACTION  - > 

*  FIELD  METHODS  CURRENTLY  USED: 

5.  ELECTROMAGNETIC  SURVEYS  -  EM  Profiling:  Very  quick 

AND  EFFICIENT  METHOD  FOR  DEFINING  SUBSURFACE 

material  changes;  primarily  laterally,  but  also 

VERTICALLY.  EQUIPMENT  MEASURES  CONDUCTIVITY  OF 

material  from  a  "focus"  depth,  dependent  upon 

TRANSMITTER  AND  RECEIVER  (cOIL)  SEPARATION  AND 
ORIENTATION  (VERTICAL  OR  HORIZONTAL  DIPOLES) 
NEED  TO  BE  CAUTIOUS  OF  BURIED  FERROUS  METAL 
OBJECTS  (i.E.,  conductors),  AS  WELL  AS  OVERHEAD  OR 
NEARBY  POWER  LINES  [6] . 
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Induced  current  flow  (hooogenous  halfspace) 
Vertical  dipole  (VD)  mode. 


Horizontal  dipole  (HD)  node. 
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ELECTROMAGNETIC  RESISTIVITY  DATA  PROFILES 


GEOPHYSICS... 
I .  LIQUEFACTION  - > 


*  FIELD  METHODS  CURRENTLY  USED: 

B.  BOREHOLE  TECHNIQUES  -  Drilling  required 

TO  PERFORM  SURVEYS: 


1.  CROSSHOLE  SEISMIC  SURVEYS:  Preferred  method  for 

DETERMINING  VERTICAL  DISTRIBUTION  OF  P-  AND  S-WaVE 
VELOCITIES.  An  ASTM  Standard  exists  (D4428  M-84) . 

The  diagnostic  geophysical  tool  used  by 
Reclamation  for  evaluating  liquefaction  potential 
OF  soils:  1989  Design  Standards  indicate  - 

•  Vg  >1200  ft/s  indicates  nonliquefiable  materials. 

•  Vg  BETWEEN  800-1200  ft/s  NONLIQUEFIABLE;  HOWEVER,  SUPPORTING 
DATA  MUST  BE  USED  IN  THE  DYNAMIC  ANALYSIS. 

•  Vg  <800  ft/s  indicates  potentially  liquefiable  materials. 
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,  GEOPHYSICS... 

I.  LIQUEFACTION  — — > 


*  FIELD  METHODS  CURRENTLY  USED: 


B.  BOREHOLE  TECHNIflllFS 

2.  DOWNHOLE  SEISMIC  TESTING:  Alternate,  less 

EXPENSIVE  METHOD  FOR  DETERMINING  VERTICAL 
DISTRIBUTION  OF  P-  AND  S-WaVE  VELOCITIES.  OnE 
DRILL  HOLE  REQUIRED,  TYPICALLY  SHALLOW  (<70  FT) 
SURVEYS  WHERE  VERTICAL  VARIATION  OF  MATERIALS 
ENCOUNTERED  IS  MINIMAL.  LeSS  RESOLUTION  AND 
"average  or  apparent"  VELOCITY  PROFILES  ARE 
OBTAINED  (compared  TO  CROSSHOLE  DATA) . 

APPARENT  S-UAVE  VELOCITY  (ft/s) 
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CROSS-SECTION 


Geologic  contact  depths  based  cm  drill  hole  log  of  DH-206. 

Seismic  interface  based  on  interpretation  of  S-wave 
velocity  cantrast. 
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GEOPHYSICS... 

I.  LIQUEFACTION  - > 

*  FIELD  METHODS  CURRENTLY  USED: 


B,  BOREHOLE  TECHNIQUES 

3.  BOREHOLE  GEOPHYSICAL  LOGGING:  Technique  involves 

LOGGING  BOREHOLES  WITH  DIFFERENT  GEOPHYSICAL 
PROBES  THAT  MEASURE  DIFFERENT  PHYSICAL  PROPERTIES 
OF  THE  SURROUNDING  MATERIALS.  EaCH  LOG  CAN  BE 
INTERPRETED  SEPARATELY,  OR  THE  "SUITE"  OF  LOGS 
FROM  AN  INDIVIDUAL  BORING  CAN  BE  INTERPRETED 
THROUGH  LOG  ANALYSIS  TO  DETERMINE:  POROSITY, 
MOISTURE  CONTENT,  RELATIVE  FRACTIONS  OF  CLAY/SILT- 

sand/oversize,  and  density.  There  are  specific 

REQUIREMENTS  FOR  DRILL  HOLE  COMPLETION  (E.G., 
STEEL  OR  PYC  CASING,  DIAMETER,  AND  GROUTING)  FOR 
EACH  BOREHOLE  LOG. 


T  .  GEOPHYSICS. . 

I .  LIQUEFACTION  - > 


*  to  watch  for  when 

assessing  liquefaction 

P2I.^!jTIAL  utilizing  GEOPHYSICAL 

METHODS: 

?.  interpreters  nightmare  -  HIDDEN  LAYERS 

OW  VELOCITY  LAYERS  BENEATH  HIGH  VELOCITY  LAYERS) : 
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’Tools  of  the  Trade" 

by  Bob  Dewey,  Bureau  of  Reclamation 
representative,  U SCOLD  Newsletter 
In  at  least  this  and  the  next  two  issues  of 
the  USCOLD  Newsletter,  the  1992 
issues,  you  will  find  a  series  of  articles 
featuring  some  very  useful,  unusual  or 
otherwise  special  techniques  available 
for  the  analysis  and  design  of  dams.  The 
purpose  of  this  series  is  to  provide  more 
exposure  for  tools  of  our  trade  that  have 
proven  to  be  highly  valuable  and  are  not 
well  known  or  sometimes  overlooked.  In 
this  issue  we  wU  focus  on  the  use  of 
geophysics  for  soil  characterization,  in 
particular  its  use  for  investigating  the 
liquefaction  susceptibility  of 
embankment  dams  and  their  foundation 
materials.  Future  articles  will  include 
the  use  of  petrographies  and  other 
state-of-the-art  tools. 

The  possibility  of  continuing  this  series 
beyond  1992  will  be  considered  if 
numerous  interesting  ideas  for  subject 
matter  are  su^ested.  Please  contact 
Bob  Dewey  at  the  Bureau  of 
Reclamation  (303-236-4196),  if  you 
would  like  to  write  about  "Tools  of  the 

Trade"  vdiich  would  be  of  engmeering 

interest  to  our  profession. 


Tools  of  the  Trade 

Geophysics  for  Soil 
Characterization  at 
New  and  Existing 
Dams 

by  Phil  S  tries.  Bureau  of  Reclamation 

Introduction 

With  the  profession’s  commitment  to 
dam  safety  there  has  been 
implementation  of  as  many  "tools  of  the 
trade"  as  necessary  to  characterize  in 
situ  soil  conditions.  Within  the  realm  of 
site  characterization,  determming  the 
liquefaction  potential  of  an 
embankment  dam  and/or  its  foundation 
is  just  one  example  of  the  engineering 
concerns  associated  vdth  Reclamation  s 
Safety  of  Dams  projects.  For  many  years 
the  primary  methods  utilized  for  soil 
characterization  were  drilling,  test 
pitting,  sampling  and  lab  testing. 
However,  over  the  past  decade 
engineering  geophysical  surveys  have 
become  a  major  component  of  field 
exploration  programs  via  the  integration 

of  the  geophysical  results  with 

site-specific  geolo^c  data  to  enhance 
characterization  of  the  soil  conditions. 
Using  a  "phased"  approach  to  soil 
characterization  permits  a  group  of 
interdisdplinary  personnel,  known  as  an 
exploration  team,  to  interactively 
evaluate  the  extent  and  most  efficient 
methods  for  the  field  exploration. 

This  article  presents  the  current 
philosophy  for  using  geophysical  "tools" 
so  that  characterization  of  soil 
conditions,  in  particular,  investigations 
for  liquefaction  assessment,  may  be 
completed  as  thorou^y  and 
economically  as  jwssible.  There  should 
be  no  pretense  that  geophysical  results 
are  to  be  used  as  definitive  or  final 
answers  during  field  explorations; 
geophysical  investigations  are 
performed  to  augment  and  facilitate  the 
exploration  team  to  conduct  their  soil 
characterization.  Information  presented 
in  this  article  demonstrates  that 
incorporating  accepted  geophysical 
field  methods  will  lead  to  better 


evaluation  of  site  soil  conditions  in  a 
cost  effective  manner. 

Typically  three  distinct  phases  are  used 
for  field  explorations  after  the  initial 
recommendations  for  conducting  field 
explorations  have  been  identified  in  a 
Safety  Evaluation  of  Existing  Dams 
(SEED)  report.  The  exploration  team 
convenes  to  discuss  available 
information,  site  conditions,  and  review 
the  recommendations  identified  in  the 
SEED  report.  During  this  meeting  the 
team  agrees  upon  methods  for 
evaluating  the  soil  conditions,  prepares 
a  field  program  for  Phase  1 
investigations,  and  puts  together  a 
strategy  for  later  phases  of  exploration. 
For  most  projects,  geophysics  will  be 
conducted  in  the  first  phase  of  field 
investigations,  but  for  some  SEED 
recommendations  or  at  some  damsites 
geophysics  may  not  be  warranted. 

Phase  1  Exploration 
Phase  1  explorations  result  in  a  general 
soil  characterization  through  literature 
review  of  existing  information,  geologic 
mapping,  and  surface  geophysical 
measurements.  The  purpose  of  Phase  1 
is  to  establish  if  areas  or  depth  zones 
critical  to  the  dam  stability  (Le., 
anomalous  conditions)  are  present  at 
the  site;  also,  whether  or  not  both  the 
embankment  and  foundation  require 
evaluation.  Objectives  of  the  Phase  1 

geophysical  program  leading  to 

identification  of  anomalous  conditions 
arc:  1)  determine  material  strengths  for 
the  embankment  and  foundation;  2) 
identify  the  types  of  materials  present, 

3)  evaluate  the  lateral  continuity  of  the 
soil  deposits;  4)  estimate  the  depth  to 
competent  soils  or  rOck;  and  5) 
determine  depth  to  the  water  table.  If 
further  field  explorations  are  required, 
results  from  the  Phase  1  geophysic^ 
surveys  assist  the  exploration  team’s 
dedsions  regarding  the  most  practical 
and  effective  methods  of  investigation 
and  where  they  should  be  located. 

To  address  Phase  1  objettives, 
geophysical  investigations  are  generally 
performed  along  the  dam  crest  and 
downstream  toe.  T^ically,  the 
investigations  consist  of  seismic 
refraction  and  electrical  resistivity 
surveys.  These  two  geophysical  methods 
are  briefly  described  below. 


Seismic  Refraction  Surveys 

Refraction  surveys  are  performed  to 
determine  the  velocities  of  subsurface 
layers  and  the  lateral  and  vertical 
distribution  or  variations  of  the 
velocities.  Refraction  surveys  are 
divided  into  two  techniques: 
compressional-  (P)  and  shear-  (S)  wave 
velocity  measurements. 

P-wave  velocities  obtained  in 
unconsolidated  soils  above  the  water 
table  provide  information  about  the 
material’s  in  situ  density  and  variability. 
Below  the  water  table,  the  P-wave 
velocity  is  affected  by  the  presence  of 
pore  water  and  saturation;  because  of 
this,  P-wave  refraction  surveying 
determines  depth  to  the  water  table. 
Additionally,  P-wave  surveys  are  used  to 
determine  the  depth  to  bedrock  or  a 
competent  soil  layer.  P-wave  velocities 
obtmned  in  competent  materials  cjin  be 
used  to  estimate  their  rippability  and 
erodibility  for  construction  purposes. 
Refraction  survey  results  are  presented 
in  cross-sectional  format  as  shown  on 
Figure  L  The  P-wave  seismic  interfaces 
between  unsaturated  materials, 
saturated  materials  (the  water  table), 
and  bedrock  are  displayed  on  the 
profile. 
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Figure  1.  P-wave  refraction  profile. 
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The  seismic  stiffness  of  a  material  is 
defined  by  its  shear  modulus,  which  can 
be  computed  after  measuring  the 
material’s  in  situ  density  and  S-wave 
velocity.  A  material’s  S-wave  velocity  is 


not  affected  by  the  presence  of  pore 
water,  hence,  S-wave  refraction  surveys 
can  determine  the  stiffness  of  the 
deposit  above  and  below  the  water 
table.  Dynamic  analyses  have  shown 
that  the  S-wave  velocities  at  soil  sites  are 
not  as  important  as  the  vertical  and 
lateral  distribution  of  the  velocity. 
Therefore,  the  purpose  of  performing 
S-wave  refraetion  surveys  is  detection  of 
the  lateral  and  vertical  variability  of 
velocity  for  the  embankment  and  its 
foundation.  Figure  2  displays  S-wave 
data  collected  at  the  same  location  as 


Rgure  2.  S-wave  refraction  profile. 


Figure  1.  Since  the  water  table  does  not 
affect  S-  waves,  the  embankment/ 
foundation  contact  was  the  only  seismic 
interface  detected.  The  embankment/ 
foundation  contact  is  slightly  different 
for  the  two  refraction  profiles  because 
"bedrock"  is  a  weathered  shale  and 
P-waves  could  not  differentiate  between 
residual  soil  and  saturated  embankment 
materials. 

In  general,  refraction  surveys  have  the 
limitation  of  sampling  large  volumes  of 
the  subsurface  which  produces  an 
average  velodty  for  each  layer.  The 
value  in  refraction  surveying  is  the 
assessment  of  soil  conditions  over  large 
areas.  Also,  determining  the  S-wave 
velocities  (high-  and  low-velocity  areas) 
permits  a  qualitative  judgment  about 
the  liquefaction  potential  for  the  soils. 

At  Reclamation,  a  conservative 
approach  is  used  for  liquefaction  of 
soils:  S-wave  velocities  less  than  800  ft/s 
indicate  materials  that  are  potentially 
liquefiable;  S-wave  velpcities  ranging 
between  800  and  1200  ft/s  are 
considered  potentially  non-liquefiable, 
but  supporting  evidence  is  required 


(e.g.,  gradation,  S  PT  blow  counts, 
percent  saturation);  and,  if  the  S-wave 
velocities  are  greater  than  1200  ft/s  the 
materials  are  considered  non- 
liquefiable. 

Electrical  Resistivity  Measurements 

The  major  factors  controlling  the 
electrical  resistivity  of  earth  materials 
are  porosity,  moisture  content  (degree 
of  saturation),  pore  water  conductivity, 
grain  size,  and  clay  content.  A  material’ 
resistivity  generally  increases  as  any  of 
these  parameters  decrease,  except  that 
resistivity  generally  decreases  with 
decreasing  average  grain  size.  Because 
of  surface  conduction  effects,  the 
presence  of  even  a  small  amount  of  clay 
in  a  soil  deposit  often  results  in  a  large 
decrease  in  the  measured  resistivities. 
Conversely,  because  the  resistivity  of  air 
is  essentially  infinite,  an  air-filled  void, 
or  a  region  having  a  high  percentage  of 
air-filled  pore  space  (e.g.,  porous  media 
above  the  water  table)  will  have 
extremely  high  resistivity.  Since  the 
electrical  resistivity  of  a  soil  or  rock 
depends  on  the  combined  effect  of  all  of 
the  parameters  listed  above,  resistivity 
values  generally  cannot  be  correlated 
with  spedfic  soil  or  rock  types  until  site 
geologic  data  are  correlated  with  the 
resistivity  data.  The  lateral  and  vertical 
changes  of  subsurface  resistivities 
obtained  at  a  site  provide  an  indication 
of  relative  changes  in  soil  or  rock  tyjjes. 

Similar  to  seismic  testing,  electrical 
surveys  are  also  divided  into  two 
methods:  dipole-dipole  resistivity 
profiling  and  vertical  electric  sounding 
(VES).  Dipole-dipole  profiling, 
sometimes  referred  to  as  "electrical 
trenching",  provides  resistivity 
information  "focused"  at  depth  along  an 
entire  survey  line.  For  engineering 
projects  the  technique  has  a  limited 
depth  of  investigation,  generally  about 
40  feet.  Dipole-dipole  measurements 
provide  a  continuous  apparent 
resistivity  profile  called  a 
pseudosection.  To  determine  the  true 
in-situ  resistivity  distribution 
corresponding  to  the  apparent 
resistivity  data  measured  in  the  field, 
two-dimensional  forward  modelling 
must  be  performed.  The  interpreted 
model  (true)  resistivities  and  area 
boundaries  are  displayed  superimposed 
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on  the  pseudosection  contours.  The 
modeled  true  resistivities  are  presented 
in  cross-sectional  format  (Figure  3). 
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Figure  3.  Dipole-Dipole  resistivity  profile. 

The  data  presented  in  Figure  3  were 
collected  over  laterally  heterogeneous 
glacial  deposits  where  coarse-  and 
fine-grained  materials  (lugh-  and  low- 
resistivities,  respectively)  are  juxtaposed 
due  to  the  depositional  environment. 
The  high  resistivity  areas  indicate 
outwash  channels,  whereas  the  low 
resistivities  are  indicative  of  point-bar 
deposits. 

VES  measurements,  sometimes  referred 
to  as  'electric  drilling*,  provide  a 
one-dimensional  earth  model  of  the 
resistivity  layering  beneath  a  single 
measurement  station.  Generally,  more 
vertical  resolution  of  the  subsurface 
layering  and  greater  depth  of 
investigation  are  obtained  from  VES 
data.  Computed  resistivity  models 
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•Figure  4.  Vertical  electric  soundings. 


represent  interpretation  of  the  field 
data  through  an  inversion  procedure. 

The  earth  models  are  typically  plotted 
in  "stick  log"  fashion,  since  the  resistivity 
data  represent  the  layering  beneath  an 
array  center  point.  Figure  4  shows  two 
VES  earth  models  obtained  on  an 
embankment  crest,  one  collected  at  the 
same  location  as  the  P-  and  S-wave 
seismic  data.  The  high  resistivity  surface 
layer  represents  unsaturated  shell 
(Zone  2)  materials,  while  the  low 
resistivity  middle  layer  represents 
fine-grained  core  (Zone  1)  materials, 
and  the  deepest  layer  of  moderate 
resistivity  represents  the  foundation 
bedrock  (shale). 

Incorporating  both  the  Phase  1  seismic 
and  electrical  results  permits  general 
characterization  of  material  strengths 
and  t3'pes,  respectively.  The  exploration 
team  then  meets  to  discuss  the  Phase  1 
results  and  decides  what  is  required  to 
verify  (ground  truth)  the  geophysical 
results,  and  also  discusses  the  optimal 
field  techniques  for  Phase  2 
explorations.  Drill  holes,  test  pits,  and 
material  sampling  can  be  located  and 
performed  in  the  anomalous  areas  as 
determined  by  the  Phase  1  geophysics. 

Phase  2  Exploration 

Phase  2  field  explorations  rarely  involve 
additional  geophysical  surveys;  rather, 
drilling^  sampling,  and  in  situ  and 
laboratory  testing  are  performed.  The 
types  of  field  investigations  used  during 
Phase  2  may  be  highly  varied  (e.g.. 
Standard,  Becker,  or  cone  penetration 
testing) .  The  phased  approach 
anticipates  using  correlations  between 
the  Phase  1  geophysical  and  Phase  2 
geological  data. 

Physically  measuring  and  testing  the 
Phase  2  material  samples  results  in 
properties  which  can  be  correlated  with 
the  geophysical  information.  This 
correlation  permits  "calibration"  of  the 
Phase  1  results  on  a  site-specific  basis. 
That  is,  by  combining  the  geophysical 
and  physical  sets  of  data,  the 
determination  of  lateral  continuity  of 
material  properties  between  drill  holes 
can  be  made. 

Material  properties  obtained  in  the  field 
during  Phase  2  must  be  considered  as 
investigations  are  progressing,  such  that 


the  exploration  team  has  flexibility  for 
revising  the  original  concept  of  the 
Phase  2  program.  If  the  correlation 
between  Phase  1  and  2  data  is  poor  then 
additional  exploration.  Phase  3,  maybe 
warranted.  An  example  encountered 
frequently  during  Phase  2  is:  the  Phase  1 
refraction  S-wave  measurements  show 
high  velocities  (  >  1200  fps)  for  the  near  • 
surface  materials  and  they  increase  with 
depth,  indicating  non-liquefiable  soil 
deposits;  but,  during  drilling  the  SPT 
blow  counts  are  moderately  high  to  high 
for  the  upper  soil  layers,  then  they 
decrease  with  depth  to  very  low  values 
indicating  the  presence  of  potentially 
liquefiable  layers  at  depth.  Tliis 
situation  does  not  imply  the  geophysics 
is  incorrect,  rather,  the  subsurface 
drilling  information  has  confirmed  that 
a  low-strength  layer  was  undetected  by 
refraction  surveying  and  Phase  3 
explorations  are  required. 

Phase  s  Exploration 

When  representative  geological  data  are 
obtained  during  Phase  2  and 
correlations  with  the  Phase  1 
geophysical  data  are  satisfactory,  Phase 
3  field  explorations  may  not  be 
required.  If  the  data  obtained  during 
Phases  1  and  2  are  not  sufficient  to 
characterize  the  soil  conditions'  at  the 
site,  supplemental  information  involving 
additional  geological  and  geophysical 
testing  must  be  acquired.  For  the 
example  previously  discussed,  where 
low-strength  soils  are  present  beneath 
high-strength  layers  (which 
geophysically  would  be  referred  to.  as  a 
"ghost"  or  "hidden"  laye^,  one  of  two 
geophysical  techniques  is  used  to 

measure  the  in  situ  S-wave  velodty 
distribution  with  depth:  CTOSshole 
testing  or  speciral  anal^is  of  surface 
waves  (SASW).  For  soil 
characterization,  each  geophysical 
method  has  its  advantages  and 
limitations.  Crosshole  testing  is 
considered  state-of-the-practice  and  is 
used  more  frequently;  however,  since 
crosshole  testing  requires  drill  holes 
there  are  c'lrcumstances  where  time 
constraints,  site  logistics,  or  the  project 
budget  does  not  permit  installation  of 
borings  and  SASW  testing  is  acceptable 
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40  percent  increase  in  S-wavc  velocity 
after  compaction,  with  an  effective 
depth  of  treatment  to  about  30  feet. 

Summary 

Implementing  a  phased  approach  to 
geotechnical  exploration  programs  for 
soil  characterization  permits  the 
exploration  team  to  acquire  the 
maximum  amount  of  site  information 
with  the  most  effective  methods. 
Typically,  three  phases  of  field 
exploration  wiU  sufficiently  characterize 
the  soils  by  evaluating  both  the  lateral 
and  vertical  distribution  of  material 
properties.  Assigning  an 
interdisciplinary  team  of  individuals  to 
supervise  all  phases  of  the  exploration 
program  enables  each  team  member  to 
bring  specialized  testmg  methods  to  the 
project  in  order  to  address  the 
engineering  concerns  of  the  site. 
Geophysical  surveys,  as  presented  in 
this  article,  should  be  used  during  the 
first  phase  of  exploration  to  obtain  an 
overall  site  assessment.  Phase  1 
geophysical  results  are  presented  to  the 
exploration  team  in  order  to  assist  in 
planning  the  Phase  2  investigations. 
Specified  field  testing  can  be  located 
and  performed  in  areas  or  depth  zones 
as  defined  during  Phase  1  which  are 
potentially  critical  to  the  dam  safety.  If 
spil  conditions  at  the  site  are  not 
sufficiently  characterized  after  Phase  2 
exploration,  or  the  correlation  between 
Phase  1  and  2  data  is  poor,  additional 
geological  and  geophysical  data  should 
be  acquired  in  Phase  3.  This  type  of 
approach  to  field  investigations  makes 
use  of  state-of-the-art  geophysical  "tools 
of  the  trade"  and  requires  cooperation 
among  geologists,  geophysicists,  and 
engineers  to  utilize  them  as  effectively 
as  possible.Et 


IDB  Loan  to  Nicaragua 

The  Inter-American  Development  Bank 
has  approved  a  $19  million  loan  to  help 
Nicaragua  rehabilitate  its  national 
electric  power  system  to  meet  current 
demand  and  allow  for  future  expansion. 
The  program  will  be  carried  out  by  the 
Nicaraguan  Power  Authority  and 
consist  of  four  components: 

1.  The  rehabilitation  of  the  spillways  of 
the  Asturias  and  Apanas  dams  and 
other  cml  works  of  the 
Asturias-Centroamerica-Carlos 
Fonseca  hydroelectric  complex. 

2.  The  rehabilitation  and  upgrading  of 
37.5  kilometers  of  transmission  lines  and 
substation  equipment  with  a  total  of  385 
megavoltamperes  of  transformer 
capacity. 

3.  The  procurement  of  meters,  service 
connections,  metering  equipment, 
vehicles  and  related  equipment  to 
regularize  service  to  illegal  users, 
reduce  losses,  build  feeder  lines  to  low- 
income  areas  of  Managua  and  improve 
INE’s  maintenance  capacity. 

4.  The  contracting  of  the  consulting 
services  required  to  strengthen  the 
environmental,  financial  and  planning 
procedures  of  INE. 

The  total  cost  of  the  program  is 
estimated  at  $36.5  miliion.  It  will  be 
co-financed  by  the  Central  American 
Bank  for  Economic  Integration  with  a 
$9.8  million  loan  that  includes  $8.9 
million  from  the  Mexican  government, 
and  by  the  Nordic  Development  Fund 
with  a  $3  million  loan.>2 
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Abstract 


This  paper  introduces  parametric  modeling,  commonly  known  as  system 
identification,  to  estimating  dynamic  soil  properties.  These  methods  are  opti¬ 
mized  for  short  and  non-stationary  signals.  The  Wildlife  Site,  CA,  subjected  to 
two  large  earthquakes  (Elmore  Ranch  and  Superstition  Hills)  on  November  24, 
1987,  was  used  in  an  example  analysis  since  the  associated  data  is  the  only 
publicly  available  record  of  buried  and  surface  motions.  During  the  second, 
larger,  earthquake  (Superstition  Hills  earthquake)  the  site  soils  liquefied.  The 
pseudo-stiffness  of  the  soil  system  between  the  two  accelerometers  is  defined  in 
order  to  present  the  changes  of  system  behavior  through  time  in  a  convenient 
form.  The  pseudo-stiffness  decreased  with  rising  pore  pressure  and  correlated 
with  the  evolution  of  the  pore  pressure  and  shaking  history.  Both  stationary  and 
recursive  parametric  modeling  methods  are  applied  to  the  unique  input-output 
data  set  from  the  two  earthquakes.  The  results  show  weakening  of  the  soils 
system  as  pore  pressure  increases,  with  estimates  of  natural  frequency,  damping 
ratio,  and  participation  factor  given. 


Introduction 


An  important  goal  for  the  geotechnical  engineer  is  the  ability  to  estimate  soil 
properties  without  intruding  into  the  soil  mass.  For  the  engineer  interested  in 
seismic  behavior  of  soils,  the  dynamic  properties  of  the  soil  are  of  interest,  par¬ 
ticularly  large-strain  properties.  The  archetypal  large  strain  field  excitation  is 
earthquake  strong  motion.  Ideally,  both  soil  motions  into  the  soil  layer  of  inter¬ 
est  and  soil  motion  on  the  surface  above  the  layer  are  recorded.  Given  this 
known  input  propagating  upward  from  depth,  and  the  output  at  the  top  of  the  soil 
column,  the  behavior  of  the  soil  can  be  modeled  by  inverse  theory.  If  a  suitable 
model  is  chosen  to  represent  the  system  of  interest,  the  estimated  model 
parameters  will  correspond  to  important  mechanical  parameters  of  the  system, 
such  as  damping,  natural  frequency,  and  stiffness.  This  estimation  of  parameters 
is  commonly  known  as  system  identification  (SI). 

The  traditional  method  of  geotechnical  analysis  of  dynamic  soil  motions  is 
through  the  Fourier  transform.  However,  serious  problems  arise  when  this 
method  is  applied  to  short  data  streams,  and  to  signals  changing  through  time  - 
non-stationary  signals.  This  study  was  undertaken  to  show  the  effectiveness  of  a 
different  type  of  model,  a  parametric  model  commonly  used  in  automatic  control 
and  geophysics,  which  avoids  many  of  the  limitations  inherent  with  calculating 
the  system  transfer  function  by  Fourier  techniques.  An  important  aspect  of 
particular  parametric  models  is  the  theoretical  link  between  the  estimated  system 
parameters  and  the  mechanical  parameters  of  a  lumped-mass  oscillator.  The 
parametric  model  allows  estimates  of  system  dynamic  properties  to  be  to  be 
made  if  an  input-output  data  set  is  available. 
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This  paper  applies  both  stationary  and  time-varying  parametric  models  to  the 
input-output  data  sets  recorded  during  the  two  (Elmore  Ranch  and  Superstition 
Hills)  1987  earthquakes  at  the  Wildlife  Site  in  the  Imperial  Valley,  California. 
Since  the  inverse  problem  (identifying  a  system  by  working  from  recorded  data) 
is  non-unique,  it  is  necessary  to  understand  the  mechanistic  behavior  of  the 
system  being  modeled  in  order  to  facilitate  selection  of  the  correct  model  and 
solution.  To  this  end  detailed  geotechnical,  geological,  and  seismological 
analyses  of  the  Wildlife  site  and  earthquake  events  are  presented.  To  gain 
clearer  understanding  of  the  relation  between  force,  displacement,  pore  water 
pressure,  and  material  stiffness,  a  new  method  of  combining  complex  time- 
dependant  variables  -  the  pseudo-stiffness  -  is  introduced. 

This  paper  will  start  by  introducing  and  explaining  parametric  modeling  in 
the  context  of  SI.  The  Wildlife  site  will  then  be  introduced  and  physical  prop¬ 
erties  examined.  The  geotechnical  behavior  of  the  site  during  strong  shaking 
will  be  examined,  in  particular  the  relation  between  strong  shaking,  pore  pres¬ 
sure,  and  soil  stiffness.  Finally,  the  SI  techniques  will  be  applied  to  the  Wildlife 
site  data  and  system  parameters  estimated. 
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System  Identification 


Parametric  Modeling 

The  goal  of  system  identification  is  to  model  a  system  in  a  manner  that  pro¬ 
vides  needed  mechanical  information  about  that  system.  The  most  common 
techniques  have  evolved  from  electrical  and  mechanical  engineering,  and 
involve  solving  the  inverse  problem  for  the  system  transfer  function.  Each 
method  has  limitations;  in  the  words  of  G.  E.  P.  Box,  “All  models  are  incorrect, 
but  some  are  more  useful  than  others.” 

The  process  of  inversion  allows  the  estimation  of  the  system  response  func¬ 
tion  (filter)  if  the  input  and  output  signals  are  known.  A  simple  model  for  char¬ 
acterizing  a  system  is  as  a  parametric  relationship  between  system  input  and 
output.  Such  a  model,  referred  to  as  an  autoregressive-moving  average  (ARMA) 
model,  is  based  on  discrete  time  series  analysis: 


+  ...  +  bx  +  bx  .  +  ...  = 


nb 


Vy=0  ’ 


*=1 


(1) 


where  y,  is  the  actual  output  data  sequence,  x,,  is  the  input  sequence  (assume 
white  noise  for  simple  spectral  estimation),  na  and  nb  are  the  AR  and  MA 
orders,  respectively,  and  the  subscript  is  the  time  step  counter.  The  output  is 
seen  as  a  combination  of  the  input  history  acted  upon  by  the  "b  ”  coefficients 
plus  the  past  outputs  acted  upon  by  the  “a”  coefficients.  The  input  series, 
involving  the  coefficients,  is  a  causal  moving  average  (MA)  process 
(convolutional).  The  series  involving  weighted  past  output  values  (“a” 
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coefficients)  is  a  noncausal  autoregressive  (AR)  process.  The  lengths  of  the  AR 
and  MA  processes  (model  order)  must  be  explicitly  chosen  so  that  the  model 
best  represents  the  process. 

Applying  the  shifting  theorem  to  Equation  1  yields  the  Fourier  transform 
(Bracewell  1978) 


+  a^e  + 


-)  = 


+  b^e  + 


...) 


(2) 


where  /  is  ^  and  m  is  circular  frequency.  Applying  the  Z-transform  (Brace- 

well  1978),  where  z'‘  =  and  rearranging,  yields  the  frequency  domain 
transfer  function  H„ 


H.. 


+  b^z  '  +  b^z^  +  ... 

1  -ajz'  -a^z^  +  ... 


(3) 


The  ARMA  model  is  very  powerful  in  that  it  can  easily  model  sharp  drops, 
sharp  peaks,  and  smooth  spectral  behavior.  It  is  also  the  most  parsimonious 
estimator  (Robinson  1982),  describing  a  complex  process  with  very  few  param¬ 
eters  calculated  from  a  small  length  of  data.  Parametric  modeling  avoids  many 
of  the  difficulties  inherent  in  the  traditional  Fourier  methods,  discussions  of 
which  can  be  found  in  many  books  and  journals  (e.g.,  Glaser  1993;  Johansson 
1993;  Pandit  1991).  Extensions  of  this  model,  e.g.,  ARMAX,  ARX,  Box- 
Jenkins,  allow  input,  system,  and  output  noise  to  be  expressly  modeled  (Ljung 
1987).  In  particular,  the  ARX  model  includes  the  effect  of  uncertainties  and 
noise  as  a  white  noise  term. 

The  ARMA  model  has  special  significance  since  it  can  be  derived  directly 
from  the  differential  equation  of  motion  for  an  A-degree-of-freedom  (DOF) 
system,  with  the  damping  ratio  and  resonant  frequency  as  the  model  parameters 
(e.g.,  Gersch  and  Luo  1970).  A  2n-2n  ARMA  model  is  therefore  a  valid  model 
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for  a  layered  soil  system,  or  soil-structure  interaction  problem.  The  damping 
ratio  and  resonant  frequency  of  the  /^-degree-of-freedom  oscillators  are 
contained  in  the  In  AR  parameters.  Phase  relations  are  preserved  in  the  MA 
parameters.  The  modal  frequencies  percent  of  critical  damping  (Wj,  (Ghanem 
et  al.  1991)  and  power  participation  factor  pj  (Pandit  1991;  Safak  1988)  are 
calculated  from  the  system  poles  and  residues  found  from  partial-fraction 
expansion  of  Equation  2.  The  modal  parameters  are  defined  as 


0).  = 
1 


(4) 


(5) 


p.  =  -r.conj(zp  -  z.conj(rp 


(6) 


where  X  =  Arg  (zj) ,  dj  =  -(0.  5)  In  |Zj|^  Zj  is  the  pole  for  modey,  is  the  residue 
for  mode  j,  and  At  is  the  digitization  rate. 

Traditional  methods  of  system  estimation,  both  parametric  and  non- 
parametric,  are  strictly  valid  only  for  stationary  data.  A  stationary  signal  is  one 
whose  statistics  do  not  change  with  time.  The  commonly  invoked,  loose 
definition  of  stationarity  requires  that  the  variance  of  the  signal  be  constant  over 
any  and  all  time  windows.  Inherent  in  a  stationary  transformation  is  the 
averaging  of  the  signal  components  over  the  sampling  duration.  The  difficulty 
with  non-stationary  signals  is  that  these  energies  are  changing  during  this  period. 
If  the  frequencies  present  are  changing  over  this  time  window,  the  resulting 
estimation,  regardless  of  method  used,  will  be  a  smeared  average  as  if  all  the 
frequencies  with  energy  were  active  without  change  throughout  the  period. 
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The  field  of  adaptive  filtering  was  developed  to  model  non-stationary  pro¬ 
cesses.  As  the  statistics  of  the  signal  change  through  time,  the  filter  “adapts”  to 
the  changing  variance  with  new  parameter  values  that  reflect  the  structure  of  the 
system  at  that  point.  The  predicted  value  for  the  next  time  step  can  be  compared 
with  the  actual  value,  and  the  difference  will  give  a  measure  of  how  well  the 
filter  is  doing  its  job  (Ljung  1987). 

A  popular  direct  adaptive  filter,  or  process  model,  is  the  Kalman  filter  (Kal¬ 
man  1960).  Sorenson  (1970)  points  out  that  the  Kalman  approach  is  a  direct 
descendant  of  Gauss’s  least  squares  mapping,  except  now  neither  the  signal  nor 
the  noise  model  is  stationary.  The  extended  Kalman  model  (Ljung,  1979)  con¬ 
stantly  updates  the  estimates  of  the  dynamic  process  by  examining  the  error 
between  actual  and  estimated  values  for  the  previous  time  step.  The  dynamics 
can  be  due  to  a  changing  input  or  noise  process,  or  it  can  be  due  to  the  system 
itself  changing.  The  effect  is  a  linearization  between  single  time  steps,  but  if  the 
system  is  changing  slowly  compared  to  the  time  step  used,  the  linearization  is 
invisible”  and  the  non-linear  behavior  is  well  modeled. 


Analysis  Methods  Applied  to  the  Soil  System 


For  the  problem  at  hand,  estimating  the  soil  parameters  during  the  Supersti¬ 
tion  Hills  earthquake  at  the  Wildlife  site,  both  the  system  input  and  output  are 
known.  Analysis  of  the  records  was  undertaken  using  the  standard  routines 
contained  in  the  MATLAB  System  Identification  Toolbox  (Math Works  1991). 
The  data  were  initially  analyzed  to  determine  if  they  could  be  modeled  as  sta¬ 
tionary  segments.  A  recursive  segmentation  scheme,  which  attempted  to  break 
the  data  into  segments  with  a  chosen  maximum  variance  (Ljung  1987),  did  not 
work  well  since  it  was  not  possible  to  determine  the  “correct”  variance  a  priori. 
Instead,  a  more  direct  method  was  used  the  output  simulated  by  the  calculated 
system  had  to  accurately  model  the  actual  measured  output. 
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The  input-output  data  record  was  broken  into  segments  based  on  a  mecha¬ 
nistic  understanding  of  the  seismic  event  (discussed  in  following  sections).  Ini¬ 
tially  is  was  assumed  that  the  various  segments  were  basically  stationary.  If  the 
stationary  model  could  not  accurately  and  parsimoniously  simulate  the  segment 
output,  a  non-stationary  recursive  model  was  used.  In  addition,  the  appro¬ 
priateness  of  the  model  was  checked  by  insuring  99  percent  confidence  in  both 
the  whiteness  of  the  residual  autocorrelation  function  and  the  cross-correlation 
function  between  the  input  and  output  residuals  (Bohlin  1987).  A  final 
verification  on  model  suitability  was  made  by  examining  the  ability  of  the 
parameters  to  model  independent  input-output  data  from  the  90°  Superstition 
Hills  and  the  complete  Elmore  Ranch  time  histories. 

The  stationary  algorithm  uses  a  least  squares  estimation  for  the  ARX  model. 
It  is  necessary  to  estimate  the  number  of  parameters  to  be  calculated,  which  is 
essentially  estimating  the  degree  of  freedom  of  the  soil  system.  There  is  no 
obvious  answer  to  the  degrees  of  freedom  of  the  system,  so  several  verification 
techniques  were  employed  to  insure  that  a  proper  model  order  was  estimated. 
Model  order  was  increased  in  2n-2n  steps  corresponding  to  an  additional  degree 
of  freedom  each.  The  simulated  output  of  the  model  was  then  compared  with 
the  actual  output  for  congruence,  and  the  fewest  parameters  needed  to  accurately 
characterize  the  system  was  chosen  as  the  system  model  order.  Examination  of 
the  pole  and  zero  plot  insured  that  excessive,  overlapping  parameters  were  not 
included  (Ljung  1987).  If  the  segments  proved  non-stationary,  they  were 
analyzed  using  a  recursive  Kalman  filter  technique  that  expressly  accounts  for 
non-stationarity  (Ljung  1987). 


System  Identification 
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introduction  and  Geography 

Two  large  earthquakes,  the  Elmore  Ranch  and  Superstition  Hills,  occurred  on 
November  24,  1987.  The  shaking  was  strong  enough  to  induce  liquefaction  at 
the  Wildlife  site,  which  suffered  more  than  2  m  of  lateral  spread  towards  the 
Alamo  river,  and  sand  boils  occurred  over  at  least  33  hectares  (Holzer  et  al. 
1989).  Accelerograms  from  buried  and  surface  transducers,  as  well  as  pore 
pressure  time  histories,  were  captured  for  both  the  Elmore  Ranch  and  Supersti¬ 
tion  Hills  earthquakes  (Brady  et  al.  1989).  This  is  the  only  full  input-output  data 
set  known  to  exist  for  a  site  experiencing  liquefaction,  and  the  only  set  of 
records  complete  enough  to  allow  an  unambiguous  system  identification 
analysis. 

The  U.S.  Geological  Survey  (USGS)  chose  to  install  an  instrumentation  array 
at  the  Wildlife  site  for  two  major  reasons:  it  liquefied  during  the  1 98 1 
Westmoreland  earthquake  (shear  magnitude,  M5,  6.0),  and  its  location  in  a  very 
seismically  active  area.  The  Wildlife  site  is  located  in  the  western  Imperial 
Valley,  13  km  north  of  Brawley,  Imperial  County,  California,  as  shown  in 
Figure  1 .  The  Imperial  Valley  is  one  of  the  most  seismically  active  regions  of 
the  United  States  (Hudnut  and  Sieh  1989).  The  Wildlife  site  is  located  in  the 
Imperial  Wildlife  Management  Area,  on  an  incised  flood  plain  of  the  Alamo 
river  (Bennett  et  al.  1984)  approximately  32  km  east  of  the  epicenter  of  the 
Superstition  Hills  earthquake.  The  deposits  consist  of  various  flood  plain,  flu¬ 
vial,  and  lacustrine  materials  with  7  distinct  units  in  the  first  26  m.  The  top  three 
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units  were  investigated  thoroughly  since  they  were  believed  to  be  the  seat  of  any 
liquefaction  (Holzer  et  al.  1989). 

Besides  surface  and  buried  (8  m)  accelerometers,  an  array  of  pore  pressure 
transducers  was  installed  at  various  depths  at  the  site,  making  it  one  of  four 
piezometric  arrays  in  the  U.S.  (Holzer  et  al.  1989;  Brady  et  al.  1989).  The  plan 
view  of  the  array  is  shown  at  the  top  of  Figure  2,  while  the  bottom  of  Figure  2 
presents  the  physical  relation  of  the  piezometers  to  the  soil  units.  Two 
piezometers  were  installed  in  subunit  Bl,  three  in  B2,  and  one  in  the  dense 
clayey  silt  of  Unit  D. 

Geotechnical  Description 


The  pertinent  geotechnical  properties  of  the  first  13  m  of  soil  at  the  Wildlife 
are  shown  in  Figure  3.  The  water  table,  controlled  by  the  near-by  Alamo  river, 
was  approximately  1.2  m  below  the  surface  when  the  Superstition  Hills 
earthquakes  struck.  Unit  A  is  a  very  soft  sandy  to  clayey  silt,  which  grades  into 
the  silty  clay  at  the  bottom  of  Unit  B  below.  The  upper  meter  of  Unit  B,  Bl ,  is  a 
very  loose,  poorly  graded,  silty  sand  which  phases  towards  the  denser  silty  sand 
of  subunit  B2.  In  general,  the  materials  in  Unit  B  grade  from  coarsest  at  the 
bottom  (7  m)  to  the  finest  at  the  top.  Brady  et  al.  (1989)  note  that  a  wood 
fragment  recovered  from  a  depth  of  6  m  was  estimated  to  be  about  230  years  old 
by  Carbon- 14  analysis.  Unit  C  ranges  from  a  medium-to-stiff  clayey  silt  to  a 
very  stiff  silty  clay.  Unit  D  is  a  dense,  well-graded  silt,  cemented  towards  the 
top  of  the  unit. 

The  uses  undertook  a  thorough  site  investigation  at  the  time  the  instrumen¬ 
tation  array  was  installed  (Bennett  et  al.  1984),  and  other  researchers  performed 
geophysical  investigations  to  determine  the  in  situ  shear-wave  velocity 
(Bierschwale  1984;  Nazarian  and  Stokoe  1984).  Bierschwale  reported  a  cross¬ 
hole  S-wave  velocity  of  about  100  m/s  for  subunit  Bl,  while  Nazarian  and 
Stokoe  report  a  spectral-analysis-of-surface-waves  S-wave  velocity  of  about 
85  m/s.  Both  standard  penetration  tests  and  cone  penetration  tests,  as  well  as  the 
soil  description,  indicated  that  subunit  Bl  would  be  the  most  likely  to  liquefy  in 
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the  event  of  strong  dynamic  excitation.  Piezometer  P5  and  P2  were  installed  in 
this  subunit. 


Geology 


The  Imperial  Valley  is  a  structural  depression  caused  by  the  active  spreading 
of  the  Gulf  of  California  (Magistral  et  al.  1989).  The  flanks  of  the  Valley  are 
the  original  continental  crystalline  plutonic  and  metamorphic  rocks,  about  2  km 
deep  on  the  west  side,  while  the  middle  of  the  central  Valley  is  comprised  of 
metamorphosed  sedimentary  rocks  buried  beneath  5  km  of  sediments 
(Magistrale  et  al.  1989). 

The  surface  sediments  in  the  areas  of  interest  are  Holocene  lacustrine  silty 
sand  and  claystone,  interbedded  with  alluvial  deposits,  below  which  are  slightly 
consolidated  Pleistocene  silty  and  clayey  lake  deposits  with  sand  and  gravelly 
units  (Sharp  et  al.  1989;  Bennett  et  al.  1984).  The  impedance  contrasts  (i.e., 
stiffness  ratio)  of  the  geology  are  such  that  the  vertical  motions  due  to  a  local 
earthquake  are  primary-waves  (P-waves),  while  the  horizontal  motions  are  man¬ 
ifestations  of  the  passage  of  shear- waves  (S-waves)  (Wald  et  al.  1990).  This 
assertion  is  borne  out  by  the  phase  shift  that  occurred  between  the  buried  and 
surface  horizontal  motion  records  and  by  the  lack  of  shift  in  the  vertical  records. 

The  Superstition  Hills  fault  is  a  northwest-trending  fault  that  can  be  seen  as  a 
continuation  of  the  San  Andreas  and  San  Jacinto  fault  zones  (Figure  1).  The 
region,  from  the  west  end  of  the  Imperial  Valley  to  the  Brawley  seismic  zone,  is 
marked  by  short,  northeast-trending  faults,  believed  to  correspond  to  step-backs 
in  the  crystalline  basement,  perpendicular  to  the  main  faults. 


Seismology 
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On  the  evening  of  November  23,  1987,  an  earthquake  occurred  along  the 
previously  unknown  Elmore  Ranch  fault,  with  a  Shear  Magnitude  (Mj)  of  6.2 
and  a  hypocentral  depth  of  about  1 1  km  (Magistrale  et  al.  1989).  The  fault  was  a 
north-east  trending,  left-lateral  displacement  made  up  of  a  multitude  of  small 
displacements  or  fault  zones.  The  fault  zone  appears  to  have  followed  along  a 
basement  step-back  and  continued  to  the  Brawley  seismic  zone.  The  maximum 
measured  surface  horizontal  displacement  was  125  mm,  which  occurred  during 
the  event  itself  There  has  been  no  measured  movement  along  the  zone  for  over 
2  years  after  initial  displacement  in  1987  (Sharp  et  al.  1989). 

Twelve  hours  after  the  Elmore  Ranch  earthquake,  a  larger,  long-duration 
earthquake  struck  on  the  Superstition  Hills  fault.  This  6.6  Mg  seismic  event  was 
caused  by  a  northwestward-striking  right  lateral  displacement  coincident  with 
previous  earthquakes  in  1968  and  1979  (Magistrale  1989).  The  two  1987 
earthquakes  Elmore  Ranch  and  Superstition  Hills,  are  associated  with  a  well- 
defined  conjugate  fault  system  (Hudnut  et  al.  1989).  The  left-lateral  motion 
relieved  normal  stress  on  the  existent  major  northwest-striking  fault  system, 
allowing  the  shear  stress  from  the  Pacific  plate  to  be  relieved.  Both  faults  appear 
to  be  geometrically  constrained  by  the  interface  between  the  crystalline 
basement  rocks  to  the  West  and  the  metamorphic  rocks  to  the  East  (Sharp  et  al. 
1989). 

The  Superstition  earthquake  was  a  complicated  event  characterized  by  an 
extended  period  of  strong  motion  (Wald  et  al.  1990).  Figure  4  shows  the  accel¬ 
eration  time  history  in  the  north-south  direction,  recorded  at  a  depth  of  8  m. 
There  appear  to  be  arrivals  of  energy  at  approximately  3  and  9  seconds 
following  arrival  of  the  initial  S-wave.  Thorough  and  independent  study  of  the 
strong  motion  and  teleseismic  records  by  several  researchers  indicates  the  over¬ 
all  earthquake  event  was  made  up  of  three  distinct  subevents  that  had  to  have 
occurred  over  an  extended  area  (Hwang  et  al.  1990;  Wald  et  al.  1990;  Frankel 
and  Wennerberg  1989).  The  physical  location  of  these  subevents  is  given  in 
Figure  5  in  both  the  map  and  plan  view.  While  all  three  studies  reached  essen¬ 
tially  the  same  general  conclusions,  the  study  by  Wald  et  al.  was  the  most 
exhaustive  and  their  quantitative  results  will  be  used  in  this  paper. 
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Wald  et  al.  (1990)  calculated  the  arrival  times  relative  to  initial  energy 
arrival,  seismic  magnitude,  and  horizontal  extent  of  the  three  subevents,  which 
are  presented  in  Table  1.  While  the  of  the  three  subevents  are  very  similar, 

the  energy  released  by  subevent  3,  characterized  by  the  seismic  moment,  is 
almost  3  times  larger  than  the  combined  energy  release  of  the  first  two  sub¬ 
events.  The  source  durations  of  the  first  two  subevents  were  short  (less  than  2  s 
for  subevent  2),  while  it  was  estimated  that  subevent  3  required  about  7.5  s  for 
the  displacement  to  propagate  the  18  km  (Wald  et  al.  1990).  The  calculated 
arrival  times  of  the  three  subevents  are  shown  in  Figure  4,  exactly  matching  the 
visually  perceived  arrival  of  packets  of  energy.  While  the  propagation  velocity 
of  each  subevent  was  greater  or  equal  to  the  S-wave  velocity,  the  very  slow 
rupture  velocity  for  the  entire  earthquake  implies  an  overall  rupture  mechanism 
reminiscent  of  progressive  failure  (Frankel  and  Wennerberg  1989). 


Pore  Pressure  History 


Analysis  of  pore  pressure  histories  recorded  during  the  Elmore  Ranch  and 
Superstition  Hills  earthquakes  yield  additional  understanding  of  the  mechanistic 
behavior  of  the  Wildlife  site.  Complete  pore  pressure  records  were  made  during 
both  temblors  at  five  different  depths  (Brady  et  al.  1989).  Records  show  that  an 
increase  in  pore  water  pressure  was  measured  during  both  the  Elmore  Ranch  and 
Superstition  Hills  earthquakes.  The  Superstition  Hills  records  have  been  of 
greater  interest,  since  for  this  excitation  the  pore  pressures  increased  enough  to 
cause  liquefaction.  This  was  proven  by  the  appearance  of  sand  boils  at  the 
Wildlife  site  after  the  Superstition  Hills  temblor  (Holzer  et  al.  1989). 

This  paper  will  focus  on  the  pore  water  pressure  history  of  piezometer  P5, 
which  was  buried  at  a  depth  of  2.9  m,  the  top  of  soil  layer  B 1 .  This  transducer  is 
being  singled  out  for  several  reasons;  it  most  clearly  indicated  that  the  excess 
pore  pressure  exceeded  the  total  overburden  pressure;  it  was  located  at  the  top  of 
the  layer  expected  by  the  experts  at  the  USGS  to  liquefy  during  the  Superstition 
Hills  shaking  and  therefore  (by  some  theories)  would  be  at  the  point  of  initial 
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liquefaction  (Scott  1987;  Florin  and  Ivanov  1962);  and  it  is  believed  that  P5  is 
the  transducer  in  the  array  most  likely  to  have  performed  correctly  (Hushmand 
etal.  1992,  1991). 


Figure  6  compares  the  surface  north-south  acceleration  (Figure  3;  Brady  et  al. 
1989)  with  the  pore  pressure  ratio  r„  (ratio  of  excess  pore  pressure  u  to  initial 
effective  overburden  stress  d^)  history  for  piezometer  P5,  Superstition  Hills 

earthquake.  The  pore  pressure  history  is  addressed  here  to  enhance  mechanistic 
understanding  of  the  recorded  strong  motion  signals  and  thus  simplify  and 
strengthen  later  interpretation.  There  is  some  controversy  whether  the  pore 
pressure  transducers  at  the  site  were  operating  properly  (Hushmand  et  al.  1992, 
1991),  however,  it  was  decided  to  utilize  the  must  uncontroversial  information 
available. 

The  initial  segment  of  the  P5  curve  (Figure  6)  is  a  sharp  “spike”  believed  to 
be  due  to  the  powering  up  of  electronics  upon  triggering.  This  spike  is  present  in 
the  records  of  all  5  piezometers  from  both  the  Elmore  Ranch  and  Superstition 
Hills  earthquakes.  At  time  t,,  (subscript  =  time  from  first  arrival  in  seconds)  the 
initial  S-wave  energy  from  the  first  subevent  arrived,  and  a  very  slight  increase 
in  pore  pressure  was  recorded  by  P5.  At  approximately  time  tj ,  energy  from 
subevent  2  excited  the  accelerometer,  and  the  pore  pressure  started  to  increase  at 
a  steady  rate.  Finally,  at  time  tj  g  strong  shaking  was  triggered  by  the  arrival  of 
subevent  3,  and  the  pore  pressure  ratio  at  P5  increased  at  a  rapid  rate  to  a 
maximum  greater  than  1.2.  After  tj  g  the  rate  of  increase  of  ru  was  relatively 
constant  for  the  duration  of  source  motion  (==7.5  s),  until  it  reached  0.6  (a  level 
often  considered  as  functional  liquefaction).  The  rate  of  pore  pressure  increase 
slowed  until  approximately  time  tjj,  at  which  time  r„  was  approaching  1  and  the 
rate  of  increase  reaches  equilibrium. 


The  Pseudo-Stiffness 


As  the  pore  pressure  ratio  increases,  the  effective  stress  acting  among  the 
sand  grains  decreases,  and  the  shear  modulus  (G)  decreases.  Since  it  was 
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impossible  to  accurately  measure  stress  and  strain  inside  the  soil  body,  the  ratio 
of  applied  force  to  displacement  (stiffness)  is  used  to  characterize  material 
behavior.  A  simple  model  for  the  soil,  system,  which  allows  calculation  of  a 
dynamic  pseudo-stiffness  for  the  liquefying  soil  layer  between  the  two  accele¬ 
rometers,  was  chosen  for  several  reasons.  The  soil  profile  of  the  Wildlife  site 
above  layer  B1  is  relatively  uniform  and  the  particle  motion  due  to  S-waves  is 
geologically  restrained  to  the  horizontal.  In  addition,  the  pore  pressure  history 
indicated  location  P5,  at  the  top  of  layer  Bl,  was  the  point  of  initial  liquefaction 
(expected  from  Scott  1986j  Florin  and  Ivanov  1961).  Once  the  pore  pressure  in 
layer  Bl  begins  to  rise,  the  sand  in  this  layer  is  very  soft  in  comparison  to  the 
rest  of  the  system.  Therefore,  the  motion  of  the  surface  accelerometer  relative  to 
the  buried  accelerometer  essentially  reflects  the  behavior  of  the  soil  above 
layer  Bl  (i.e.,  layer  A).  Prior  to  softening,  the  stiffness  will  be  dominated  by  the 
weakest  soil  between  the  two  geophones,  which  is  within  the  top  3.5  m. 

The  here-proposed  pseudo-stiffness  is  calculated  by  the  following  algorithm. 

a.  All  horizontal  motions  from  the  surface  (J  accelerometers  are  accounted 
for  by  taking  the  Euclidean  norm  of  the  north-south  and  east-west  surface 
acceleration  records:  . 


b.  All  horizontal  motions  from  the  buried  (d  accelerometers  are  accounted 
for  by  taking  the  Euclidean  norm  of  the  north-south  and  east-west  buried 


acceleration  records:  a. 


l^b.(N-S)  ^b,(E-W)- 


c.  The  relative  acceleration  of  the  surface  to  buried  transducer  is  determined 
from  the  step-by-step  difference  between  the  surface  and  buried  horizontal 
accelerations:  a  ,  =  a  ,  -a 


d.  The  relative  accelerations  are  converted  into  forces  per  unit  volume  by 
multiplying  the  differential  accelerations  acting  at  the  top  and  bottom  of 
the  layer  of  interest  by  the  soil  unit  mass,  which  should  change  relatively 


little  during  the  softening  process  (e.g.,  Udwadia  1985):  F.  =  a 


'■"'(0  ^  ■ 
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e.  The  running  total  of  the  amount  of  force  is  computed  for  the  system  during 
each  time  step:  sum(t)=sum(t-l)+[force(t)-force(t-l)J  -  yielding  the 
cumulative  horizontal  unit  force. 

f  Steps  1,  2,  3,  and  5  are  applied  to  the  displacement  records  to  yield  the 
cumulative  horizontal  displacement. 

g.  The  ratio  of  the  resulting  cumulative  horizontal  unitforce  to  displacement 
are  two  independent  variables  that  combine  to  give  a  pseudo-stiffness. 

The  use  of  a  single  transducer  to  provide  acceleration  and  displacement  mea¬ 
surements  is  not  problematic  since  the  derived  displacement  (Brady  et  al.  1989) 
is  generally  assumed  to  be  identical  to  what  would  be  recorded  by  an 
independent  displacement  transducer  placed  next  to  the  accelerometer.  The 
calculations  do  not  distinguish  between  active  (external,  inertial)  and  restorative 
(internal,  stiffness)  forces  acting  on  the  soil. 

The  pseudo-stiffness  behavior  due  to  the  Superstition  Hills  earthquake  is 
shown  in  Figure  7.  Each  point  on  the  curve  is  associated  with  a  particular  time 
step  of  the  earthquake  history.  Some  important  times  associated  with  changes  in 
the  earthquake  source  mechanism,  motion  record,  and  pore  pressure  history  are 
highlighted.  The  pseudo-stiffness  due  to  the  Elmore  Ranch  earthquake  was  also 
calculated,  and  is  shown  in  comparison  in  Figure  7.  The  pore  pressure  rise 
during  the  Elmore  Ranch  shaking  was  not  large  enough  to  force  the  sand  out  of 
elastic  behavior,  and  the  system  exhibited  the  same  stiffness  up  to  that  pore  pres¬ 
sure  (3.1  kPa,  associated  with  an  identical  pseudo-force)  12  hr  later  when  it  was 
excited  by  the  Superstition  Hills  temblor. 

The  interrelation  between  Superstition  Hills  strong  motion,  pseudo-stiffness, 
and  pore  pressure  ratio  is  shown  in  Figure  8.  Layer  B1  had  an  initial  pseudo¬ 
stiffness  of  3.2  MPA/m,  disturbed  slightly  by  increasing  pore  pressure  due  to 
arrival  of  energy  from  subevent  2.  With  the  arrival  of  subevent  3  at  time  t9  8  the 
pore  pressure  rapidly  increases,  and  the  pseudo-stiffness  shows  the  soil  dra¬ 
matically  weakening.  The  pseudo-stiffness  decreased  for  approximately  the 
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duration  of  the  subevent  3  motion,  becoming  a  constant  0.07  MPA/m  at  approx¬ 
imately  time  t25,  coincident  with  r„  =  1 . 

The  pseudo-stiffness  is  a  convenient  representation  of  the  system  behavior, 
combining  load,  time,  and  displacement  information  in  a  simple  curve.  The 
pseudo-stiffness  shows  that  the  soil  indeed  did  soften  in  relation  to  pore  pressure 
buildup.  In  addition,  the  increases  in  pore  pressure  and  decreases  in  stiffness 
occurred  when  additional  energies  arrived  from  subevents  2  and  3.  The 
understanding  of  the  physical  behavior  of  the  Wildlife  site  gained  in  these  sec¬ 
tions  will  allow  recognition  of  reasonable  estimates  in  the  following  analysis. 
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System  Identification  of  the 
Wildlife  Site 


The  input-output  data  record  was  initially  broken  into  segments  based  on  the 
mechanistic  understanding  discussed  in  previous  sections.  Segment  1  (see  Fig¬ 
ure  6)  is  the  segment  from  the  start  of  the  time  history  to  the  arrival  of  sub¬ 
event  2  at  approximately  t2.8  sec  (171  data  points).  Segment  2  runs  from  the 
arrival  of  subevent  2  to  the  arrival  of  subevent  3  at  approximately  tg  g  sec 
(150  data  points).  The  segment  corresponding  to  the  shaking  due  to  subevent  3 
is  assumed  to  be  very  non-stationary  due  to  the  nature  of  the  signal  and  the  rapid 
rise  of  pore  pressure,  and  was  not  individually  analyzed.  Segment  3 
encompassed  the  coda  of  the  signal,  after  the  shaking  stopped  at  approximately 
t,7  seconds,  to  the  end  of  the  record  at  93  sec  (1,870  data  points). 

The  first  two  segments,  corresponding  to  the  beginning  of  the  record  up  to 
the  arrival  of  subevent  3,  can  be  accurately  modeled  by  a  stationary  ARX  model. 
Segment  1  is  represented  as  a  3-Degree-of-Freedom  (DOF)  system,  with  a  root 
mean  square  (RMS)  error  between  simulated  and  actual  velocity  of  0.26. 
Segment  2  proved  to  be  representable  as  a  3-DOF  system,  with  an  RMS  error 
between  simulated  and  actual  velocity  of  1 . 1 .  A  comparison  of  the  spectral 
estimates  of  the  first  two  segments  is  shown  in  Figure  9,  along  with  the  estimates 
of  modal  parameters.  Table  2  presents  a  summary  of  the  ARX-calculated  modal 
parameters,  where  it  is  seen  that  the  natural  frequency  diminished  slightly  after 
the  arrival  of  subevent  2,  an  indication  of  system  softening.  This  is  to  be 
expected  since  the  pore  pressure  begins  to  increase  during  this  time  segment. 

The  change  in  the  soil  system  was  so  slight  during  the  initial  8.8  s  of  excita¬ 
tion,  that  it  proved  possible  to  model  the  combined  segment  accurately  with  one 
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4-DOF  ARX  model  (RMS  error  =  0.59),  with  the  characteristics  of  the  second 
subevent  overwhelming  the  lower  energy  first  arrival.  The  suitability  of  this 
model  is  verified  by  comparing  the  actual  surface  velocities  from  the  Elmore 
Ranch  360°  horizontal  record  with  the  simulated  output  for  the  Elmore  Ranch 
data.  The  pseudo-stiffness  and  pore  pressure  rise  from  Elmore  Ranch  was  iden¬ 
tical  to  that  of  the  Superstition  Hills  segment,  indicating  that  the  physical 
behavior  of  the  soil  was  the  same  for  this  degree  of  pore  pressure  increase.  The 
congruity  between  simulation  and  real-life  is  shown  in  Figure  10.  The  signals 
are  virtually  identical,  with  an  estimation  RMS  error  of  only  0.72.  The 
identified  model  performed  equally  well  on  all  the  other  records  available  for 
this  site. 

In  order  to  illustrate  the  time-varying  behavior  of  the  soil  (due  to  the  time- 
varying  pore  pressure)  the  initial  20  s  of  the  Superstition  Hills  record  was  ana¬ 
lyzed  with  a  recursive  technique  employing  a  Kalman  filter  model  with  3-DOF. 
The  calculated  model  was  able  to  simulate  the  actual  surface  velocity  history 
with  an  RMS  error  of  1 .27.  The  spectral  time  history  of  the  top  8  m  of  the 
Wildlife  site  for  the  initial  20  s  of  shaking  is  shown  in  Figure  1 1 .  In  this  figure 
the  smoothed  spectral  time  history  is  shown.  The  softening,  indicated  by  the 
downward  shift  in  natural  frequencies  and  increase  in  damping  corresponds  with 
the  initial  increase  in  pore  water  pressure  to  about  3.1  kPa.  This  behavior  also 
corresponds  with  decreasing  pseudo-stiffness  as  illustrated  in  Figure  8.  The 
smoothed  recursive  estimate  of  the  relevant  soil  parameters  is  given  in  the  top 
portion  of  Table  3.  On  the  whole,  the  recursive  estimates  for  important  times  in 
the  seismic  history  compare  well  with  the  stationary  ARX  estimates  from 
Table  2. 

The  system  behavior  during  the  subevent  3  strong  shaking  is  obviously  non¬ 
stationary  and  nonlinear  due  to  the  rapid  build-up  of  pore  water  pressure  and 
irregular  shaking.  The  point  at  which  a  stationary  model  could  no  longer  model 
the  system  behavior  corresponds  to  the  time  when  the  pore  pressure  ratio 
reached  approximately  0.6,  at  t,65.  The  spectral  estimates  from  a  Kalman  filter 
model  showed  that  the  natural  frequency  diminished  and  damping  ratio 
increased  as  the  effective  stress  in  soil  layer  B  decreased. 


System  Identification  of  the  Wildlife  Site 


Initial  visual  inspection  of  the  final  segment,  corresponding  to  the  inertial  soil 
motion  after  the  end  of  source  energy,  indicated  that  the  input-output  signal  had 
a  regular  oscillatory  pattern  and  would  therefore  be  easy  to  model  as  a  low-order 
ARX  system.  However,  this  segment  proved  to  be  quite  non-stationary  due  to 
the  constant  increase  in  phase  shift  between  buried  and  surface  signals,  and  the 
soil  could  not  be  accurately  modeled  as  a  time  invariant  system.  This  section 
lent  itself  well  to  rigorous  statistical  segmentation  by  defining  segments  such 
that  the  squared  error  between  model-predicted  and  actual  values  were  less  than 
0.  1.  Application  of  a  3-DOF  Kalman  filter  captures  the  essence  of  the  soil 
response  (RMS  error  =  0.76),  with  the  spectral  time  history  shown  in  Figure  12. 
As  indicated  by  the  time  dependent  values  for  coj,  pj  summarized  at  the  bottom 
of  Table  3,  the  soil  system  is  slowly  changing  (weakening)  as  the  pore  water 
pressure  comes  to  equilibrium.  The  energy  transmitted  through  the  system  is 
also  much  smaller  once  the  pore  ratio  approached  1 .  The  increase  in  damping  as 
the  sand  changes  to  a  viscous  fluid  is  responsible  for  dissipating  a  large  part  of 
the  energy  at  this  point. 


Discussion 


The  parametric  models  used  in  this  study  are  powerful  tools  for  the  qualita¬ 
tive  and  quantitative  estimation  of  in  situ  soil  properties.  The  recursive  methods 
allow  the  behavior  of  the  soil  system  to  be  monitored  throughout  the  excitation 
process  and  liquefaction,  as  shown  in  Figures  1 1  and  12.  These  two  figures 
(Figure  1 1,  two  initial  segments;  Figure  12,  coda)  show  the  changing 
characteristics  of  the  soil  system  through  time  by  plotting  the  spectral  density 
evolution  through  time.  This  unique  insight  into  the  liquefaction  process  eluci¬ 
dates  the  reaction  of  the  soil  system  to  incoming  energies,  and  eventual  softening 
due  to  the  rise  in  pore  water  pressure.  Confidence  in  the  estimated  models  is 
increased  by  comparison  with  independent  data  from  the  E-W  direction,  and  the 
Elmore  Ranch  temblor.  The  estimates  all  converged  on  the  values  reported. 

A  question  must  be  raised  about  the  meaning  of  the  damping  values  calcu¬ 
lated  in  this  study.  A  major  problem  is  that  the  damping  is  being  modeled  as 
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VISCOUS  damping,  which  is  known  to  be  a  poor  model  but  the  most  amenable  to 
calculation  and  the  one  most  commonly  used  (Pandit  1991).  Even  with  a  perfect 
technique,  the  damping  values  estimated  by  the  model  will  not  truly  represent 
what  is  physically  occurring  in  the  soil.  The  values  calculated,  reasonable 
compared  to  other  modal  studies  (e.g.,  lemura  et  al.  1990)  and  many  laboratory 
studies  (Vucetic  and  Dobiy  1991),  are  better  seen  as  “effective”  damping 
coefficients  which  capture  all  forms  of  frequency-dependent  mechanical  losses, 
rather  than  an  ‘intrinsic’  material  property. 


The  damping  estimates  are  often  very  sensitive  to  subtle  changes  in  the  mod¬ 
eling  of  the  system,  while  the  natural  frequency  estimates  are  more  robust 
(Gersch  1974).  This  is  especially  true  for  recursive  estimates  which  have  very 
different  damping  values  changing  rapidly.  In  this  case  a  limited  number  of  data 
points  enter  directly  into  the  calculation  and  variance  is  inversely  proportional  to 
the  square  root  of  the  number  of  data  points. 

The  physical  interpretation  of  the  instantaneous  mechanical  values  is  also  not 
immediately  clear.  The  ARMA  parameters  recursively  calculated  at  any  given 
time  define  the  filter  needed  to  transform  that  time  step  of  input  data  into  the 
next  output  value.  The  mechanical  characteristics  are  mathematically  extracted 
from  the  filter  based  on  a  limiting  set  of  assumptions,  and  an  “instantaneous 
frequency”  or  damping  has  no  physical  meaning.  The  ARX  estimates  are  made 
with  more  input  data  points  than  the  Kalman  estimates,  but  the  Kalman 
assumptions  better  model  the  changing  system.  There  is  no  obvious  choice 
between  the  ARX  or  Kalman  model  parameter  estimates  for  the  first  two 
segments  (-4  to  8.8  s).  A  fundamental  problem  remains — at  this  point  the  geo¬ 
technical  engineer  does  not  really  know  what  the  “correct”  answer  should  be! 
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Conclusions 


The  Wildlife  Site,  CA  was  subjected  to  two  large  earthquakes  on  Novem¬ 
ber  24,  1987.  During  the  second,  larger,  earthquake  (Superstition  Hills)  the  site 
liquefied.  The  data  collected  are  the  only  available  records  of  buried  and  surface 
motions  publicly  available.  In  addition,  pore  water  pressure  records  were 
recorded  at  several  depths  during  the  temblor. 

The  Elmore  Ranch  and  Superstition  Hills  Earthquakes  were  analyzed  from  a 
seismological  viewpoint  to  gain  understanding  of  the  site  behavior.  It  is  shown 
that  the  Superstition  Hills  earthquake  was  actually  three  distinct  subevents.  The 
strong  motion  and  pore  water  pressure  records  were  examined,  and  showed 
excellent  correlation  with  the  subevents.  A  significant  rise  in  pore  pressure  that 
resulted  in  liquefaction  of  layer  B  occurred  at  the  onset  of  subevent  3. 

This  paper  introduces  the  concept  of  pseudo-stiffness,  integrating  the  change 
of  system  stiffness  in  time  into  a  convenient  form.  The  pseudo-stiffness 
decreased  with  rising  pore  pressure,  and  correlated  very  well  with  the  evolution 
of  the  pore  pressure  and  shaking  history.  At  present  the  pseudo-stiffness  is  a 
qualitative  tool,  but  work  with  other  records  should  allow  its  quantitative  use. 

Both  stationary  and  recursive  parametric  modeling  methods  are  applied  to  the 
unique  input-output  data  set  from  the  two  temblors.  These  unique  results  show 
weakening  of  the  soils  system  as  pore  pressure  increased,  with  the  recursive 
estimates  being  summarized  in  Table  3.  While  there  are  some  important 
questions  to  be  raised  about  the  damping  estimates,  reasonable  estimates  were 
presented.  Further  work  using  system  identification  to  estimate  in  situ  soil 
properties  is  warranted. 
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Figure  1 .  Location  of  the  Wildlife  Site,  Imperial  Valley,  California,  with  the  locations -of  the  Elmore 
Ranch  (M^  6.2)  and  Superstition  Hills  (M,  6.6)  earthquakes  indicated  by  *  (Sharp  et  al. 
1989) 


Figure  2.  Summary  of  the  geotechnical  properties  of  the  Wildlife  Site  (Bennet  et  al.  1 984) 
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Figure  3.  Layout  of  the  instrumentation  at  the  Wildlife  Site  (Bennet  et  al.  1984) 
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Figure  4.  Acceleration  record  at  a  depth  of  8  m,  N-S  direction,  Superstition  Hills 
earthquake 
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Figure  6.  Pore  pressure  ratio  history  compared  to  top  N-S  acceleration.  Note  warm-up  detail  in  inset 
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Figure  7.  Relative  force-displacement  of  the  soil  system  during  the  Elmore. Ranch  and  Superstition 
Hills  earthquake 
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Figure  9.  Spectral  estimates  for  the  first  and  second  segments,  with  estimates  of  natural 
frequency,  damping  ratio,  and  participation  factor 
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Figure  10.  Comparison  between  the  Elmore  Ranch  velocity  history  simulated  by  the  second 
segment  Superstition  Hills  filter,  and  the  actual  velocity  history 


Figure  1 1 .  Soil  system  evolutionary  time  history  for  the  initial  20  s  of  the  Superstition  Hills 
earthquake.  The  changing  system  is  shown  in  spectral  representation 


Relative  Amplitude 


Table  1 

Parameters 

1990) 


of  the  3  Subevents  of  the  Superstition  Hills  earthquake  (data:  Wald  et  al. 


*  Shear  Magnitude  M,  |  Seismic  Moment  (x10“) 


Deiay  (s)  Source  Extent  (km) 


Table  2 


Summary  of  System  Parameters  Estimated  for  the  First  Two  Subevents  of  the 
Superstition  Hills  Earthquake  Using  the  ARX  (Stationary)  Model 


2.8:8. 8  I  2.4 


1st  Mode 


(%)  crh. 


2nd  Mode 


(%)  crrt. 


-0.30  10.7 


3rd  Mode 


(%)  crit. 


Tables 

Summary  of  Recursively-Estimated  System  Parameters  Evolving  Through  Time  During 
the  Superstition  Hills  Earthquake _  - 


l5t  Mode 


(%)  crit. 


2nd  Mode 


(%)  crH. 


3rd  Mode 


(%)  crit. 


72 


0.21 


45 


0.08 
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A  Nondestructive  Method  for 
Estimating  Properties  of  Soil  Subjected 
TO  Earthquake  Strong  Motion 


Waterways  Experiment  Station,  Earthquake 
Engineering  Site  Characterization  Workshop, 

June  27, 1995 


Steven  Glaser,  Division  of  Engineering,  Colorado  School  of  Mines 


PROBLEM:  Can  we  use  earthquake  strong  motion 
records  to  gain  understanding  of  the  non-linear 
behavior  of  soils? 


•  In  particular,  look  at  liquefiable  sands 

•  Vitally  important  knowledge  for  lifeline  safety 


SOLUTION:  System  identification,  In  Its  broadest 
sense,  can  yield  important  soil  parameters. 

•  Need  input-output  data  set 

•  Need  firm  mechanistic  understanding  of  the 

physical  events  that  occur 


study  the  Wildlife  Site 
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The  Wildlife  Site  In  the  Imperial  Valley,  CA 

•  Area  is  extremely  seismically  active 

•  Site  liquefied  during  the  1981  M  5.5 

Westmoreland  earthquake 
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SUPERSTITION  HILLS  EVENTS 


Relation  of  Earthquakes  to  Regional 

Geology 

Due  to  rapid  increase  of  velocity  with 
depth, 

-  P-waves  act  verticaliy 

-  S-  and  R-waves  act  horizontally 
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Instrumentation  and  Layout  of  The  Wildlife 

Site 

•  The  only  site  in  the  U.S.  having 

experienced  strong  motion  and  equipped 
with  buried  and  surface  accelerometers, 
and  piezometers 
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CModlfii 


Inut  £« 


November  24, 1 987  -  2  earthquakes  strike 

0154  GMT,  M  =  6.2  "Elmore  Ranch" 

1315  GMT,  M  =  6.6  "Superstition  Hills" 
The  temblors  occurred  on  conjugate  faults 


Elmore  Ranch  Temblor 


o 

lO 


Acceleration  (cm/sec^) 

o 

o  o  o 

lO  O  lO 


o 

o 


o 

lO 


{'^£>/n)  OjjBy  ejnssejd  ©JOd 


Time  (seconds) 

The  Elmore  Ranch  Temblor  is  a  Single  Event 

•  Some  Increase  in  pore  water  pressure  at  the  Wildlife  Site 

•  No  attendant  liquefaction  at  the  Wildlife  Site 


Sliikt  Slip  Dtilocfttion  (ctn) 


Event  2,  after  2.1  sec.,  M  =  6.0 

Events,  after 8.6 sec.,  M  =  6.4  (slip  «  18  km) 


Event  2 


Obvious  arrivals  of  the  subevents 

7.5  second  duration  of  event  3  apparent  from  record 


Pore  Pressure  Ratio  History  at  Piezometer  #5 

•  Note  3  distinct  segments 


(D/n)  ojiey  0jnss0Jd  9J0d 


•  Three  subevents  are  obvious  in  the  seismic  and  pore  pressure  records 

•  Change  in  pore  pressure  relates  directly  to  soil  motion 


Formation  Stiffness  is  a  Very  important 

Parameter 


Is  there  a  Simple  Way  to  Calculate  Stiffness 

Through  Time? 


Note  That  - 

•  Only  horizontal  motion  is  needed 

•  Acceleration  is  proportional  to  force 

•  Once  pore  water  pressure  starts  to  build, 

the  system  can  be  modeled  in  a 
simplified  manner 


Postulate  a  Pseudo-Stiffness  Curve 


1)  Take  Vector  sum  of  the  orthogonal 

horizontal  records,  for  both  surface  and 
buried  transducers. 

2)  Take  difference  between  surface  and 

buried  records. 

3)  Take  running  sum  of  acceleration- 

sum(i)=sum(i-1)+(accel(i)-accel(i-1)) 

4)  Multiply  curnulative  acceleration  by 

representative  unit  mass  to  yield  - 

Cumulative  Horizontai  Unit  Force. 

5)  Find  the  Cumulative  Horizontal 

Displacement  through  a  similar 
calculation. 


NOTE:  This  calculation  does  not  distinguish  between 
active  and  restorative  force. 


seudo-Force  Displacement  Curve  for  Superstition  Hills 

•  For  relative  motions 

•  Time  relation  shown  on  curve 


Pore  Pressure  Ratio  (u/a)  Cumulative  Horizontal  Unit  Force  (MPa) 


Top  Accolorstion  Compdrod  to  Psoudo-Stiffnoss 

•  At  the  rate  of  pore  pressure  buildup  starts  to  decrease,  r^  >0.6 

-  Material  starts  to  seriously  weaken 

-  End  of  source  events  (event  3) 


t 


Acceleration  (cm/sec^) 


90 


No  liquefaction  during  the  Elmore  Ranch  event 

Soil  remained  in  elastic  range  during  the  Eimore  Ranch  event 

r. «  0.45  for  equivlent  Cumulative  Horizontal  displacements  in  both  events 


Earthquake  Data  is  Transient,  so 

Fourier  Estimate  of  System  is  Poor. 

(Z\  2DOF,  variance  ~  mean) 

Use  Parametric  Models 

ARMA  (Auto  Regressive-Moving  Average)  model 
Originally  Derived  for  Single-Degree-of- 
Freedom  (SDOF)  System  (Yule,  1927). 


ARMA-based  Parametric  Models  are  - 

Ideal  for  short  data  sequences  (=  50) 

Most  parsimonious  system  estimator 

Statistical  approach  allows  quantification  of  (un)certainty 


Model  System  as  Linear  Filter 


’‘t 

input 


Physical 

System 

yt 

niitniit 

y.  =  b,x,  +  +  •  •  •  +  a  j,_,  +  + 


y  =  actual  system  output  sequence 
X  =  actual  system  input  sequence 
t  =  time-step  counter 


Which  can  be  rewritten  as  - 


nb  na 

/•=0  *=1 


na  =  AR  order 
nb  =  MA  order 


Applying  the  shifting  theorem  and  Z-Transform  Yields 


l-a^z'  -a^z" - 


-  The  Frequency  Domain  Transfer  Function  of  the 

System  of  Interest 


Output  is  a  Combination  of : 

input  history  acted  on  by  "b"  coefficients 
past  outputs  acted  on  by  "a"  coefficients 

AR  =  noncausal  process;  parameters  contain  resonant 
frequency  and  damping  ratio. 

MA  =  causai,  convoiutional  process;  parameters 
contain  the  phase  relations. 


/.  N  -  DOF  System  Modeled  as  2n  -  2n  ARMA 


Solve  For  ARMA  Coefficients  by  Various  Least- 
Square  Based-Schemes 


Earthquake  Data  is  Non-stationary! 

Transient,  non-repetitive  signal. 

Excitation  statistics  (mean  &  variance)  change  through 
time. 

1"^  System  statistics  change  through  time. 

Neither  parametric  nor  Fourier  methods  are  proper  for 
non-stationary  signals. 


Solutions? 


^  Cut  data  into  pseudo-stationary  pieces  (too 
short  for  reliable  Fourier  estimates).  Use 
Berg  parametric  algorithm. 

^  Use  recursive  (e.g.  Extended  Kalman  filter) 
techniques.  Characterize  the  evolution  of 
the  dynamic  system  directly. 


Time  (s) 

Simulated  Elmore  Ranch  velocity  and  the  actual. 
-  Filter  derived  from  Superstition  Hills  data. 
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Surface  Wave  Method  -  A  Tool 

for 

Lifeline  Earthquake  Engineering 

by 

Ronald  D.  Andrus 

Structures  Division 
Building  and  Fire 
Research  Laboratory 
National  Institute  of 
Standards  and  Technology 
Gaithersburg,  MD  20899 


Introduction 


Liquefaction  of  granular  soils  is  one  of  the  major  causes  of  lifeline  damage 
in  past  earthquakes.  Screening  techniques  based  on  geology,  hydrology,  and 
soil  conditions  can  identify  areas  along  the  lifeline  corridor  requiring  more  rig¬ 
orous  analyses.  These  areas  can  extend  for  several  kilometers,  however.  The 
Spectral-Analysis-of-Surface-Waves  (SASW)  method  has  great  potential  for 
rapid  determination  of  the  layer  thickness  and  small-strain  shear  wave  velocity, 
Vj,  of  soil  deposits.  The  SASW  method  is  nonintrusive  and  nondestructive.  It 
is  based  on  the  principal  that  surface  seismic  waves  of  high  frequency  propa¬ 
gate  only  in  near-surface  layers,  and  surface  waves  of  low  frequency  propagate 
at  different  velocities  if  stiffness  varies  with  depth.  Thus,  different  portions  of 
the  soil  profile  can  be  tested  by  using  surface  waves  over  a  wide  range  of  fre¬ 
quencies.  The  general  SASW  test  configuration  of  source,  receivers,  and 
recording  equipment  is  shown  in  Figure  1.  Liquefaction  potential  and  Vj  of 
granular  soils  are  influenced  by  many  of  the  same  factors  (e.g.,  density,  con¬ 
finement,  and  geologic  age).  This  paper  evaluates  the  ability  of  the  SASW 
method  to  delineate  liquefiable  soil  using  data  from  two  sites  (Goddard  Ranch 
and  Larter  Ranch)  where  liquefaction  occurred  during  the  1983  Borah  Peak, 
Idaho,  earthquake  and  the  liquefaction  assessment  procedure  by  Stokoe  et  al. 
(1988b).  Application  of  the  SASW  method  to  lifelines  is  discussed. 


Introduction 


Examples  of  Correlating  Vg  and 
Liquefaction  Potential 


Stokoe  et  al.  (1988b)  applied  the  strain  approach  by  Dobry  and  his  col¬ 
leagues  (1982)  in  an  analytical  study  of  the  liquefaction  potential  of  sandy  soils 
in  the  Imperial  Valley,  California,  to  generate  liquefaction  assessment  charts 
based  on  measured  Vg  and  maximum  horizontal  ground  surface  acceleration, 
a^,  for  stiff  soil  (nonliquefiable  clay)  site  at  the  candidate-site  location.  The 
chart  for  15  cycles  of  shaking,  the  approximate  number  of  cycles  of  the  Borah 
Peak  earthquake  (Mg  =  7.3),  at  a  level  ground  site  with  the  liquefiable  sand  in 
the  upper  12  m  is  shown  in  Figure  2.  Liquefaction  is  predicted  to  occur  right 
of  the  shaded  region.  Within  the  shaded  region,  liquefaction  would  likely 
occur  depending  on  thickness  and  depth  of  the  liquefiable  layer.  Liquefaction 
is  predicted  not  to  occur  left  of  the  shaded  region  because  the  sand  is  too  dense 
to  liquefy.  During  the  Borah  Peak  earthquake,  the  estimated  values  of  a^  at 
Goddard  Ranch  and  Larter  Ranch  are  0.34  g  and  0.6  g,  respectively.  As  illus¬ 
trated  in  Figure  2,  liquefaction  is  predicted  at  Goddard  Ranch  in  sediments  with 
Vg  ^  110  m/s  and  liquefaction  likely  in  sediments  with  llOm/s  <Vs^ 

170  m/s.  Similarly,  liquefaction  is  predicted  at  Larter  Ranch  in  sediments  with 
Vg  ^  165  m/s  and  liquefaction  likely  in  sediments  with  165  m/s  <  Vg  ^ 

280  m/s. 

At  Goddard  Ranch,  liquefaction  occurred  as  evidenced  by  numerous  sand 
boils  in  the  low-lying  areas  (Youd  et  al.  1985).  Shown  in  Figure  3  are  sedi¬ 
ment  layers  beneath  the  gravelly  side  bar  investigated.  Also  shown  are  several 
penetration  profiles  determined  by  the  Cone  Penetration  Test  (CPT)  and  the  - 
Standard  Penetration  Test  (SPT).  Liquefaction  most  likely  occurred  in 
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Examples  of  Correlating  and  Liquefaction  Potential 


Unit  Cl,  a  loose  (low  penetration  resistance)  sandy  gravel  with  less  than  a  few 
percent  fines  (Andrus  1994).  Unit  C2,  a  loose  to  medium  dense  sandy  gravel, 
using  penetration  resistance  is  predicted  to  be  liquefiable  to  marginally  lique¬ 
fiable  material.  Unit  B  is  a  sandy  silt  with  clay,  and  is  nonliquefiable.  Three 
Vj  profiles  determined  by  the  SASW  method  are  shown  in  Figure  4.  Profile 
SA-85  was  determined  in  an  earlier  study  (Stokoe  et  al.  1988a)  before  penetra¬ 
tion  and  borehole  data  were  available.  Profiles  SA-2  and  SA-3  were  deter¬ 
mined  assuming  the  layering  observed  in  penetration  profiles.  Regions  of 
liquefaction  and  liquefaction  likely  have  been  shaded  in  Figure  4  using  values 
of  Vj,  and  are  in  good  agreement  with  the  assessment  based  on  penetration 
resistance.  Soil  type  is  needed  to  correctly  assess  no  liquefaction  for  Unit  B. 


At  Latter  Ranch,  liquefaction  caused  the  steeply  sloping  (about  34  percent) 
distal  end  of  an  alluvial  fan  to  move  laterally  downslope  about  1  m.  Numerous 
sand  boils  erupted  along  the  toe  of  the  fan.  Shown  in  Figure  5  are  sediment 
layers  beneath  the  slide.  Also  shown  are  penetration  profiles  determined  by  the 
CPT,  SPT,  and  Becker  Penetration  Test  (BPT).  Liquefaction  most  likely 
occurred  in  Unit  C,  a  loose  to  medium  dense  sandy  gravel  with  about  7  percent 
fines  (Andrus  1994).  Beneath  the  zone  of  fissures.  Unit  C  is  predicted  margin¬ 
ally  liquefiable  since  it  exhibits  higher  penetration  resistances.  Three  Vj  pro¬ 
files  determined  by  the  SASW  method  are  shown  in  Figure  6.  These  profiles 
were  determined  in  an  earlier  study  (Stokoe  et  al.  1988a)  before  penetration 
and  borehole  data  were  available.  Regions  of  liquefaction  and  liquefaction 
likely  have  been  shaded  in  Figure  6  using  Vs,  and  are  in  good  agreement  with 
the  assessment  based  on  penetration  resistance. 

These  findings  illustrate  the  ability  of  the  SASW  method  to  delineate  lique¬ 
fiable  materials  in  granular  deposits.  Soil  type  is  needed  to  correctly  predict 
behavior  in  deposits  with  layers  of  soft  clay. 


Examples  of  Correlating  and  Liquefaction  Potential 


Application  of  the  SASW 
Method  to  Lifelines 


The  SASW  method  consists  of  making  field  measurements  of  surface  wave 
velocity  at  various  frequencies  and  determining  Vj  profile  through  a  process 
called  inversion.  Several  different  source  and  receiver  spacings  are  required  to 
measure  over  a  wide  range  of  frequencies  due  to  attenuation  and  near-field 
effects.  Although  it  is  preferable  to  conduct  the  test  with  receivers  spaced 
about  a  common  midpoint  (see  Figure  1)  and  source  locations  reversed,  the 
common  source  test  configuration  shown  in  Figure  7a  is  more  appropriate  for 
automation.  Hiltunen  and  Woods  (1990)  demonstrated  that  the  sacrifice  in 
quality  is  small  with  the  common  source  configuration.  Recent  work  by 
Nazarian  et  al.  (1994)  has  shown  that  inversion  can  be  automated,  reducing  the 
time  for  testing  and  determining  Vj  profile  at  a  site  like  the  one  shown  in  Fig¬ 
ure  7a  to  about  30  min.  For  alignments  of  moderate  length,  test  arrays  could 
be  effectively  overlapped,  as  shown  in  Figure  7b,  to  produce  profiles  at  dis¬ 
crete  intervals  of  4x.  Since  most  liquefiable  soil  layers  lie  within  the  upper 
10  m  of  soil  deposits,  a  typical  value  of  x  would  be  2  m  for  profiling  depths  of 
0.5  m  to  10  m.  Moving  source-receiver  systems  are  under  development  at  the 
University  of  Texas  at  Austin  (Stokoe  1995)  for  generating  continuous  Vj  pro¬ 
files  rather  than  profiles  at  discrete  locations.  Based  on  these  developments, 
the  SASW  method  appears  well  suited  for  rapid  profiling  along  lifelines. 
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Application  of  the  SASW  Method  to  Lifelines 
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Shear  Wave  Velocity,  Vs,  m/s 


Figure  2.  Liquefaction  assessment  chart  by  Stokoe  et  al.  (1998b)  based 
on  Vg  and  peak  horizontal  ground  surface  acceleration,  at 
stiff  soil  site  for  1 5  cycles  of  shaking,  Nj. 
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Figure  3.  Cross  section  of  the  gravelly  side  bar  at  Goddard  Ranch  (Andrus  1994) 
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Liquefaction  assessment  based  on  Vs  of  sediments  beneath  the  gravelly  side  bar  at 
Goddard  Ranch 
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Figure  5.  Cross  section  of  the  lateral  spread  at  Latter  Ranch  (Andrus  1994) 
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Figure  7.  SASW  testing  with  common  source  configuration  along  pipeline 
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Introduction 


Anyone  who  has  experienced  the  process  of  site  preparation  and  site  charac¬ 
terization  for  field  experiments  shares  the  knowledge  that  a  substantial  amount 
of  time  and  resources  must  be  devoted  to  this  phase  of  research.  Often  the  cost 
of  site  characterization  exceeds  the  cost  of  other  operations  associated  with  the 
research.  For  these  reasons  atone  there  is  strong  justification  for  development 
of  experimentation  sites  which  have  been  thoroughly  characterized  and  which 
have  been  made  user-friendly  with  on-site  utilities  including;  electricity,  water, 
telephone  and  field  office  facility.  It  is  also  widely  recognized  that  full-scale 
model  or  prototype  tests  are  extremely  valuable  in  geotechnical  engineering, 
and  well  characterized,  large  field  sites  are  necessary  for  these  tests. 

For  the  above  reasons,  it  became  evident  during  the  past  decade  to  many 
geotechnical  engineers  that  easy  access  by  investigators  to  well-characterized 
field  sites  would  be  of  enormous  benefit.  These  sites  could  be  used  to  evaluate 
new  techniques  of  soil  improvement,  new  methods  of  soil  testing,  new  methods 
of  foundation  and  earth  structure  construction,  as  well  as  other  activities  not  yet 
defined.  Completely  documented  field  sites  with  a  public  database  containing 
results  of  previous  tests  and  site  characteristics  could  also  provide  appropriate 
locations  for  the  installation  of  permanent  instrumentation  for  measurement  of 
site  response  during  earthquakes  and  other  mass  movement  phenomena  where 
appropriate. 

Sites  like  this  exist  in  other  countries  like  Canada,  Brazil,  United  Kingdom, 
France,  Italy,  Japan,  and  Norway,  and  have  promoted  cooperation  and 
exchange  of  information  between  public  agencies,  universities  and  private 
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sector  groups  across  international  boundaries  in  a  cost  effective  way.  Three 
examples  are  presented  to  give  a  flavor  of  the  sites  in  other  countries. 

In  Canada,  two  experimentation  sites  have  been  created  to  study  the  charac¬ 
terization  of  sand  for  dynamic  and  static  liquefaction  potential.  One  site  is  at 
the  Syncrude  location  in  Alberta  and  the  other  is  in  the  Frazier  River  Delta  in 
British  Columbia.  Specifically,  experiments  are  underway  to:  (1)  evaluate 
high  quality  undisturbed  sampling  techniques  like  ground  freezing,  (2)  calibrate 
and  verify  in-situ  tests  particularly  with  respect  to  influence  of  high  overburden 
pressures,  and  (3)  study  triggering  mechanisms  for  both  static  and  dynamic 
liquefaction.  These  experiments  are  expected  to  lead  to  the  improvement  of 
design  of  embankments  including  tailings  dams  subject  to  earthquake  shaking. 

In  Brazil,  an  experimentation  site  was  created  in  1989  at  the  Polytechmcal 
School  of  the  University  of  Sao  Paulo  in  cooperation  with  the  Brazilian  Associ¬ 
ation  of  Foundation  Engineering.  Experimentation  at  the  sites  has  been  princi¬ 
pally  on  load  capacity  of  drilled  pier  foundations.  Extensive  laboratory  tests 
have  also  been  performed  on  samples  recovered  from  this  site.  In  addition  to 
the  pile  load  tests,  experiments  have  been  conducted  on  full-scale  high  pressure 
injected  columns,  diaphragm  walls  and  high  capacity  earth  anchors. 

In  the  United  Kingdom  there  are  several  sites  but  we  highlight  only  the  one 
at  Bothkennar,  Scotland.  This  is  a  deep  soft  clay  site  and  was  created  as  a 
result  of  a  need  recognized  in  the  UK  in  the  1980' s  by  the  Science  and  Engi¬ 
neering  Research  Council  (SERC)  for  a  well  characterized  field  site.  The 
Council  purchased  the  Bothkennar  Site  in  1987  to  ensure  permanent  access  and 
promote  broad  usage  of  the  site.  Experimentation  progressed  at  an  impressive 
rate  at  Bothkennar  and  a  seminar  was  organized  by  SERC  in  1992  to  re-focus 
and  plan  fumre  developments  at  the  site.  In  addition  a  special  issue  of  GEO¬ 
TECHNIQUE  in  June  1992  detailed  the  geologic  and  geotechnical  characteris¬ 
tics  of  the  site.  Some  of  the  work  reported  in  that  article  and  performed  since 
then  includes:  full-scale  negative  skin  friction  tests  on  a  pile  group,  installation 
of  and  load  tests  on  flight  angered  piles,  full  scale  embankment  studies,  and 
geosynthetic  reinforcement  of  several  types  of  earth  structures.  Most  work  has 
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been  performed  by  UK  researchers,  but  international  collaboration  has  also 
played  an  important  role. 
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Creation  of  National 
Geotechnical  Experimentation 
Sites 


Two  workshops  in  1988  and  1991  identified  and  quantified  the  need  for  a 
system  of  multiple-user  geotechnical  experimentation  sites  in  the  United  States 
(Benoit  and  de  Alba  1988  and  1991).  As  a  result  of  these  workshops,  the 
National  Geotechnical  Experimentation  Sites  (NGES)  system  was  established  in 
November  1992  and  now  is  in  full  operation. 

The  first  workshop,  “Designated  Sites  for  Geotechnical  Experimentation  in 
the  United  States,”  (Benoit  and  de  Alba  1988)  was  held  in  September  1988  at 
the  University  of  New  Hampshire  to  explore  the  need  for  and  potential  interest 
in  Geotechnical  Experimentation  Sites.  Conclusions  generated  by  this  work¬ 
shop  included  the  following: 

a.  There  is  a  strong  need  for  a  network  of  designated  geotechnical  experi¬ 
mentation  sites. 

b.  There  are  numerous  individually  operated  sites  which  are  not 
coordinated. 

c.  The  USA  lags  behind  many  other  countries  in  developing  multi-user 
sites. 

d.  A  system  of  sites  could  increase  cooperation  among  geotechnical 
researchers  and  lead  to  more  cost-effective  use  of  research  funds. 
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e.  Sites  are  necessary  to  quickly  and  thoroughly  evaluate  emerging 
technology. 


/.  Sites  could  be  used  to  enhance  engineering  education. 

An  initial  list  of  sites  was  created  by  Benoit  and  de  Alba  and  with  it  the 
interest  in  an  identified  site  system  grew.  From  this  list  it  was  determined  that 
a  sub-set  of  candidate  sites  for  a  nationwide  system  could  be  selected  and  to 
accomplish  this  a  second  workshop,  “Selection  and  Management  of  National 
Geotechnical  Experimentation  Sites,”  (Benoit  and  de  Alba  1991)  was  spon¬ 
sored  by  NSF  and  FHWA  in  October  1991  in  Orlando,  FL. 

The  National  Geotechnical  Experimentation  Sites  system  (NGES  system) 
grew  out  of  the  1991  Orlando  Workshop.  The  workshop  established  the  organ¬ 
izational  framework  for  the  system  and  recommended  a  set  of  sites  to  be  desig¬ 
nated  for  improvement.  It  also  set  into  motion  the  preparation  of  a  fimding 
proposal  to  NSF  and  FHWA  which  would  support  the  framework  organization, 
initiate  improvements  at  some  selected  sites,  and  establish  a  system  database. 


The  resulting  system  consists  of  five  sites  which  receive  funding  for  site 
improvements,  database  creation  and  up-dating,  and  site  management;  and 
35  sites  which  are  included  in  the  system  and  are  described  in  the  NGES  Cata¬ 
log  (Benoit  and  de  Alba  1993)  but  which  currently  receive  no  funding. 

The  NGES  sites  are  overseen  by  a  self-sustaining  Management  Board  of 
eight  members  who  direct  the  activities  of  the  System  Manager,  currently, 

R.  D.  Woods  (University  of  Michigan).  The  principal  investigators  for  the 
NSF/FHWA  NGES  Umbrella  Project  are  Professors  J.  Benoit  and  P.  A.  de 
Alba  of  the  University  of  New  Hampshire.  The  data  base  management  and  site 
subcontract  management  is  the  responsibility  of  the  principal  investigators. 

The  current  cataloged  40  NGES  sites  were  selected  from  a  group  of  about 
eighty  sites,  each  of  which  was  presented  by  a  sponsor.  Three  categories  of’ 
sites  were  established  for  the  system:  Level  I  sites  receive  the  highest  level  of 
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funding  for  site  improvements  and  user-friendly  activities,  Level  II  sites 
receive  some  funding,  and  Level  III  sites  are  listed  in  the  NGES  Catalog  and 
are  included  in  the  data  base  but  do  not  receive  any  NSF/FHWA  funds  at  this 
time. 

Funding  provided  to  the  Level  I  and  Level  II  NGES  sites  up  to  now  has 
been  used  to  establish  and  keep  up-to-date  the  databases,  bring  the  sites  up  to  a 
high  level  of  geotechnical  characterization,  to  make  the  sites  user-friendly,  and 
to  provide  continuity  of  management  with  the  attendant  preservation  of  all  valu¬ 
able  data.  All  NGES  sites  are  available  to  any  interested  and  qualified  individ¬ 
ual,  entity,  or  university.  NGES  site  users  are  expected  to  provide  their  own 
funding  for  specific  experiments  and/or  tests  through  traditional  funding 
sources  or  self  funding  in  the  case  of  industrial  users.  NGES  Sites  are  not  in  a 
position  to  provide  funding  for  research  at  their  sites  by  outside  users. 

In  the  future,  site  managers  may  have  to  make  a  charge  for  site  usage  to 
provide  for  maintenance  and  further  upgrading.  This  charge,  if  any,  will  be 
nominal  for  all  users.  Potential  NGES  site  users  are  encouraged  to  contact  the 
Site  Manager  early  in  the  planning  for  any  use  of  that  site. 
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Description  of  Sites 


The  Level  I  sites  are;  Treasure  Island  Naval  Station,  San  Francisco,  Cali¬ 
fornia  and  Texas  A  &  M  University,  Riverside  Campus,  College  Station,  TX. 
The  Level  II  sites  are:  University  of  Houston,  Houston,  TX;  Northwestern 
University,  Evanston,  IL;  and  University  of  Massachusetts,  Amherst,  MA. 
The  Level  III  sites  are  distributed  over  the  entire  continental  United  States. 


The  Treasure  Island  Site  is  located  in  San  Francisco  Bay  and  is  currently 
operated  by  the  U.S.  Navy.  Treasure  Island  an  artificial  island  formed  by 
hydraulic  filling  on  a  shoal  adjacent  to  Yerba  Buena  Island.  The  hydraulic  fill 
varies  across  the  island,  but  it  is  basically  loose,  fine  to  medium,  silty  sand, 
with  occasional  clayey  zones.  Because  of  the  likelihood  of  an  earthquake  in  the 
near  future  and  because  of  the  loose  character  of  the  soils.  Treasure  Island  is  a 
unique  site  for  studying  the  potential  for  liquefaction  of  soil  due  to  earthquake 
excitation.  Figure  1  shows  a  sand  boil  on  Treasure  Island,  indicative  of  lique¬ 
faction  of  the  loose  sand  under  the  surface,  which  occurred  during  the  Loma 
Prieta  Earthquake  in  1989. 


The  Texas  A&M  University  Site  actually  consists  of  two  sub-sites,  one  pre¬ 
dominantly  clay  and  one  predominantly  sand.  Both  are  located  on  the  River¬ 
side  Campus  of  TAMU  in  College  Station,  TX.  These  sites  have  been  set 
aside  for  geotechnical  testing  since  the  late  1970’s.  Many  organizations  both 
within  and  outside  of  the  TAMU  system  have  used  this  site. 

Previous  studies  at  the  clay  sub-site  started  with  tests  on  horizontally  loaded 
drilled  shafts  by  the  Texas  Highway  Department  in  the  late  1970's.  Other  shal¬ 
low  and  deep  foundation  studies  and  in  situ  testing  studies  have  been  performed 
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in  the  intervening  years.  Two  substantial  full-scale  soil-structure  experiments 
have  been  performed  at  the  sand  sub-site.  The  first  consisted  of  a  70  m  long 
by  5  m  high  embankment  passing  over  a  full-size,  instrumented  culvert.  The 
second  consisted  of  an  instrumented,  60  m  long  and  8  m  high  tied-back, 
ground-anchored,  wall,  which  can  be  seen  in  the  background  in  Figure  2. 

This  sub-site  was  also  the  location  of  the  SPREAD  FOOTING  PREDIC¬ 
TION  SYMPOSIUM,  1994,  sponsored  by  the  Federal  Highway  Administration 
and  held  in  conjunction  with  the  ASCE  Geotechnical  Engineering  Division 
Specialty  Conference  on  the  prediction  of  Vertical  and  Horizontal  Deforma¬ 
tions  for  Foundations  and  Embankments  over  sand.  Considerable  new 
characterization  data  were  collected  for  this  settlement  event  in  the  area  shown 
in  Figure  2.  The  Prediction  Symposium  is  described  in  detail  in  Briaud  and 
Gibbens  (1994). 

Previous  characterization  of  the  University  of  Houston  Site,  a  deep,  stiff 
clay  site  has  been  extensive  and  consists  of  many  borings  from  which  samples 
have  been  retrieved  and  subjected  to  various  laboratory  tests.  This  was  the  site 
of  the  full-scale,  nine  (9)  pile  group  load  test  sponsored  by  FHWA  and  per¬ 
formed  by  Prof.  Michael  O’Neill. 


Extensive  site  characterization  of  the  Northwestern  University  Site,  a  lake- 
fill,  loose  sand  site,  was  performed  in  preparation  for  the  1989  ASCE  Geotech¬ 
nical  Congress,  at  which  a  pile  load  prediction  event  was  held.  Four  test  piles 
and  nine  anchor  piles  were  installed  for  this  prediction  event.  Figure  3. 

The  deep,  varved  clay  deposits  at  the  University  of  Massachusetts  Site 
generally  are  in  the  age  range  of  10,000  -  8,000  years  old.  A  considerable 
body  of  field  and  laboratory  testing  data  from  a  variety  of  tools  is  available. 

Thus,  through  NGES,  a  resource  is  available  to  geotechnical  engineers  for 
testing  new  ideas  for  earthworks  construction,  in  situ  measurement  of  earth 
parameters,  and  long  term  measurement  of  ground  behavior.  The  NGES  sites 
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are  available  to  all  interested  users,  government,  university  and  private  sector. 
The  well-characterized  NGES  sites  should  encourage  a  variety  of  field  experi¬ 
ments  which  should  lead  to  better,  safer  and  more  economical  earth  structures. 


Any  interested  engineer  is  encouraged  to  obtain  a  floppy  disk  version  of  the 
NGES  Database  and/or  a  Catalog  of  NGES  Sites  from  Prof.  Jean  Benoit,  Uni¬ 
versity  of  New  Hampshire. 
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Appendix  I  -  Lecture  Slides: 
Effect  of  Stress  State  and  Soil 
Fabric  on  Seismic  Wave 
Velocities  Measured  by 
Borehole  Methods  -  Prof.  Ken 
Stokoe 
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Simplified  Sketch  of  Experimental  Set-Up 
for  Seismic  Wave  Testing  in  the  LSTC 
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b.  Shear  Waves  with  Horizontal  Particle  Motion  (termed  SH-Waves) 
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c.  Shear  Waves  with  Vertical  Particle  Motion  (termed  S  V-Waves) 
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Directions,  Planes,  and  Stress  States  Associated 
with  P-WaveVelocities  Measurements 
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Variation  in  P-Wave  Velocitias  with  Propagation  Direction 
for  Measurements  in  the  Horizontal  and  Vertical 
Planes  under  Isotropic  and  Biaxial  Loadings 
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Directions  of  Ray  and  Wave  Normal  on  the  Wave 
Surface  from  a  Point  Source  in  Isotropic  and 
Anisotropic  Full  Spaces 
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Simplified  Schematic  of  Crosshole  Seismic  Test 
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Three-Component  Velocity  Transducer  Used  in 

Crosshole  Testing 
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Comparison  of  Shear  (SV)  Wave  Velocities  Measured  in  the 
Forward  and  Reverse  Directions  Along  Path  1  at  Treasure  Island 
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Comparison  of  Shear  (SV)  Wave  Velocities  Measured  Along 
Path  1  at  Treasure  Island  with  Theoretical  S-Wave  Velocities 
Predicted  by  Hardin  and  Drnevich  (1972)  for  a  Normally- 
Consolidated  Soil  with  Void  Ratios  (e)  of  0.5  and  1.0 
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Average  Shear  Wave  Velocity  Ratio  (SV/SH)  Along 
Paths  1,  2  and  3  at  Treasure  Island 
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